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Theorie der Druckstabilität der Sandwichplatte I. 
Von H. Neuber, Dresden 


Für die Sandwichplatte (dreischichtige Platte) werden die Grundlagen der elastischen Stabilität in der 
vom Vf. 1943 aufgestellten verschärften Form?) integriert. Die Lösung klärt den Einfluß der verschiedenen 
Parameter auf die Knicklast und bestätigt einwandfrei die Existenz von zwei charakteristischen Instabili- 
Lätsbereichen. 


The basical equations of elastic stability are integrated for sandwich plates in the ewact form stated by 
the author 1943°). T’he solution clears the influence of the different parameters on the buckling load and shows 
the existence of two characteristic regions of buckling. 


Les öquations fondamentales de l’elastique stabilitE sont integrees pour des plaques i Sandwich dans 
leur form exacte, donnee par l’auteur 1943°). La solution montre l’influence des parametres differentes. 
L’existence des deux domaines d’Ecrasement est demonstre£e. 


ÖOcHoBbt yupyrol yCTOüYNUBOCTU CAHABHYHOH INIACTUHBI (TPEeXCIOÜHOH IMIACTHHEI) HHTETPY- 
pylTc# B YTOo4HeHHOii d$opMe°), Amanmoii aBTopoM B 1943 r. Pemerme BbIscHfdeT BimsHHe 
PasımyHBIX NapaMeTpoB Ha YCHAIMe IPONOBHOTO HarHOa U Öe3yCJIOBHO TIOATBEPIKNAET CYINECTBO- 
BAHHE ABYX XAP&KTePHCTHYyecKuX O0JACcTeli HeyCTOHYHBOCTH. 


1. Übersicht 


Der Begriff ‚‚Sandwichplatte‘‘ bezieht sich auf geschichtete Platten (auch Verbundplatten 
genannt); man versteht darunter ebene Platten, die aus mehreren Schichten verschiedener Bau- 
stoffe hergestellt sind. Bei geeigneter Wahl der Abmessungen der einzelnen Schichten, sowie der 
Werkstoffe zeichnen sich solche-Platten durch besonders geringes Gewicht bei verhältnismäßig 
hoher Steifigkeit aus. 

Vorliegender Bericht, der als Fortsetzung meiner früheren Arbeiten über das Verbundstab- 
problem!) entstanden ist, gibt im technisch besonders wichtigen Falle der Druckbeanspruchung 
über die günstigsten Abmessungen solcher Platten Aufschluß. Hierbei macht die naheliegende 
Frage der Konkurrenzfähigkeit mit entsprechenden Baugliedern der Blechbauweise eine besonders 
sorgfältige Behandlung des Stabilitätsproblems erforderlich, um so mehr, als eine über den Gegen- 
stand bereits vorhandene englische Arbeit?) weitgehende Vernachlässigungen enthält, derart, daß 
sich dort u.a. der Grenzübergang zum Euler -Verhalten des Gesamtquerschnittes, der für 
große Knicklängen erfüllt sein muß, für die dreischichtige Platte noch nicht herausstellte. Um 
dieser Forderung zu genügen, mußte eine möglichst strenge Lösung des Problems der Druck- 
stabilität der Verbundplatte angestrebt werden. Wie im vorliegenden Bericht nachgewiesen wird, 
gelingt es im Anschluß an meine Arbeiten über die Grundgleichungen der elastischen Stabilität 
und ihre Integration®), die strenge Lösung des Problems aufzustellen (genaue Formulierung der 
Randbedingung beim Übergang von einer Schicht zur benachbarten, Berücksichtigung der Druck- 
vorspannungen auch in der leichteren Schicht usw.) Die Rechnung wird für die beiden wichtigsten 
Fälle der Zweistoffbauweise, die zweischichtige und die dreischichtige Platte durchgeführt; sie 
liefert jeweils eine Beziehung für den kritischen Wert der Zusammendrückung (Eigenwert) und 
damit für die Knicklast. Die Auswertung bestätigt einwandfrei die Existenz eines für große 
Knicklängen geltenden Eulerbereiches des Gesamtquerschnittes. Zur. Unterscheidung vom Falle 
der selbständigen Knickung oder Beulung der dünnen Außenschicht (‚‚erster Eulerbereich‘‘) wird 
der erwähnte Bereich als ‚‚zweiter Eulerbereich‘ gekennzeichnet. 

Zur bequemen Ermittlung der optimalen Knicklast bzw. der günstigsten Wandstärken 
wurden besondere Diagramme entwickelt, in welchen ein der Knicklast entsprechender dimensions- 
loser Kennwert über der reziproken Länge aufgetragen ist. Dieser Teil des Berichtes geht wesent- 
lich über die bereits in einem früheren Bericht?) gebrachten Ergebnisse hinaus. Eine im Anschluß 
an jenen Bericht erschienene Arbeit von Flügge und Marguerre?) gibt zur Ermittlung 


1) H.Neuber, Jb. dtsch. Luftfahrtforsch. (1941), 8. I 491; ferner: Techn. Berichte 11 (1944), Nr. 2, 
ID 022. 
| 2) Gough-Elam-de Bruyne, J. Roy. Aeron. Soc., April 1940. 
-3) H. Neuber:, Z.angew. Math. Mech. 23 (1943), 8. 321. 
4, H.Neuber, Druckstabilität der Verbundplatte; Vorbericht (LFA) 1944. 
5) Flügge-Marguerre: Untersuchungen und Mitteilungen, Nr. 1360. 
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der optimalen Knicklast ein Näherungsverfahren an. Anläßlich eines Gastvortrages, den Herr 
Marguerre über das vorliegende Problem an der TH. Dresden im September 1951 hielt, teilte 
mir Herr Marguerremit, daß eine neue Veröffentlichung dieses Dimensionierungsverfahrens 
in Frankreich beabsichtigt ist. Jenes Verfahren hat zum Ziel, die Vielzahl der auftretenden Para- 
meter durch Näherungsbetrachlungen unter Ausnutzung funktionaler Zusammenhänge soweit 
zu verringern, daß sich eine für praktische Zwecke geeignete Methode ergibt. Die Grundlage der 
hierbei benutzten Betrachtungsweise bildet jedoch die von mir bereits im Jahre 1944 aufgestellte 
strenge Lösung des Problems?), insbesondere der exakte funktionale Zusammenhang zwischen der 
Stauchung und der sogenannten reduzierten Wandstärke (Verhältnis der z-fachen Wandstärke 
zur Knickhalbwelle). 

Diese Lösung, welche in den Jahren 1945 und 1946 weiter vervollkommnet wurde, wird nun- 
mehr mit der vorliegenden Arbeit veröffentlicht. Dabei sind die verschiedenen Grade der Ideali- 
sierung klar herausgestellt. Bei der 1: Idealisierung handelt es sich um die Ausnutzung der Tat- 
sache, daß die Stauchung wesentlich kleiner als List. Bei der 2. Idealisierung wird die reduzierte 
Wandstärke der Außenhaut kleiner als 1 und kleiner als die reduzierte Gesamtwandstärke vor- 
ausgesetzt. Bei der 3. Idealisierung wird außerdem auch die reduzierte Gesamtwandstärke kleiner 
als 1 angenommen. Die exakten Formeln werden im vorliegenden Bericht innerhalb jedes Ideali- 
sierungsbereiches explizit aufgestellt. Für die zweischichtige und dreischichtige Platte wird der 
Zusammenhang zwischen Stauchung und reduzierter Wandstärke einerseits, sowie zwischen der 
Knicklast und dem Verhältnis von Wandstärke zu Baulänge in Diagrammen wiedergegeben. In 


den Knicklastdiagrammen sind außerdem die Kurven für konstantes Baugewicht eingezeichnet, 


so daß damit bereits eine exakte Optimallösung gefunden werden kann, die für die praktische An- 
wendung von grundlegender Bedeutung ist. 
2. Die Grundgleichungen des Problemes 


Als Ausgangsgleichungen des Druckstabilitätsproblemes dienen die Gl. (50) der Arbeit des 
Verfassers über die Grundgleichungen der elastischen Stabilität?). In kartesischen Koordinaten 
x, y, zohne Massenkraft lauten diese Gleichungen 


N] oK,; o®V, 2 02V ) 
TER NENNE RE! en D En 
> > Ta K,,- rd KlsV. sur N eh: 


r,1=%.Y,2 


Hierin ist K,, der Spannungstensor nach der Deformation und V, der elastische Verschiebungs- 

vektor. Die x-Achse möge entsprechend der auf Bild lersichtlichen Skizze parallel zur Platten- 

ebene und in Richtung der Druckkräfte liegen, 

während die y-Achse senkrecht zur Platten- 
ebene orientiert ist. 

Die ersten Ableitungen des Verschie- 

bungsvektors liefern den Formänderungstensor 


s = an .. (2), 


Ge Are Dr 
welcher durch das Hookesche Gesetz mit 
dem Tensor t,, der bei der Knickdeformation 
geweckten Spannungen verknüpft ist: 


+ 
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Bild 1. Bezeichnungen bei der Bild 2. Bezeichnungen bei a G d, 6, —]): 3 . 
dreischichtigen Platte der zweischichtigen Platte Bi n us e. Er MC 2 ( ) 


Hierbei ist 


N min tnieeniern 
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dem Kroneckerschen Symbol entspricht. 

Für schmale, d.h. stabartige Platten gilt der ebene Spannungszustand mit t,,—= 0, wobei 
die Spannungskomponenten nur von z und y abhängen, wenn von kleinen Zusatzgliedern mit 2? 
abgesehen wird. Aus Gl. (3) folgt hier 


d, = — 


die Volumdehnung, während 
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Wird die Summe der Dehnungen in der Plattenebene mit e* bezeichnet, d.h. 
erde Qu, SA EN ee 2, (2), 


so ergibt sich aus Gl. (6) nach Elimination von d,,: 


e e> s 
—— = ——— 0... ern ne. BB). 
m—2 m—l|l (8) 
Danach lautet das Ho okesche Gesetz für die ebenen Spannungskomponenten 
e* Ä 
= 2 «(a +6; N) N N ET EC) F 


Die entsprechend Gl. (2) an sich noch auftretenden Schubkomponenten t£,, und t,, werden hierbei 
vernachlässigt. 

Für breite Platten, die hier den ieh Gegenstand der Untersuchung bilden, gilt der 
ebene Formänderungszustand, bei welchem vorausgesetzt wird, daß die dritte Verschiebungs- 
komponente verschwindet und die Spannungen im übrigen wieder nur von x und y abhängen. 
Die Schubspannungen ?,, und £,, verschwinden dann exakt, während für die ebenen Spannungs- 
komponenten Gl. (3) mit e* statt e bestehen bleibt: 


* 
20 (d.+ ne für 1T,S—=L,Y u Se a N ee (9a). 


Die beiden Formulierungen des Hookeschen Gesetzes (Gl. (9) bzw. (9a)) können bei Einführung 
einer Größe c entsprechend der Festsetzung 


1 
zn für breite Platten (ebener Formänderungszustand), 
bzw. (10) 
1 . 
c=-—-— für Stäbe und schmale Platten (ebener Spannungszustand), 


m-+1 


in der gemeinsamen Form 


e* 
2 (dt) (11) 


geschrieben werden, die nachstehend weitere Verwendung finden soll. 


In Gl. (1)ist zu beachten, daß K,,den Spannungstensor nach der Deformation kennzeichnet. 
Da beim praktischen Problem stets nur der Vorspannungstensor vor der Deformation bekannt ist, 
welcher hier mit T',, gekennzeichnet werden soll, so muß noch eine Beziehung zur Ermittlung von 
K,, eingeschaltet werden. Der exakte Weg würde®) eine eingehende Untersuchung der Vor- 
deformation erforderlich machen, wobei es auf sehr kleine, kaum meßbare Größen ankäme. Für 
praktische Zwecke ist es zu empfehlen, diese Schwierigkeit durch Einschaltung einer physikalisch 
sinnvollen Annahme zu umgehen, wie es auch in der Literatur (meist unbewußt) geschehen ist. 


Verf. empfiehlt folgende Beziehung 
| K,,= re lee — > (4: Ts: +d,.T,.) N A WR 1 Se (12), 
=t,y,2 


welche auf der Tatsache beruht, daß die mit dem Knickvorgang verknüpfte Deformation für ein, 
die (quasi-starre) Drehung des Materials mitmachendes Koordinatensystem einer Transformation 
entspricht, welche aus den Ableitungen des Verschiebungsvektors hervorgeht, wenn die Drehungen 


oV, N) . (13) 


KemE | Ss ı Or 
unberücksichtigt bleiben, wenn also in der Gleichung 


oV, 
es. Sure (1A) 


nur das erste Glied der rechten Seite Berücksichtigung findet. Dann setzt sich X, , additiv aus dem 
ursprünglichen Vorspannungstensor 7,,, dem bei der elastischen Deformation geweckten Tensor 
t,, und den infolge der Dehnungstransformation hinzukommenden Anteilen von T7,, zusammen, 
und man erhält Gl. (12). 


®) Siehe Anm. 3). 
22* 
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Werden die Komponenten des Vorspannungstensors 7,, als von den Koordinaten unab- 
hängig vorausgesetzt, d.h. 
ME 


Ta. (15) 


so ergibt sich durch Einsetzen von Gl. (12) in Gl. (1) folgende Gleichgewichtsbedingung: 


91,5 09, 009: & 
DD PB 7.)=0 Vu RN EN 
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r,1=%,y,2 
oder 
RE FE | Be ER 
22 vr "29% (Bi r) ae as) 7.) =0 RENNER 


_ Bei dem hier zu untersuchenden ebenen Problem ist nur über x, y zu summieren (Drehungen aus 


der x, y-Ebene heraus bleiben unberücksichtigt) und der Vorspannungstensor T',, hat die einzige 
Komponente 
Tax 2 @ SEE 1 ee Be er SE 


Die Einführung der Größe S ist in dieser Form bei der weiteren Rechnung vorteilhaft. Die Gl.(17) 
gehen so über in 


Oyn | Oz Kal 0? ei FRE (a ak ie) 

2% U | @ S K Y IE u FIAE an TU Y u Y er. a == 

0% 4 yon ay®, er 0% - oYy Rn 072 .1,0%0% a! 
Aus den Gl. (11) ergibt sich 


oV, 08° ) je Ber ) Be oV, ) 
ne 2 2 yuT— 2 z u ’ == —Z 2 2 . 
to: EG er ; 2) a oy 1—2e u ® °y = 0x ) 
Durch Addition der ersten beiden Gleichungen folgt 
2@G 
any = Eee ee Far de A N Se ee EEE (21), 
so daß die Auflösung nach den Dehnungen auf die in anderer Schreibweise bekannten Beziehungen 
oV, oV. 
26 Su ben — EC (tea + byy)> ae — tl) 
führt. Demnach wird 
LE | EL EL ö | 
20 MG age ee ee) 


und entsprechend 


a ea 1? 1 Be Diss 0) 
a ee oe ee 1) NE Se 
Nunmehr lassen sich in den Gl. (19) alle Verschiebungskomponenten eliminieren, und es folgen 
Oz t, Q 
(1 — cS)— + (1— sy er + (1), un = (IM 
0% oy 0% x 
ig 


Obyy Ol 
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Diese Gleichungen stellen die endgültigen Gleichgewichtsbedingungen des Problems dar. 

Zur vollständigen Behandlung des Problems ist noch eine weitere Gleichung erforderlich, 
die sog. Kompatibilitätsbedingung. Sie folgt allein aus den Formänderungsbedingungen und wird 
z. B. erhalten, wenn die erste der Gl. (22) zweimal nach %, die zweite zweimal nach & differenziert 
wird und hiervon die mit 2 multiplizierte und nach x und % differenzierte dritte der Gleichungen 
(20) subtrahiert wird, und lautet 


— (1—0)8 


N 
+ +48) 


ET 021, 

zul 2. DUgET #3 — Bd: 
ag 522 a Alta tty)=0.....% (26). 
Hierbei ist A der Laplacesche Operator in der z, y-Ebene: 
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3. Die Integration der Grundgleichungen 


Die Grundgleichungen (25) lassen sich durch Einführung einer Spannungsfunktion inte- 


 grieren. Werden die Spannungskomponenten nämlich entsprechend folgendem Ansatz aus einer 


Spannungsfunktion hergeleitet, so sind die Grundgleichungen identisch erfüllt: 


i j 0 F a N) 
DS : 8) - _- = \ SE 
We 0er or 
a —e8) + (1—c)$S (1 — 8) ER rs l28): 
02 F 
.,=|—1 — 1-+2e) 8? - 


Durch Einsetzen dieser Ausdrücke in die Gl. (25) ist die Richtigkeit des Ansatzes ohne weiteres 
ersichtlich. 

Andererseits führt die Kompatibilitätsbedingung (26) durch Einsetzen der Ausdrücke für 
die Spannungskomponenten auf die Differentialgleichung der Spannungsfunktion F und damit 
des Problems. Unter Umgehung einer kurzen Zwischenrechnung ergibt sich 


2:02 E Kor 
/ S) — u (2 
Al S) 23 + (1 > ayr. 0 (29) 
oder mit Einführung der Größe 
1+S . . 
er EN RE KURORT RR >10) € 
2 2 MM 
ne An al). 
0 Aley)’ 


Hieraus geht hervor, daß F sich aus zwei Arten harmonischer Funktionen aufbaut, und zwar 


. solchen des &, y-Systems und solchen des z, (oy)-Systems; d.h. es kann gesetzt werden: 


ED ED, (2, Ne een 6 (3ka), 
Ad, y)=0 und Ads )=0 ...... 0.0... (82) 


gilt. Es könnte auch eine Substitution durch komplexe Funktionen vorgenommen werden, 
was hier jedoch nicht lohnt, da schon einfache Funktionen zur Lösung des Problems ausreichen. 
Es genügt nämlich zur Lösung bereits der folgende einfache Ansatz, wobei zu beachten ist, daß 
für jede Schicht der Platte eine besondere Spannungsfunktion definiert werden muß. 


Für die dicke Schicht (Füll- oder Stützschicht): 

F=sin (pa) [A Sin (py)+ B&in (opy)+CCof (Py)-+DEoflepy) - - - : (33). 
Für die Außenschicht („Blech“, Index 1): 
F,= sin (px) [A, Sin (pyı) + Bı Sin (ep Yyı) + Co] (pyı) + Di&ol (apyı)l- - - (34). 


wobei 


Mit Bezug auf Gl. (18) ist 


U a 3 
= —[1+-)e ee ee er (30), 


Hierbei stellt e die für alle Schichten gleich große negative Vordehnung dar. Demnach nimmt auch 
S für jede Schicht einen anderen Wert an, sofern die Poissonsche Konstante einen anderen 
Wert besitzt. Da dieser Einfluß jedoch nicht wesentlich ist und nur die Rechnung erschweren 
würde, soll zur Vereinfachung 1/m = 1/m, und damit c=c,, S=S,, 0= o, angenommen werden. 
Die Koordinate y, wird bei der dreischichtigen Platte (Bild 1) zweckmäßig mit y—h/2 festgelegt, 
während sie bei der zweischichtigen Platte (Bild 2) mit yidentisch ist. Die Funktion sin (px) be- 
deutet eine sinusförmige Durchbiegung der Platte beim Knicken. Wird die Knickhalbwelle mit I 
bezeichnet, so entspricht pl der halben Periode der Sinusfunktion, d.h. es gilt 9=z/l. 


S= 


4. Aufstellung der Randbedingungen 


An der Übergangsstelle (Verbundstelle) von einer Schicht zur benachbarten halten sich die 
an der Grenzfläche angreifenden Spannungskomponenten das Gleichgewicht. Zugleich ist bei 
absolutem Haften der Flächen aneinander eine Übereinstimmung der Verschiebungskomponen- 


Er 
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ten zu erwarten. An den freien Oberflächen verschwinden andererseits die Randspannungen. Es 
gelten daher folgende Randbedingungen: 
l. Ku - Lu, = 
Il. Kyy— Kyy =0, | 
An der Verbundstelle: ER De a Re 
| 1 en ai] | 


Freie Oberfläche der Außenschicht: 


ek 0: N EG NE En. (37). 
Freie Oberfläche der Stützschicht: 
NILTR,Fe05 NIIESK,), = 8 A So 


Bei der dreischichtigen Platte genügt infolge der Symmetrie die Beschränkung auf eine Seite des 
Systems, so daß nur eine Außenschicht betrachtet wird; die Gl. (38) treten nicht auf. Ferner 
sind beim gegensymmetrischen Knickfall die Konstanten (' und D, beim symmetrischen Knick- 
fall Aund Bin Gl. (33) Null zu setzen. Die restlichen 6 Konstanten sind alsdann aus den sechs 
Bedingungen (I bis VI) zu ermitteln. Bei der zweischichtigen Platte sind sämtliche Gleichungen 
zu berücksichtigen und aus ihnen die 8 Konstanten zu bestimmen. Da das System der Bestim- 
mungsgleichungen in jedem Falle homogen ist, existiert nur dann eine von Null verschiedene Lö- 
sung, wehn die Determinante verschwindet. Aus dieser Bedingung folgt mithin jeweils der kri- & 
tische Wert von S bzw. e. 

Der Rechnungsgang besteht zunächst darin, die in den Randbedingungen auftretenden 
Spannungs- und Verschiebungskomponenten auf die Spannungsfunktion zurückzuführen. 

Für XK,, und K,, gilt mit Bezug auf Gl]. (11), (12) und (18) 


Kt - la le) Us SD 

Aus den Gl. (28) folgt dann weiter 
a 
RD: DE SRRRERALN, 


0° F 
KK, = (18)\(1 ee Apr} N Ss (1—8) öy? 
Um die in den Randbedingungen III und IV auftretenden Verschiebungskomponenten gleichfalls 
auf die Spannungsfunktion zurückführen zu können, ist es erforderlich, nach x zu differenzieren 
(dies ist zulässig, da die Randbedingungen für einen konstanten y-Wert gelten). Mit Bezug auf 
die Gl. (22) und (28) folgt für die erste Ableitung von V,;: 


Day, 2.) 8] 0? F 
HH HIST HAN AI 
Entsprechend ergibt sich für die zweite Ableitung von V,: £ 

"Fr 2a[2 (de Me) We) c 


GE —2G 
0% 0y 02 0y (42). 


0% 
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Nunmehr können die Randbedingungen unmittelbar auf die Ansätze (33) und (34) Anwendung 
finden. Bei Unterdrückung gemeinsamer Faktoren ergeben sich mit den Abkürzungen pgh=2ß, 
Yh,=P, und E,/E=G,/@=H folgende Gleichungssysteme: 


A..Für die dreischichtige Platte mit gegensymmetrischer Auslenkung: 


(Gl. I, II, III, IV mit y=h/2, y=0; Gl. V, VImit y‚=h,; Gl. VII und VIII entfallen; 
C=D=0; gegensymmetrischer Knickfall). 


I. ACojß—A,+o(B&ofoß—B,)=0, { H 
II. (1+(1—2c) 82)(—A6& Sin6 +) +A—)(-BEin oß+D,)=0, 
II. (14 (1— 2c) 8?) (AH &in B—C,) + (14+8)(1-+(1— 2c)8)(BH Sinoß—D,)=0, 
V. (1+(1—26)8?) (AH Cofß—A,) + (1—5)(1— (1—2c) 8) o(BH Cofe P—B)=0, 
V. A,60j6,+B,o&ojoßf, +0, ©inf, +D,oSinoß,—0, 
I. (14 (1—20) 82) (A, Sinß, +0, Cofß,)+(1—8?) (B, Sinoß, +D,Cojoß,)=0. | 


(43). 
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C. Für die zweischichtige Platte: 
(GT, II, II, IV mit y=0, ,=0; Gl. V, VI mit „=--A,; Gl. VIL VII mit y=h). 
I. A—A,+o(B—B)=0, 
IL: (1+(1—20)82)(C—0,)+(1—8?)(D—D,)=0, 
II. (1+(1—2c)82)(CH—C,)+(1+8)(1-+(1—2c)S)(DH—D,)=0, 
IV. (1+(1—2c)82)(AH— 4,)+(1—S)(1—(1— 2c)8)o(BH—B,)=0, 
V. 4,&0j/,+ Bıo&ofoß,— O, Sinß,—D,o Einoßı—=0, 
VI. (14+(1—20)82) (—A, Sin ß,+0,C0ff,)+(1—8?)(— B, Sin oß, + D,Cofoß,)=0, 
VII. A&Coj28-+ Bo&oj206-+CESin2ß+DoSin2oß=0, 


(44). 


_ VI. (1-+(1—20)82) (ASin28+ 0Co]28)(1—8?)(B Sin208-+DEof2of)—0 


Für. den symmetrischen Knickfall der dreischichtigen Platte (Fall B) sind die Konstanten A und 
B (statt C und D), Null zu setzen. Das zugehörige Gleichungssystem geht aus (43) hervor, 
wenn Coj£ mit Sin£ß und ECojoß mit Sin oß vertauscht wird. Es genügt daher, die Rechnung 
nur für eine der beiden Knickformen durchzuführen, da das Endergebnis mit Hilfe dieser 
Vertauschungsregel auch für die andere Knickform benutzt werden kann. Bei der zweischichtigen 
Platte tritt eine Unterscheidung verschiedener Knickformen nicht von vornherein auf. 


5. Die Stabilitätsbedingung 
Die weitere Rechnung besteht in der Auflösung der Gleichungssysteme (43) und (44), bzw. 
in der Aufstellung der aus dem Verschwinden der zugehörigen Koeffizientendeterminanten her- 
vorgehenden Eigenwertgleichungen. Hierbei ist die Einführung folgender Abkürzungen von 
Vorteil: 


1— 82 el dB) 
| BEN EREERTTE TER, RE ee, | 
bzw. een (4); 
SH A20)81 , 
1+(1-20)8 RE 


Bei der dreischichtigen Platte, die zunächst behandelt werden soll, ergibt sich nach Di- 
vision der Gleichungen II, III, IV und VI durch (1--(1—2c) 8?) folgende Koeffizientendeter- 
minante: 


(A) (B) es (Bi) (G) (Di) 
I Cojß oßojoß —1 —o 0 0 
II | — ESinß — yE&inoß 0 0 1 m 
III |HSinß H2@—u)Sinoß 0 0 sn ea 20 (äh). 
IV | 760] £ Huo&o)oß —1 —u0 0 0 
V 0 0 Cojß, oßofoß, Sind, oSinoß, 
VI 0 0 Sinß, ySinoß, Lojf, Yyeojop, 


Hierbei sind die Spalten durch die zugehörigen Konstanten gekennzeichnet, während die Reihen 
entsprechend den Gleichungen mit römischen Zahlen versehen sind. Da die Auflösung der De- 
terminante auf einen Ausdruck führt, der X bis zur zweiten Potenz enthält, ist es vorteilhaft, von 
vornherein die Berücksichtigung des Parameters HZ dadurch zu erledigen, daß eine Zerlegung in 
die drei Fälle H=0, H=1 und H= © vorgenommen wird. Wird der Wert der Determinante 
mit Ö bezeichnet, so erfüllt der Ansatz 


DEN EH)d, THSLH(H 19. 2: 20.202002... (47) 
die Bedingung, daß H nur bis zur zweiten Potenz auftritt; zugleich wird 
” | ph 
(d) a0 do» Br dir); a) SER Ru 
H=o» 


Der Vorteil dieser Zerlegung besteht in der einfachen Berechnung von d,, d, und Ö,, da der Para- 
meter H hierbei nicht mehr stört. 


[% 
a 
v 
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\c Be Für H=0 folgt ein Zerfall der Determinante in ein Produkt von zwei Determinanten: r S 1. 
e De de RR DEE RN Re Br AN . 
2 wobei EN TE RE 
AR (4) (B) | | | | Ke 


" 
se 2 


Co)% oCojoß 


& 
3 
| 
u} 


er — 0 Sind Cof of —yEofß Sinoß.. . . . (50) 

Rx Il —-Sin? —ySinoß * 

‚ und | | | | E 
R: EEE na ann E 

ER a) ) 2 RER Ne | : 

‘ IV es ER u ) 0 : RU 
Be: ‚dr = he er en (De 

vEH V|[ &jß, okojop, Sinf, oESinoß, a“ 5 


YVrl’Sind yo Bor EP 


wird. Zur Berechnung von 9,, wird zweckmäßig folgende Umformung vorgenommen, bei welcher 


En ar der Wert von d,, unverändert bleibt: | 
Ei B)—uea)= BD), (DL Id-DN). 2.2... 
En; Wie erwähnt, sind die Spalten durch die zugehörigen Koeffizienten gekennzeichnet, so daß z.B. RE 
En:  (4,) die Spalte bedeutet, die zu den Koeffizienten von A, gehört. (B,) und (D,)gehendurchden 
Br _ Rechnungsgang (52) in die neuen Spalten (BT) und (D}) über, und die Determinante Oy redu- 
BU; aa sich infolge der auftretenden Nullen auf die zweireihige Form: a 
(Bi) (Di) 
v o&ofoß, — uobojfı oSingfi + (u—2) Sinp, 
vg = RER N RENT. N 
VII ySingf, —uoSinp;  YEofoh: ” a a28ofß, so ee 


Die Ausrechnung ergibt: N R ER 
Pu—ev—on (u—2)+ (eu—Yy (u —2)) Sins, Eine te ny-t HM — Yoorkıeg oh 0. Are 
Für H=1 ergibt sich aus Gl. (46) nach Anwendung der Rechnungsgänge: N var BE 

(4) -+(0,) Sinß = (4m), rc (D,) Sieh — — BI: Eine = : (55) & = Bir 


die vierreihige Determinante: 


Erde ua. 0Cofoß en) e- rs 
IV|  Cojß > weCojoß ’ De % ee 


— (u =[ Vz 2). BT ER Fk, 
V| Sin? ink, oSin oß ein eh of A: ne 


vil Einf &ofa, a, | ee 2 
Nach den Rechnungsgängen: hr A WAHRE 
(4*) + (41) off = (4**), (B*) + (B,) Cof je (B), 


kann auf zwei Reihen reduziert werden, Ana es folgt. x 


leidet) 
el) laty 2). | = 


SR bzw. nach Ausrechnung 


Aal Se y)[e Sin (B Hay & e D + 


3" d=e 


EyaRıy: 


AT E Pe DV & e. 2 ie } ge uf £ 1 2 fl 2 t & "ey 
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“ 


Für H= ist Gl. dem Arsch durch A? zu dividieren, indem die Spalten (A) und (B)j je durch 


H dividiert werden. Dann folgt mit H = » ein Zerfall in eine zweireihige und eine vierreihige 
Determinante: 


en : Bd late TS ea RE PIE TO LS (60), 
_ wobei 
(4) (B) Pier 
III) Sin? (2— u) Sinoß 
du ee = uoSinßbofof + (u— 2)Coff&ineß . . . (61) 
IV | &oj£ uo&ojoß 
und 
(Ah) (B}) (C)) (D,) | 
II—1..—e 0 0 
E20 ) 1 y 
du= RN; BAER a9 tu, (02) 
V| 6&off, oßojoß, Sind, oSinoß, 
vil Sins yEinoh Cofhı  YCofeh, 
wird. 


Mit Hilfe der Operationen 
: (B)— e(4,)=(B}), D)—ylC)=(D) - 02.0.0. (63) 


reduziert sich d,, auf die zweireihige Form 
| | (Bi) ar) 
| V| eGofon—ecein,  eSinen ven] E 
vı ySingß,—eSinp, yEojoh,— y&ojß, | | 
2 und erhält den Wert 
SE = 20y— 2opCofhCofoht(e ty) EinASineh --:--- (65), 


\ Im Falle der zweis chichtigen Platte folgt aus dem Gleichungssystem (44) die acht- 


Dir reihige Koeffizientendeterminante: 


om Bon am. m m. m | 
.o i a 0 e ; 0 Eye Be 0 0 | 
0 1 FR, y ae Zu —_y 
on 0.0 a 


090. Ch obojeh, — Einh, —eEjneh, 
0.0 — End —yEineb Cofhı YCofoh, 
Ina | gSin2gß Bor. 0 
| | 0 0 Basen (66); 
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wobei 
(A) (B) (0) (D) | 
I 1 [ 0 0 
II 0 0 1 Y 
TEN TER, EB RER ae Pe Po (68) 
VI | @0j2#  oCoj20ß Sin2ß oSin2oß 
VII | Sin28# vSin2oß Coj2ß YyEoj2oß 
und 
(A) (Bı) (Ch) (D,) 
Il 0 0 —l u—2 
IV —l — u0 0 0 
Iu= u ee (69) 
V| 6jß, oßofoß, — Sind, —oE©inoß, 
VI | — Sinf, — ySinoß, Colß, yEojopı 
wird. Für %,, führt der Rechnungsgang 
(B)—o(A)=(B*), (DY-pO)=l DIE De re (70). 
auf die zweireihige Form 
(B>) (D*) 
VII] o&of20ß— o&oj2ß oSin2oß— ySin2ß 7) 
en RE 1 REBEL ee ae a ee : 
VII|ySin2oß—o&in2ß YEof2oß— yEoj2% | 
während d,, durch die Operation 
(B)) — uo(A,)= (B}), (Di (Be Ce Dre (72) 
zweireihig wird: (BF) (DY) 
VI o&ojoA—uoboff, —eSinoß,+l2—u) Sind, 3) 
In EEE ET. \ 
VIl —y&inoß,+uoSind, YEojoß+(u—2)Coiß, 
Die Ausrechnung liefert: 
do = 20Y (— 14 E02 ßE&o|20P)— (0? + y?) Sin2ß a ER IT CL SR 


In ep + 2u— u) te (u—2—up)Eofß,Cofeß + (ent 


Für 7=1 kann man durch die Rechnungsgänge 


(A)+(A)= (A*), 


(B)+(Bı)=(B*), 


(HC) (CH), 


(D)+(D)=(D*) (76) 


gleichfalls einen Zerfall in zwei vierreihige Determinanten erreichen. 


Es wird 
Be RR A (77) 
mit 
(4*) (B*) (C*) (D*) 

V| &jß, oßojoß, — Einf, —o&inoß, 
VIl—-ESinf, —yESinoß, off, YEojoß, (78) 

vl &0f2ß 0Coj208 Sin2%  oSin2eß 

vın| &in2g y&in2oß Koj2ß YyLoj2oß 


| 
| 
4 


= 
= 
[ei 
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und 
A ee 
I —1 —o 0 0 
II 0 0 1 —y 
De N 2 er (79) 
III 0 0 —1 u—2 
IV —1 — u0 0 0 


d,, kann ae die Umformungen 


(Be)—o ehr 40)— (Ber), (C*)+Tgbı : (A*)=(C**), (D*)+Xgoß, (B*)=(D**) (80) 


Col, 
auf die dreireihige Form 
Be) (er) (D*) 
— P6&inoß, 
VI 1 y 
+o&ojop,Tgß, Sofhı Cofoh, 
VII o&oj2oß Sin2ß oSin2gß 
I = Cofß, ee Eoj2# +Tgp,Coj2f  -+oTgeß,Coj20 en 
1 
yE&in2oß 60] 2£ yEoj 20% 
VIII Eofoß 
N Gh Sin2ß +Tgb,Sin2f +yTgop, Ein2oß 
übergeführt werden, welche sich durch den weiteren Rechnungsgang 
(B**) + (D**) Sinoß, Cof of — (C**)o Sinf,&of oß, = (B***), 
& . (82 
Dr) pn (Or) = (Die) | = 
in die zweireihige Form 
(D***) % EBAE®) 
ee x » 
Hi Cofoh, Sine@e+P) | eCojop,&oje(2f+%ı) 
— Go op, OR HA)| —eCof oBıCof2P+R) 
1 
De ee (83) 
Se ß, Cojo(28 +8) | yEojop, Sine (28 +6ı) 
VII y | 
7 Cofoß, Sof (2&-+ 1) | — oo] of, Sin(2£ + ß,) 
umwandeln läßt. 
Nach Ausrechnung folgt 
= 20p (— 14 ECof(28+1)Eojo(28+4,)) — (o?+ Y*)Sin(2P + Pı)Sing (2 +P,) (84). 
Für 9,, ist folgende Umformung zweckmäßig: 
(Bı)— e(4)=(B}), D)—ylC)=(D). -:...... (85). 
Dann geht GI. (79) in die zweireihige Form 
(B}) (D,) 
IE 0 a N Pe RE En (86) 
9 


IV] ew-—1) 0 
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h mit der Lösung SH an, | 
E Br 9 et Ve yU)TFN He on BR Sun 
Er über. | EN 
a Im Falle H= © werden in Gl. (66) zunächst die Reihen III und IV durch A dividiert; FR 
Be) dann folgt für H= w wieder ein Zerfall in zwei vierreihige Determinanten: | 
Sa L » D) i e 
= 0 m) de Dar er Ne ee (88), 
er wobei 
B: (4) (B) 2: (0) . (D) 
k BET TREND ER 77% 
” IL) 0 1 2—u | | 
h e IV 1 10 0 0 ER 
Bi! Go Eh | FF TR ee (89) : 5 
Br. VII] &052£ o&oj 20% Sin2£ o&in2oß 
R vın]| Sin 2£ ySin2oß Cof2ß  YyEof2oß a 
iv. und i x 
(A,) (Bı) (0,) (Dı) A 
| -1 0.8 0 0 es: 
11 0 0 —1 —y BES 
— Sal el) 
‚Y]I 680%, okojeh —- Einf —e Sin oßı rd 
wird. : AN re RR 
Zur Vereinfachung von %,, dient die Umformung aan ie 
(B—ne(A)= (Br),  (D+W-IO-DN) : . 2...) 
N, welche auf die zweireihige Form ne ur 
(B*) ze (D*) alas Rt 
VII] 0&oj2 08 —uo Co] 2% oSin2oß+w—2)Ein2B| .:::M) 
vn] yEin2oP — ug Sin2ß Bar er u PER SS Are a 
führt. Die Ausrechnung liefert Be WEL 
| I = ep +2u—n)+E (u —2— up) Cof 26 Sof? ed. (e* Ba y) Sin2ß Einzen ws. 
Für Gl. (90) ist die Umformung | a U er 
(B)— e(A)= (B}), By —y (0 - = 0) 
vorteilhaft. Dann folgt die zweireihige Form: 
(Bi) | (Di) 


v|o&ofeh + 0Cofh —e &imeh ty ein, 5 


Di) Fe 
4 VI| y Sin of, — eo Sin A [v Gof eb, — ar r 


mit der Lösung 


d, = 2oy(—1+ Coj Pf, &of oß) — (ei +w) Sin 
Damit sind alle auftretenden Größen bestimmt. Man erke 

bis auf das Vorzeichen bei der zwei- und dreischichtigen PI, tie 

durch Einführung des Kr ABU a RS Zelane) in. 


= 


Cof 2 sie ge 
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Der zweite Faktor stellt dieselbe Funktion dar, die bereits beim gegensymmetrischen Knickfall 
der dreischnittigen Platte in Ö,, und Ö, auftrat, während der erste Faktor mit der entsprechenden 
Funktion bei symmetrischer Knickverformung identisch ist (Vertauschung von &in mit Co)). 
Dadurch lassen sich die im Endergebnis vorkommenden Funktionen auf die folgenden fünf 
Basisfunktionen reduzieren: 


D,(P) = Sin Cof of — co ß Sin oß, 


®,(ß) = Cof ß Sin oß— - Sin ß Col oß, 


D,(ß) = (2 — u) Col B Sin oB — uo Sin ßCof of, (98). 
D,(b) = 2 — u) Sin P Co) of — wo Co) P Ein 0%, 
d,(ß) = (# (u — 2) — u 0) Sin2ß Sin20oß + u? —2Au—y 
+(u (a — 1) + 2) Co] 2 BCof 2 06 
Nunmehr läßt sich das Stabilitätskriterium auf übersichtliche Weise darstellen. Wird noch bei 


der dreischichtigen Platte durch (—o?), bei der zweischichtigen durch 4 0% dividiert, so kann das 
Endergebnis in folgender Form geschrieben werden: 


Stabilitätskriterium der druckbeanspruehten Verbundplatte 


A. Dreischichtige Platte mit gegensymmetrischer Verformung: 


MB) + Hu u+r— Bl + A) 
(99). 
Pı Bi } 

+H(H—1)40,(% "o.()o.n=0| 

B. Dreischichtige Platte mit symmetrischer Verformung: 

a9, 0,2) + Hut) +B) 
Ah, | (100). 

+ HH—1)10,(, ')®.(5 .)®. o=0| 


» 


C. Zweischichtige Platte: 


A-Dur ale ++) | 
\ (101). 
+ Hana, (2) = 


Wie erwähnt, entsteht Fall Baus Fall A durch Vertauschung von ©in £ mit Co] £ und Sin oß 
mit &oj oß. Man erkennt leicht aus der Darstellung, daß in den drei Extremfällen 7=0; 1; oo stets 
die gleichen Funktionen nämlich ®, und ®,, das Endergebnis beherrschen: Für 7= 0 (Stützstoff 
trägt allein) mit dem Argument /, für A=1 (Platte aus einheitlichem Werkstoff) mit dem 
Argument (#-+f,) bzw. (ß-+ß,/2), für A =  (Außenschicht trägt allein) mit dem Argument 
ß,1/2. Andererseits ist ersichtlich, wie die Stabilitätsbedingung der zweischichtigen Platte bei Ver- 


ß 
tauschung der Schichten, d.h. Vertauschung von H mit und £ mit 5 in sich selbst über- 


Fr 
geht, was schon anschauungsmäßig zu fordern ist. Auch diese Tatsache zeigt — ebenso wie andere 
Kontrollen —, daß die gewonnene Lösung in der Tat für beliebige Schichtdicken gilt. 


Teil II folgt im nächsten Heft dieser Zeitschrift. 
Eingegangen am 20. November 1951. 
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Das Problem der Jungfernquelle ” 
Von H. Wendt in Bonn 


Es sei ein unendlich ausgebreitetes und unendlich tiefes Meer gegeben, das unter der Wirkung der Schwer- 
kraft sieht. Unterhalb der Meeresoberfläche befindet sich im endlichen Abstande h>0 eine räumliche 
Quelle, die nach allen Richtungen gleichmäßig stationär pro Sekunde die Wassermenge e erzeugt. Im Un- 
endlichen bleibt dabei das Meer in Ruhe. Es wird nach der Gestalt der Meeresoberfläche gefragt, wenn die 
Strömung stationär, inkompressibel, reibungs- und wirbelfrei vorausgesetzt wird. 


The gravity is acting upon an ocean of infinite width and length, and a spherical source is situated at 
the finite distance h>0 below the surface, steadily producing water of the amount e per second. At infinity, 
the ocean is still, und the liquid is incompressible and free of viscosity. Questioned is the form of ihe surface, 
if the flow is stationary and irrotational. 

Supposons une mer eiendue & l’infini et de profondeur infinie, ewposee & l’influence de la pesanteur! 
Au-dessous de la surface de la mer & la distance limitee h>0, se trouve dans l’espace une source, qui 
prodwit la quantite d’eau en toutes directions egalement et stationairement par seconde. Dans l’infini la mer 
resie en repos pendant ce procede. Il est question de la forme de la surface de la mer au cas que le courant 
est suppose d’etre stalionnaire, incompressible, sans friction et sans tourbillons. 


Ilycts umeem ÖdecRoHeyHoO IPocTmpamıneechH U ÖeCKOHeYHO TILYÖOKOE MOPpe, HAXOoNALMeecH 
ION Melicrguem Cut TAKecTu. llon Mopckoü HOBePXHOCTB® 3aAaH HA KOHeYHOM PACCTOAHHH 
h > 0 NPOCTpaHcTBeHHkIä HCTOYHHK, IPON3BONAINHUH PABHOMEPHO IO BCEM HAUPAB/IeHHAM H CTA- 
IIMOHAPHO B CeKYHAY KOJIHYEeCTBO BoNbl. Ilpy 3ToM, Ha ÖecKOHeYHOCTH Mope OCTAeTcH B HOKOE. 
CnpammsaertcH. KAKOBa hurypa MOPCKOü IMOBePXHOCTH, ecJIH TeyeHMe IIPeAMoraeTch, CTAIHO- 
HAPHBIM, HeC5KHMAeMBIM U CBOÖOAHBIM OT TPeHuA U BHXpeü. 


1. Mathematische Formulierung 


Das oben formulierte Problem ist rotationssymmetrisch. ‘Wir brauchen die Strömung 
nur in einem Vertikalschnitt, der die Quelle enthält, zu betrachten. Diese Ebene beziehen wir 
auf ein rechtwinkliges kartesisches rz-Koordinatensystem. Sehen wir einen Augenblick von der 
Quelle ab (Ergiebigkeit e = 0), so haben wir ein ruhendes Meer vor uns, dessen Oberfläche mit 
der r-Achse des Koordinatensystems zusammenfallen möge. Hat die Quelle die Ergiebigkeit e + 0, 
so soll die 2-Koordinate der Meeresoberfläche für r— + ® nach Null streben. Die z-Achse 
weise aus dem Meere heraus. Ohne die Fragestellung einzuschränken, können wir die Quelle 
der Ergiebigkeit e im Punkte Q (r = 0, z—= —1) der negativen z-Achse anbringen (Bild 1). 

z Die Geschwindigkeitskomponenten eines Wasser- 
teilchens parallel zur r- bzw. z-Achse bezeichnen wir 
mit v=u(r,2) bzw. w= w(r, 2). wund w lassen sich 
nun sowohl mittelsdes’'Geschwindigkeitspotentials a(r, 2) 
als auch durch die Stromfunktion y (r, z) darstellen: 


1 1 
ul, SEN Vz» w(r, )=-n= TH 3) (1). 


Eliminiert man 9 aus diesen Gleichungen, so gewinnt 
man die Differentialgleichung für die Stromfunktion 


1 OL 
um ler) 0: @ 
oder 


1 
Ye m rt PTR (2a). 


Wir vermerken noch die Bedingung, daß das Meer im 
Unendlichen ruht: 


u=w=(0 (gültigim Unendlichen) . . (3). 


" Bild 1. r,2-Koordinatensysten 


Die Meeresoberfläche ist zugleich eine Stromfläche und eine freie Oberfläche. Auf ihr hat 
also die Stromfunktion y(r, 2) in allen Punkten einen festen Wert, den wir gleich Null setzen. 
Bezeichnen wir die Ordinaten der Punkte der Meeresoberfläche in der rz-Ebene mit 27, so lautet 


ihre Gleichung: 
Br, Zu) = 0. ae Re Pe 
*) Die Arbeit ist angeregt durch eine Naturerscheinung in der Sächsischen Schweiz, auf die mich Herr 


C. Weber aufmerksam machte. j 
!) Wir bezeichnen gelegentlich partielle Ableitungen nach r und z durch Indizes. 
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Von dieser Gleichung selzen wir voraus, daß sie nach 27, auflösbar ist: 
zu=/(r) (Gleichung der Meeresoberfläche) . . .. . . . . . (da). 


Weiter ist auf der freien Meeresoberfläche der Druck p= 9 (r, z) in allen Punkten gleich dem 
Drucke p, = const der aufihr lastenden Atmosphäre. Der Druck läßtsich durch dieBernoul- 
lische Gleichung bestimmen. Bei Anwesenheit einer äußeren Kraft $ (r,z) (pro Volumen- 
einheit), die ein Potential P (r, z) besitzt: 


= — grad P 
hat diese die Gestalt: 


HS (w ae LA RR LAS RE APR 


Or 
— 


Mit o wurde die Dichte bezeichnet. ound Esind Konstante. 

Bei unserem Problem wirkt auf alle Wasserteilchen die Schwerkraft in Richtung der nega- 
tiven z-Achse. Das Potential lautet also P = o gz, wo g = 981 cm sec”? die Erdbeschleunigung 
bedeutet. Damit wird aus (5) die Gleichung 


Pp+5 @+w)+oge=E BI N AR rn CO), 


Zur Ermittlung der Konstanten E gehen wir auf der freien Meeresoberfläche (p = p,) ins Un- 
endliche. Wegen 


u(r, 2u) = 0, w(r, 2y) = 0, 21 = 0, p= 9, (gültig im Unendlichen der Meeresoberfläche) (7) 
findet man sofort: E=m,. 
Damit hat man als Gleichung, der alle Punkte der freien Meeresoberfläche genügen müssen: 


[u(r, 2u)]? + [w (r, zu)? = — 29 2, (gültig auf der Meeresoberfläche). . . . (8). 


Da die linke Seite von (8) eine Summe von Quadraten ist, so folgt, daß die Punkte z > 0 außer- 
halb der Meeresströmung liegen. Insbesondere befindet sich also der Punkt Q,(r=0,2=1) 
außerhalb des Meeres. 

Weiter lehrt uns die Gleichung (8), daß der Punkt r = 0, z= 0 der Meeresoberfläche an- 
gehört. Aus der Symmetrie der Strömung folgt nämlich, daß die innerhalb des Meeres liegende 
z-Achse eine Stromlinie ist. Schreiten wir vom Quellpunkt Qin Richtung wachsender z-Werte 
auf ihr aufwärts, so durchstoßen wir die Meeresoberfläche in einem Staupunkt (v= w=0). In 
diesem Punkte muß nach (8) notwendig z = Osein. 

Mittels der Gleichungen (1) bis (8) können wir das am Anfange der Arbeit gestellte 
Problem mathematisch formulieren. 

Die durch die Quelle in @ erzeugte Meeresströmung ist durch eine Stromfunktion y(r, 2) 
mit folgenden Eigenschaften gegeben: 

- I. Innerhalb des Meeres mit Ausnahme des Punktes @ (0, — 1) genügt % (r, z) der Differential- 
gleichung (2). 
II. Im Punkte Q zeigt y (r, z) die Singularität der Quelle der Ergiebigkeit e. 
III. % (r, z) erfüllt im Unendlichen die Bedingung (3), d.h. das Meer ruht im Unendlichen. 
IV. Die Meeresoberfläche 27, = f(r) hat für r > + die r-Achse zur Asymptote. 
V. Die Meeresoberfläche ist eine Stromfläche. Ihre Punkte erfüllen die Gleichung (4). 
VI. Die Meeresoberfläche ist eine freie Oberfläche. Ihre Punkte erfüllen de Bernoullische 
Gleichung (8). 


2. Die Funktionen 29, 2), on 

Ehe wir das Lösungsverfahren für das Problem der Jungfernquelle darstellen, wollen wir 
uns einen ersten Überblick über die Funktionen verschaffen, die im folgenden immer wieder 
auftreten. 

a) Die Funktion %g (r, 2): 
Wir betrachten zunächst die Funktion 

1+2z 
Nat, 2) = Ge A Sn ET SR (9). 

Die im Nenner und auch in sonstigen Ausdrücken auftretenden Wurzeln sollen stets positiv 
genommen werden. (4(r, z) genügt der Differentialgleichung (2). Es ist die Stromfunktion 
einer im Punkte @ (0, — 1) angebrachten Quelle der Ergiebigkeit 4, die stationär nach allen 
Richtungen gleichmäßig sprudelt. Im Unendlichen induziert die Quelle keine Geschwindigkeiten. 


ur 
2 gr 
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b) Die Funktionen Q®%(r,2), n=0,1,2,-:- 
Ebenso ist 
29 (y, 2) = 1—z 

[r+ A — ze 
eine Lösung von (2). Es ist die Stromfunktion einer Quelle der Ergiebigkeit 4x im Punkte 
Q,(0, +1). 29 (r, 2) geht aus 24(r, z) offenbar dadurch hervor, daß man z in —z verwandelt. 
Man liest aus (2) ab, daß mit 2%(r, z) auch alle partiellen Differentialquotienten nach z Lösungen 
von (2) sind. Wir schreiben sie 


et A 


on AO (r, 2) 


Om (2) = aan nr 2 ee 
Es sind also z.B. 
DO & DT srl —2) 
2 (7,2) PrA—zepRr 23 (1,2) = — PL (A —z)>pR 7,112) 
1m, | Lösungen von (2). 


Die Funktionen 2® (r, 2) 


n-5 (n =.0,7122,7220,) Bstellenydie 
3 Stromfunktionen mehrfacher 
Quellen in Q, dar, die das Un- 
endliche in Ruhe lassen. Sie 
weisen innerhalb der Meeres- & 


strömung keine Singularität auf. 


Zu einer Rekursionsformel 
für die Q®(r, 2) gelangt man, 
wenn man von der quadrierten 
Formel (12) für 2% ausgeht: 


[+ 3] 190, ger 
und beide Seiten partiell nach 
z differentiiert: 
—3(1—2) 2% (r, 2) 

+[r? +(1—2)2] 2% (r,2)=0. 
Bildet man von dieser Gleichung 


die n-te partielle Ableitung 
nach 2 nach der bekannten 


"B 


Y 
{ 


\ 
I 
| 


| 


. 
N Rt 
= 


Leibnizschen Regel, so ent- | 
steht die Rekursionsformel: j 
y | [r? + (1 Ns 2)?] Kolu; ze, 2) } c i 


— (2n+3)(1 — 2) 

2R+D (7,2)+n (n+2), (13). 

2, 2)=0 
BED 


Mit 2%r, z)und QW(r, 2) lassen 
sich damit die Q®(r, 2) für 
5 n—=2,3,4,... successive leicht 
berechnen. 


» 1 
BER r N 
| 2 


c) Die Funktionen o,(r?): 


Für 2=0 bezeichnen wir ab- 
gekürzt 2%(r, 0) mit o,„(r?): 


Bild 2. Die Funktionen o,(%) n=0,1,2,3,4,'5) 


9,(r?) = 2W(r, 0) = 02 or 
Die &, sind offenbar Funktionen von r?. Man hat z.B.: | TE 
1 2 A RE 
N VEN TR BUN ” ) 3r n ne 
Se) ta a) = Torre al) = er ae (iD) Bee 


ar er a a A u Tal Dee 4 
Se ES > r DR, | 
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Setzen wir in (13) für zden Wert Null, so erhalten wir die Rekursionsformel der @,„(r?): 
(+7?) @n+2 (7?) — (2 + 3) @n+1 (7?) + n (n+ 2) w, (r?) = 0 n=0(,1,2,... (16). 
Mit @,(r?) und o,(r?) sind damit die Funktionen o,(r?) (n=2,3,.. .) leicht berechenbar. Auf 


- Bild 2 ist der Verlauf der w,(r?) für n=0,1,2,3, 4,5 dargestellt. 


Aus (15) in Verbindung mit der Rekursionsformel (16) folgen sofort: 
Hilfssatz 1: lim @,(r?) = 0 At Fa BEP ÄRREREEN 
r>+% 
Hilfssatz 2: Die o,(r?)(n=1,2,.,.) gehen mit r — 0 wie r? gegen Null. 
Unter Verwendung der @,(r?) gelten folgende Entwicklungen nach steigenden Potenzen 


von z [| <Y1-+r2] 2): 


zZ (r2 2 
2 r, 2) = 2, ee ) = wur?) + @ı(r?) z en 2: | 
co" Är2 : y2 SEN (17a) 
ar, 2) er ER ( en zi = ol?) — al) z+ 2 1 ER | 
i=0 h N 
00) a) I" @3:-1(r?) TEN s N) ) 
2alr,2)- Q%Xr, z) Paare = 2) 2+ en 2n..d 32. 2.02 (175) 
29 (r, 2) - SH N ee. 
i=0 2 
(17e). 


= or?) + 0541 (?)a+ 
032. 


3. Das Näherungsverfahren 
Das zu Beginn der Arbeit formulierte Problem der Jungfernquelle soll für kleine 
Quellergiebigkeiten behandelt werden. Seine Lösung soll nicht mit der Hodographenmethode 
vorgenommen werden. Es wird ein Singularitätenverfahren entwickelt, das innerhalb der 
Strömungsebene arbeitet. Es gestattet, Näherungen beliebig hoher Ordnungen zu berechnen. 
Wir setzen dabei voraus, daß die Stromfunktion y = y (r, 2) oder genauer y = y (r, z, e) sowohl 


eine Entwicklung nach wachsenden Potenzen von z als auch nach der Ergiebigkeit eloder <= ae 


der Quelle gestattet. Von den auftretenden unendlichen Reihen wird angenommen, daß man 
mit ihnen wie mit Polynomen rechnen darf. Auf Konvergenzfragen soll nicht eingegangen werden. 

a) Die nullte Näherung: 

Um einen ersten Schritt für das Näherungsverfahren zu gewinnen, setzen wir die Ergiebig- 
keit e der Quelle in Q gleich Null. Dann ist 

i vr N. N er LEN 
die nullte Näherung für die Stromfunktion. Die nullte Näherung für die freie Meeresoberfläche 
hat dann die Gleichung: 

MR Sn ne eek: (19), 

b) Die erste Näherung: 

Wir bringen nun im Punkte Q (0, —1) eine Quelle der Ergiebigkeit e-#0 an. Damit die 
Meeresoberfläche (19) der nullten Näherung eine Stromlinie wird, überlagern wir noch die an 
der r-Achse gespiegelte Quelle. Als erste Näherung y®'(r, z) für die Stromfunktion setzen wir an: 

Dee 3er 2 (20), 


Seite 339. Bei unserer Wahl von y(r, z) ist jedoch nicht mehr zu erwarten, daß z() = 0 auch 
noch eine freie Oberfläche ist, daß also auf ihr de Bernoullische Gl. (8) gilt. Um (8) zu 
kontrollieren, berechnen wir für die erste Näherung (20) die Geschwindigkeitskomponenten, die 
wir mit u®(r, 2) und w®V(r, z) bezeichnen. Wegen (17b) findet man 


2 
nl) — = 0) 


um (r,d)= 2 v®=0(e2). 


2) Das Konvergenzgebiet ist am Einfachsten aus 4, zu ersehen. 
3) Unter dem Symbol O(zm) sollen Größen verstanden werden, die mit dem Faktor « oder einer höheren 


Potenz von x behaftet sind. 
23 
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Mit den Ausdrücken für «) und w“ gehen wir in die linke Seite der Bernoullischen Gl. (8) 


ein: 
Man find [u®(r, zu)E + [ur u) = — 2924. 
lan findet 


3 lr2)\ 2 
(a) 40.@u)S3r2gar. en ee 


Wir sehen, daß die nullte Näherung für die Meeresoberfläche z = 2) — 0 dieser Gleichung ge- 


nügt, wenn man von Größen in e® und höheren Potenzen von e& anne 
Ersetzen wir 27, in (21) durch 
2 Bi ‚(r? en 
U) — — ES DR NEN DO ERS RU ER 
M g r 02), 
so wird diese Gleichung genauer erfüllt, als es bei der nullten Näherung der Fall war. Der Aus- 
druck (22) genügt (21), wenn man Terme in e* und höhere Potenzen von e vernachlässigt. 
Die Gl. (22) definiert die erste Näherung für die Meeresoberfläche. Wegen Hilfssatz 1 ist 
die Forderung IV exakt erfüllt. Man rechnet leicht nach, daß man (22) auch schreiben kann: 


Re a). 


Vermerkt werde noch die aus (20), (22) und 
| (17b) fließende Beziehung: 
yd [% a) —- (0) (e?) wie Say (23). 6 
Wir fassen die Eigenschaften der 


ersten Näherung zusammen: 

Die durch (20) und (22) definierte 
erste Näherung für die Stromfunktion 
y®(r,z) und die Meeresoberfläche z4) er- 
füllt die Forderungen Ibis IV exakt. 

vV®: Die Gleichung 

vn) = 0 
ist erfüllt, wenn man alle Glieder vernach- 
lässigt, die mit Potenzen e°, e°, ... behaftet 
sind. (Dies folgt sofort aus (23)). 
VI®: Berechnet man die Geschwin- 


digkeiten 
1 
=, 
1 d 
vw (r,z2)=—ypV(r, 2), | 
NE, * 
so wird die Bernoullische Gleichung & 


[ud (7, zd)]? + [u (r, zw) = — 292) 


erfüllt, wenn man alle Glieder in e#, &, ... 

vernachlässigt. 
In den Koordinaten r und z lauten die 

Gleichungen für die erste Näherung 


y(r, 2) 
| 


ei ee ee ee 


ri 1% 1—? | 
er FR Ir + (1+ 2)? [r2+ (1 — 2)? 12] 
| | (20b), 

Bild 3. Die Funktionen g,(r? dg ( . = 
Nähe Kane ar n edle 207 nt ns = k (22b). 

a | 


Die Funktion g 4,(r?) = — zeigt Bild 3. 


2:72 
(1 +.72)8 
c) Die zweite Näherung: 
Die erste Näherung (20) und (22) unseres Problems enthält das erste Glied der Entwicklung Rh 
der Stromfunktion und der Gleichung der Meeresoberfläche nach Potenzen der Ergiebigkeit ee 
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e=4ne. Bei der zweiten Näherung kommen die zweiten Summanden dieser Entwicklung dazu. 
Die Stromfunktion % (r, 2, e) unseres Problems wird aus Symmetriegründen nur ungerade Po- 
tenzen von e enthalten, die Gleichung der Meeresoberfläche nur gerade Potenzen von e. Wir 
setzen daher als zweite Näherung an: 


yalr,z) = y (r,2) + L;(r, 2) e? A ee EEE DA), 
Be) a ES N a ne (25); 


Die zunächst noch freien Funktionen L, (r, z) und g, (r?) werden durch Forderungen an die 
Stromfunktion y® (r, 2) und die Gleichung der Meeresoberfläche festgelegt. Die Stromfunktion 
y'® (r, z) wollen wir so bestimmen, daß die Forderungen I, II, III von Seite 339 exakt befriedigt 
werden. Die Gl. (25) der Meeresoberfläche soll streng der Forderung IV genügen. Die Be- 
dingung V und VIsollen dagegen nur näherungsweise erfüllt werden, aber genauer, als es bei der 
ersten Näherung der Fall war. Anschließend an die Eigenschaften V® und VI® der ersten 
Näherung verlangen wir: 

V@®: Die Gleichung 


y2 (r, 22) = 0 
ist richtig, wenn man alle Glieder vernachlässigt die mit Potenzen e°, &’,... behaftet sind. 
VI@: Berechnet man die Geschwindigkeiten 
D) 1 p 1 D) 
ua (r,2)= Br vr, 2), wR(r, 2) = r vr, 2), 


so wird die Bernoullische Gleichung 
[u'% (r, za? + + [uw (r, 29)]? = — 29 2® 


erfüllt, wenn man alle Glieder in e®, &®, - - - wegläßt. 

Wir wenden uns zunächst,der Konstruktion von L; (r, 2) zu. Dazu gehen wir von der For- 
derung V® aus. Wegen (24), (25), (20), (22) und (17b) erhält man für y9% (r, 20)) = 0 die Ent- 
wicklung 

ei AT ne. DEE +0 (&) — 
9 


Führt man die in V'® angegebenen Vernachlässigungen aus, so ergibt sich als Bestimmungs- 

gleichung für Z, (r, 0): 

© (r?) 
y2 


4 
x Li(r,0)= — = (26). 


or?) 


Man kann nun leicht nachrechnen, daßsich —, in der folgenden Weise als Linearform der 
©; (r?), ©3 (r2), @, (r?) schreiben läßt: $ 


au 1° nn [4 [OP (r?) + 4; (r?) 04 (r?)] $) Er LER (27). 
Damit hat man 
L, (r, 0) = Toast ERIAEN N Ca am nl A aa) Pe (28). 


Die durch (24) definierte nn EAN (r, z) soll die Forderungen I, II, III exakt erfüllen. 
Die Funktion Z, (r, z) muß dann eine Lösung der Differentialgleichung (2 )sein, die im Strömungs- 
gebiet (z <0) keine Singularitäten aufweist. Weiter muß sie längs der r-Achse der Randbe- 
dingung (28) genügen. Dies leistet offenbar die Funktion: 


Le, 2)= [40% (r, 2) — 4.29 (r, 2) + 2%(r, 2)] 


3 
Die Summanden der rechten Seite sind nach früherem Lösungen von (2). Sie besitzen eine Sin- 
gularität im Punkte Q,, der außerhalb des Meeres liegt. Für 2 = 0 geht wegen (14) die Gleichung 
(28) hervor. Nach einem bekannten Satze aus der Theorie der elliptischen Differentialgleichungen 
ist L, (r, z) eindeutig bestimmt. 

4) Wegen der Indizes bei den Funktionen L(r,z) und g(r?) s. (58) und (55). 

5) Siehe Hilfssatz 4, S. 350. 
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Damit lautet die zweite Näherung für die Stromfunktion: \ 
de 3 


y% (7,2) = & [2 2) — 2%, 2)] + 429,2) — 4029 (7,2) +R® Kr, 2)]. . (29). 


1059 
Sie erfüllt die Forderungen Ibis Illexakt. 
Wir behandeln nunmehr die Forderung VI®, Die Berechnung der Geschwindigkeiten 


1 1 
(= — ya, Wen) = u) 


ergibt 
ra * = (7) 
re 


w&r, 2) =0O ne 


Setzen wir für zden Wert (25) ein, so erhält manfürdieBernoullische IHRE 


[u (r, zn]? 7 KRIEG = — 292 


uw (r,2)=2e 


= Ste ni +0 (622) +0 (e32) 


die Entwicklung 
16 &* [4w, (r?) ©; got (r?) @, (r?) Be (r?) ©; (r?) VEN); 2) 24 
1059 r? y2 pR |+0@) - SI 


Auf beiden Seiten haben sich die mit e? multiplizierten Glieder aufgehoben. Vernachlässigen 
wir die mit &®, &, - + behafteten Terme, wie es VI®) verlangt, so bestimmt sich die Zusatzfunk- 


tion g,(r?) zu 
8w,(r?) [Aw;(r?) Aw,(r) © u 
FEN SE 1 3 4 EN ER RN at : 

Is (r ) 105 g? y2 Is y2 y2 (30) 


Die eckige Klammer läßt sich auch als Linearform der &, bis &, umrechnen ®): 
40; (r?). Aw, (r?)  @(r? 

IN +5) + 160) +00) + Ro) + Raul) 
Damit lautet die zweite Näherung für die freie Meeresoberfläche: 


= gl) + Klr)et = 2m, (r?) {o,(r?) +@o(r?)} : | 


31). 


AT o; (r?) [o; (r?) +50, (r?) + 16 @; (r?) +40 0; (r?) + 72 0, (r?) + 72, Ge |. 


Sie erfüllt die Forderungen IV, V®, VI®. In den Koordinaten r und zheißt die zweite Näherung: 


nd)=e ne la 
FR U I 
121— 487? (1— 2)? ee SCH | 
elle +1 — 27” 
PA ar N 2 Ar Te 
ee ewo, =) ee) 


Die Funktion g? q,(r?) ist in Bild 3 dargestellt. 
d) Die dritte Näherung: 
Obwohl klar sein dürfte, wie das Näherungäverfahren weiter läuft, wollen wir noch auf die 


dritte Näherung eingehen. Ausgehend von der zweiten Näherung (24) und (25) machen wir für 
die Stromfunktion und die Gleichung der Meeresoberfläche folgenden Ansatz: 


y%r, 2) = y&lr, 2) + Lg(r, 2) & A). De Re > (32), 

Dee leere (33). 

Wie bei der zweiten Näherung sollen die noch freien Funktionen Z,(r,z2) und g,(r?) so 
bestimmt werden, daß die Forderungen I, II, III, IV exakt erfüllt werden. Die Forderungen V 


und VI sollen besser befriedigt werden, als es bei der zweiten Näherung war. An die Stelle der 
Bedingungen V@ und VI® treten: | 


V@®: Die Gleichung 
vl) 
6) Vgl. (48), 8.349. 


€ 
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ist erfüllt, wenn man alle’Glieder vernachlässigt, die mit Potenzen &”, &°, - + - behaftet sind. 
VI®: Berechnet man die Geschwindigkeiten 


1 
nd, = 2). 
so wird dieBernoullische Gleichung 
[ui9e JE + Tun = — 2a 


erfüllt, wenn man alle Glieder in &®, e!P, - - - vernachlässigt. 

Die Konstruktion der Funktionen L,;(r,z) und g,(r?) erfolgt in analoger Weise wie 
bei der zweiten Näherung. Zur Festlegung von Z/, (r, z) gehen wir von V'% aus. Eine einfache 
Rechnung liefert für pr, 29) = 0: 


Mod), LIE) IE, ren 
re a HL) +0) =0. 


Die Glieder mit e? haben sich wegen (27) aufgehoben. Die in V® genannten Vernachlässigungen 
bestimmen 2; (r, 0) zu: 


n.8 | 49(N)w;(r) 4wilr?)w;(r?) |, @i (r?) @; (r?)\ 
bs (r, 0) = 35 g? | 72 — ars - H 72 | le Mare (34). 
Wie bei der zweiten Näherung, so kann man auch jetzt die Zusatzfunktion Z, (r, 0) als Linear- 
form der ®, (r?) schreiben. Es ve nämlich die Rn Identitäten 5): 


ı(r? .(r? 5 1 
o(r n (7?) = 502 u 0) +50 De a 2.) — 50 () 


© (r?) w, (r?) 200 "2 112 136 25 
ri = 3003 @=(r ) + 30050 a) tz) | 


y2 
! 
7 (r? DD); 
— 1570) (35) 
o (r?) ©; (r?) 64 „2 64 2 8 2 8 y2 12 2 
y2 ıu= 1001 sg (r ) ir 1001 @3 (r )+ 143 © (r + 231 0; (r ) + 1001 &g (r ) “+ 


@g (r?) 


fi 
2 tar 
oe) varagas 


Damit wird aus (34): 
8 ee 1984 1984 2648 664 


Par re ee re, 2 2 
ur 0) = I Il rer | 
20 39 1 EoP): 
+ Nr ee ) 1 
Genau wie bei der zweiten Näherung kann man auf die genaue Gestalt von 2; (r, z) schließen: 
8 SE 1984 2648 664 
(2) (3); a (4) .ı u: (5) (y 
Bel) za es nr ne ers 69 
. (37). 
(6) 
r n ee 5 Be ) ) 


Wir erkennen, daß L,(r, z) eine Lösung von (2) ist, die keine Singularitäten im Strömungsgebiet 
besitzt. Die Stromfunktion y(r, z) ist damit bestimmt. Sie erfüllt die Forderungen I, II, III 
unseres Problems exakt. 

Zur Ermittlung des Zusatzgliedes 910(r?) in (33) stellen wir die Bedingung VI®) auf. Man 
findet für die Geschwindigkeiten: 


un, 2) = ydlr, 2) = Ka z ande) +4) —osl®)] 
ale) Ha Ba) + to) re] 
8 © [1984 2648 664 
tr AN let rl) 


a) + Se) +0 em) +0) +0Le 
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1 2 € do, (r?) 4e° 4 do,;(r?) 4 do;(r?) 1 do,(r?) 
ae BT Det 2 m 
u) tr (r 2) r dr ru 1059 FE UR DEAL r dr 


+0(e) +H+HO(E)+O(E). 


Wir setzen nun 


= 2: =) + gelr?) ei guo(r?) € 


ein. Für 9, (r?) wird der Faktor von e? auf der rechten Seite von (22) benutzt. Für 9, (r?) ver- 
wenden wir den Ausdruck = 


u (r, 20) = 28 T 05, Hol) +4) a4 
Eee +) 
+ En 0) — “u 0,(1) a 0) + | +, 
ee ee, 
DieBernoullische Gleichung der Bedingung VI®) bestimmt dann g,0(r?) zu: ‘ 


ENT ae | ee (r?) @; (r?) 2.:08(r2) 
Pie , ; Bee 12 N 3 Hi 4 >) 
TYıo (7 )= Akt: 5er y2 { 40; 72) + 1 @,(r? =) 0; ‚(r? )} 5 105 ri { 0) (r ) 
Br B 2 @,(r?) [1984 8 2.648 
+4 0; ur) Ne @.(r?)} 35 z E00 (@3 (#2) u © (r?)) en 15015 ©; kn) 


664 20 32 1 & 
Fr] 
[o;(r?) do,(r?) 1 E 4 dos(r?) 4 doz(r?) , 1 N 


ur dr | 105 Y dr Y dr or dr )1° 


Damit ist die dritte Näherung (33) für die Gleichung der freien Meeresoberfläche gefunden. In 
den Koordinaten r, z lautet sie: 


( (1?) . 
le) = 7 ar 


11904 r!6 + 120640 r!* + 304 256 r!? + 43855 rl" 
— 49443515 r?® — 73 667 636 »° + 13 163060 r? — 210624 r?]. 


Bild 4 zeigt den Verlauf von 9° g,0(7?). 

Vergleichen wir die Darstellungen 3 und 4 für g9;(r?), 92 g6(r?) und g?gjo(r?) so er- € 
kennen wir, daß wegen g—= 981 cm sec”? die Funktionen g4n + 2(7?) (n = 0, 1, 2) mit wachsendem 
nimmer rascher mit der r-Achse zusammenfallen. Für kleine r-Werte zeigen gg, 9g und gj0 ZU- 
nächst Schwankungen um die r-Achse mit kleiner werdenden Amplituden. Die hohen r-Po- 
tenzen im Nenner von 93, 44 910, lassen diese Funktionen für größere r-Werte asymptotisch Null 
werden. Nach Bild 3 ist die r-Achse für q, (r?) bzw. gg (r?) etwa bei r = 6 bzw„r = 3,25 erreicht. 
910 (7?) fällt mit der r-Achse etwa von r = 2,5 ab zusammen. Setzen wir die Konvergenz des 
Verfahrens voraus, so werden wir vermuten können, daß die Strömung für nicht zu große Quell- 
ergiebigkeiten e für etwa r = 3,25 durch die zweite Näherung gegeben ist. Die dritte Ns 
dürfte die Strömung für r > 2,5 bestimmen. 

In den Bildern 5 bis 7 ist die erste bis dritte Näherung der freien Meeresoberfläche für die 
Quellergiebigkeiten e — 85, 120 und 170 g sec”! gezeichnet. Im Falle e — 85 dürfte die gesuchte 
Meeresoberfläche durch die dritte Näherung bestimmt sein. Für die beiden letzten Fälle sind 
höhere Näherungen erforderlich, um die bei r — 0,75 liegende Sogspitze zu erfassen. 


4. Gestalt der @, (r?) 


Es erhebt sich die Frage, ob das bisher geschilderte Verfahren zur Konstruktion der ersten 23 
bis dritten Näherung beliebig fortgesetzt werden kann. Wir werden den Nachweis in 7. durch 
Induktion erbringen. Zuvor ist es in 4. und 5. nötig, genauer auf die Funktionen @,(r?) 
einzugehen, 
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99 (r?) ee 


Bananen 1.Näherung 
nn he ER : 
3 D 


Bild 6. Erste bis dritte Näherung der Meeresoberfläche im Falle 
2 
r =0,1. Quellergiebigkeit e = 120 y see! 


Bild 4. Die Funktion q,,(r?) derdritten Näherung der 
Meeresoberfläche (g = 981 em/sec”?) 


Bild 5.*Erste bis dritte Näherung der Meeresoberfläche im Falle Bild 7. Erste,bis dritte Näherung der Meeresoberfläche im Falle 
2 


2 . 
= = 0,05. Quellergiebigkeit e = 85 g sec” = = 0,2. Quellergiebigkeit e - 170 gsec* 
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Die Funktionen o,(r?) lassen sich leicht explizit mit Legendreschen Polynomen P. (0) ir 


\ angeben. Wir entwickeln dazu : | aa 
RE 1—2z 1— 3 1 - 

Er (0) 2 ; 

e“ 29 (r, z) + (1 — 21127 mr -L ge f RR, 1 ee IE ( DIR ni 

& | N1V+rWN+r VERF 

I: | nach steigenden Potenzen von vı ni unter Beachtung der Beziehung j Ser 

BR Nyt$ 

34 S 3 % DEE i 

ar PsX2)h® hl<Miniz +/2—1l- 
: ESEEN, =>, ie) a] < Min |o Ya —il 


Setzen wir P_,(x)=0, so findet man: 


Re na 
er 
2% 


BE Su 7 1 | Er 
$ = ——J—-t1-+r Pal | $ 
er Fade 5 rs | a: 

® N | 
= Koeffizientenvergleich mit der Reihenentwicklung für 2% (7, z) in (17a) re 
Bi: . 1 1 Bi & 
= Be ze ler. Be ) DEE 
% 5 J 2 
Ri); : o,(r?) —n! y1 Au ne nie 1 BR P E Aare Fl (39). j h 2 j 
(1 +r)2 


Der Definition der &,(r?) entnimmt man leicht, daß man sie auch in der Form 


= 
RE 


Z,0r2 Gt x 
an a N al 
1i+r) 2 ER 
schreiben kann, wo Z,(r?) ein Polynom in r? ist. Nach (39) ist RR Zaherpolynom bestimmt BR 


durch 


EEE RE MEN HUN TERR 


1 1 20 Dr 
Zr) =n!i+r2 ®[p Pu oe Bl 2 Be 
el) RR an) 
Führt man in die rechte Seite die bekannten Entwicklungen für die Legendreschen Poly- > 
nome ein ?), so findet man nach etwas mühsamer, völlig elementarer Rechnung: { R 
Zan(r?) N Iop = Hi m=6, Do, | 2 22 
Mr | re: Ä . (42) 
Zum-ı() = PR Z- e n ur ap $ m—l, 2, De BR y 
RS ER Karl ee. 

ZD(— 1) = (— 1) mr ar (1 ) EEE RAR ey. 
| a m 

- Dabei bedeutet hier und später | i u As ia BR RL N Ye 


(= 1, (= +1) «+ 2-. at 
& beliebig, k positive ganze Zahl>1. / u 
Aus den Formeln (42) und (43) folgt: Die: Polynome z, en). n=0, ie 


£ 1 
vom Grade er = ] . Mit + e ist die größte Rn Zahl bezeichnet, 2 


 Zu(2) = ZW (x) an der Stelle «= — 1. Für die höheren Ableitungen 
zo (— 1 RER 


?) Siehe z. B. Magnus- Oberhettinger, Formel 
matischen Physik, 2. Aufl; (1948), 8.69, 4. Zu d. j Forsäe 


P) 
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5. Spezielle Entwicklungen | 
Wir benötigen im folgenden die Entwicklung gewisser Klassen von Funktionen nach den 
@,. Zunächst beweisen wir: 


= do, (2) ST @;(2) a 
ENT de en n=0,1,2... (44). 


Die Berechnung von a aus (39) ergibt unter Benutzung der Rekursionsformel der Legen- 


dreschen Polynome: 


1 
er 
2 do,(t) _ SEAN + x 
(r+1)! da ar 
+9)? 
Bildet man die rechte Seite von (44) mittels (39), so erkennt man die Richtigkeit der Behauptung. 
Die Gl. (44) findet in der folgenden Form Verwendung: 


N 


1 do,(r?) So; () ) Ir 
n = ! j a 45). 
s iR (n +1) er (45) 
Weiter gilt 
il N 
re meND ro, (e n=0,1,2,... (46). 
(1-+x2) 2 ind 
Die Vorzahlen Br auf der rechten Seite bestimmen sich durch: 
E i j 1 2n—2 £ 
Be GH): i=0,12%,...n (47). 
2nm(- ) on 


1 
1.3.5»... (2n—1) 
Induktion. Man überzeugt sich sofort, daß (46) für n = 0 gilt. Gilt die Formel für die Zahl n, 
so differenziere man (46) nach & und setze auf der rechten Seite die Ausdrücke (44) ein. 

Wir wenden uns nun der folgenden Aussage zu: 


= Yejoy(a) n=1,2,3,... (48) 


(n>0). Der Beweis erfolgt durch 


Insbesondere it =} = 1, Bf = 


c 
=( 


mit 


| i REN, 
N Me J 5 a 1=0,1,2,..,n—1 (49). 
Man rechnet sofort nach, daß (48) für n = 1 gilt. Nimmt man die Formel bis zum Index n als 
richtig an, so gelangt man zu der für die Zahl n + 1 unter Verwendung der aus (16) fließenden 
Relation: 


ER = @n+1(%) — (2n +1) —— a +Rr—1)r +1) m b 


en 


Die Gleichungen (45), (46), (48) zeigen uns, daß man die Ausdrücke 


1 do, (r?) 1 ©,(r?) 
Bee din, © al r? 
(+) > 
als Linearformen der &,„(r?) mit konstanten Koeffizienten darstellen kann. 
Die bisherigen Ergebnisse setzen uns in den Stand, eine Reihe von Hilfssätzen abzuleiten. 
Hilfssatz 3: Für a, ß, y>0 gilt: 
“+ß+y+1 
0 (2) ©5(2) ©,(2) — 2 a,0;(8). 
er2, ; 
Die0,1] = 0,15: .,0 Y #-+y-+-1) bedeuten Konstante. 
Für «+#-+y21 verschwindet a,. 
Für +28-+9»22 verschwinden die Koeffizienten a, und a.. 


’ 


Be a EN 
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Zum Beweis beschreiten wir nicht den einfachsten Weg. Wir wollen — über den Hilfs- 
satz hinausgehend — die explizite Darstellung von ®, w, w, als Linearform der &; angeben. 


Wegen (40) und (42) können wir für &,(z) schreiben: 


eazocı) 
0,(8) = | 2 +1 . 
a or 
Mit E 2 Se a(a=0,1,2,...)und (43a) läßt sich diese Gleichung umformen in: 
| 155° 

ar Ze 1) 

9, (R) = D FE 26 +1 
dt) ® 


Für das Produkt ®, ©; w, findet man dann: 


4 h : a+ß+Yy vjk Be ee s) = 
9,(2) @3 (X) ©,(z) = 3 ® By a]; (1-+a) \ 2 >00): 
so 
Dabei steht K j 3 für die Summe: 
ö aßyls h U) 
(si) _ ZEN ZPENZW1) 
ee ! k! 
Nunmehr setzen wir (46) in die rechte Seite von (50) ein. Man erhält . 
Dee tel. 
9,(%) (X) w,(x) = > 0 BD ee a ER 
t=0 
mit 
a+ßtyti/jgk x 
dy = >, [ B ‚). RR» . . . . . . . . . . . . .(52) 


‚=t 


Damit ist der erste Teil des Hilfssatzes 3 bewiesen. Ist « +ß+y2z]1, so ist mindestens einer 
der Indizes > 1. Bilden wir (51) für = 0, so bleibt wegen Hilfssatz 2 übrig a, — 0. — Um zu er- 
kennen, daß für « +ß-+49y>2auch a, = ist, bringe man (51) auf den Hauptnenner und ver- 
gleiche die Koeffizienten der höchsten Potenzen von (1 + x) auf beiden Seiten. R 
Hilfssatz 4: Für a+$22, y2>1 gilt: 
&HB+y +1 


©, (2) ©E(2)@,(®) _ > b;ow;(x) 6b; Konstante mit = b, & 


T 


j=2 
Wir dividieren (51) durch x und setzen auf der rechten Seite (48) ein. Dann bekommt man: 


a+ß+y+1 


ae Dil be well ee a 


o=0 


+ > +1 
T=06 
Man sieht sofort, daß 5, verschwindet. Setzen wir nämlich in (53) für x den Wert Null ein, so 


bleibt auf der rechten Seite db, @, (0) = b, übrig. Die linke Seite verschwindet, weil auf Grund 3 
unserer Voraussetzungen «+ ß>2,y>1 mindestens zwei der Indizes a, ß, y größer als Null 


mit 


sind. Der Ausdruck Zu geht wegen Hilfssatz 2für —Onach Null. 


Um die Behauptung b, — b, zu bestätigen, bringe man (53) auf den Hauptnenner und | 
vergleiche die Vorzahlen der Potenzen (1 + z)*+?+r+2 und (1 + »)**P+r+1, Der Hilfssatzz 4 
kann zur Aufstellung der Formeln (27) u. (35) verwandt werden. Ausgehend von Hilfssatz 3 Be 
beweist man leicht durch Induktion: eh 


ÜBESFER: er 
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Hilfssatz 5°): Das Produkt &,(®) ©; (2) ©5,(@) * - - .  @p,,(#) aus einer ungeraden Anzahl 
2k-+1 von Funktionen ©,(2), ©,(2) @=1,2,..., 2%) läßt sich als Linearform der Funk- 
tionen © (2), ©, (8), . . - ; Oa+p + +B,, rk) Mit konstanten Koeffizienten e; schreiben: 

at&ß;tk 
(2) ©; (2) @8,(%)* .. .* @p,,(2) = 22 6;0; (8). 


6. Abkürzungen 


Wir haben es in Zukunft häufig mit Linearformen der &,(x) zu Lun, wobei der Wert der 
konstanten Koeffizienten der ®,(x) keine Rolle spielt. Um diesen Tatbestand auszudrücken, 
wollen wir ein besonderes Symbol einführen. 

Wir verstehen unter /,(x) (* = 0, 1, 2; u > i) eine Linearform der ©;(&), ©i+1(%), ®i+2(2), 

.,.®,(&) mit konstanten, beliebigen Koeffizienten «a,, die zum Teil oder auch sämtlich den 
Wert Null annehmen dürfen: 


RER [% 
1,(%) er 2, q; 9;(8). 
jei 


Mit dieser Definition lassen sich die Gleichungen (45), (46), (48) folgendermaßen schreiben: 


I 1 do, (r?) | n ‘ 
var NL [9 y2 Hz AR ER 4! a), 
r Ir nr) en = OT, (45a) 
1 0 9 N 9 46: 
ar Kae re 1% (r?) WM (, iR PAR EIN ( 5a) 
(+42): 
or: r ß 
n en IE In (r?) U 13 2, 3, N en Be te de (48a). 


Nach unserer Verabredung können wir die rechten Seiten von (46a) und (48a) durch das gleiche 
Symbol bezeichnen. Durch /, wird nur angegeben, daß eine Linearform der @y, . . . , @, vorliegt. 
Über die Koeffizienten wird nichts ausgesagt. 


Eine Folge von (48a) ist offenbar 


EEE RETE DE 
u) a ae Aueläb): 


Die Hilfssätze 3 und 5 geben zu folgenden Formeln Anlaß: 
0,2) @5() @,(8) = 1ars+y+1(%) ,ß,yZ0 
> ©,(2) ©5,(2) @p,(2) ° .. . " @p2 (0) = lara +... 48,44% (®) 20. ... (94). 


Außer den Linearformen ], verwenden wir später noch quadratische Formen q,(x) der o;(2). 
Diese seien durch 

uM)= 3, ui; @i(@) o;(®) MHZ2L 2.2... (8) 

1si+rjsu 
definiert. Wie bei dem I so soll es auch hier auf die reellen Konstanten «,,;,; nicht ankommen. 
Sie dürfen zum Teil oder sämtlich den Wert Null haben. In jedem Summanden der rechten 
Seite tritt ein Produkt aus genau zwei Faktoren ©; (x) @;(x) auf. Für die Indizes :, jgilt:« +71, 
‘i+7<u. Die Funktion @,(x) tritt in einem Symbol q,(x) also höchstens linear auf. Wir 
können auch sagen: Die Funktionen g, (x) gehen mit 2 — 0 mindestens wie znach Null. 
Wir vermerken noch die aus Hilfssatz 3 zu gewinnende Rechenregel: 


94 (%) 9,(2)= luro+1(%) EA a Le WERNER (56). 
Mit (48b) gilt weiter: 
ent Fo dan 4 e(*) RE I Et I EEE ne (48e). 


®) Die Formeln (44), (46), (48) und die Hilfssätze 3, 4, 5 sind Spezialfälle eines allgemeineren Entwick- 
lungssatzes, auf den wir hier nicht eingehen. 


.. V1® und VII® dadurch entstehen, daß man n durch n + 1eersetzt. 
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Quadratische Formen q,(r?) sind uns bereits in den Gleichungen (22a) und (30) für die erste und 
zweite Näherung der Gleichung der Meeresoberfläche begegnet. Sie wurden dort der Definition 
(55) entsprechend bezeichnet. Ebenso läßt sich der Ausdruck g,, (r?) von S. 346 als quadratische 
Form (55) mit « = 10 umrechnen. Man muß dazu die in 5. entwickelten Formeln benutzen. 


7. Die n-te Näherung 


Wir wenden uns nun wieder unserem Problem der Jungfernquelle zu und definieren, was 
wir unter der n-ten Näherung verstehen wollen. Betrachtet man die in 3. gewonnene erste, 
zweite und dritte Näherung, so liegt es nahe, als n-te Näherung Ausdrücke der folgenden Form 
anzusetzen: 

a) für die n-te Näherung der Stromfunktion: 


yR (r, 2) = E[lAo(r, 2) — 29% (r,2)] + e L(r,) +» N (57) 
+1 sl )H Her Innen 2) J : 


Die L,„_4(, 2) (u=2,3,...,n)sind dabei Linearformen der 2% (r,2), 2% (r, 2),..., Q4rZN (7, 2) 
mit konstantem Koeffizienten &4.—4,,», die auch den Wert Null annehmen dürfen: 
a 


Le Dr Fe) u 29. RE NOIR 


v—_= 23 
b) für die n-te Näherung der Gleichung der freien Meeresoberfläche: 


ze gl?) Ele) + ganze (fr?) ER F gms lr)er: 22 (59), 


Die 94. —2 (r?) bedeuten quadratische Formen in den @; (r?), wie wir sie in (55) definiert haben. 
Die Ansätze (57) und (59) erfüllen die Forderungen I, II, III, IV unseres Problems exakt. 
Wir wollen die so definierten Funktionen y®) (r, 2) und 2) eine n-te Näherung nennen, wenn 


folgenden Bedingungen genügt wird: 
V®): Die Gleichung vr (r, 2) = 0 
ist erfüllt, wenn man alle Glieder vernachlässigt, die mit Potenzen e?” +1, e®”+3,,,., behaftet sind. 


VI®): Berechnet man die Geschwindigkeiten 


1 1 
dl, u), 
so wird die Bernoullische Gleichung 
r) 2 Pr EL n 
[ur Cr, 20937? + ot, =]? = — 2a.2h 
erfüllt, wenn man die Glieder in &?*+?, g’%+#, ... vernachlässigt. 


VI1I®: Die Funktionen Z4„—4(r, 0) (u = 2,3, ..., n) gehen mit r— 0 mindestens wie r? 
nach Null. 

Man prüft leicht nach, daß die bisher aufgestellten Näherungen die erste, zweite und dritte 
Näherung im eben definierten Sinn darstellen. Daß für diese auch die bisher nicht gebrauchte 
Forderung VII gilt, erkennt man mit Hilfssatz 2 aus den Gleichungen (26) und ee: 


8. Der Induktionsbeweis 


Zum Schluß der Arbeit zeigen wir, daß das in 3. entwickelte Näherungsverfahren unbe- 
beschränkt fortsetzbar ist. Der Beweis wird durch Induktion erbracht. Für ce (n + 1)- ste 
Näherung der Stromfunktion und der Meeresoberfläche setzen wir an: 


yatD(r, 2)= y” (r, 2)+ L(r, 2) gen+l N PEN RE ea (60), 
Zur = ger) ent? E & a EW LS DE TIER (61). 


Die Funktionen Z (r, z) und q (r?) sind so zu bestimmen, daß (60) und (61) die Forderungen I], II, 
III, IV exakt befriedigen. Weiter müssen Vr+D, VI®+D und VII®+D erfüllt sein, die aus V®, 


Zur Festlegung von L(r, z) gehen wir von V®+2 aus. Danach muß für die (n + 1)-te 
Näherung die Gleichung 
yr+D (r, zu +D) — 0) 


a REN WR EI RE WR MN ©... 


u en ee 


N 
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erfüllt sein, wenn man alle mit e2”+3, e2”+5,,.. behafteten Glieder fortläßt. Mit (60) und (61) 
lautet (62): 
Yin Ir, m) + g(r 2).g 2n +2) + L(r, zin NE? RL N), Baer A (63). 


Wir entwickeln den ersten Summanden nach q(r?) e?"*?, den zweiten nach z(% +. Wegen 
a 2 ". 
+ = O(e2) folgt: 


En B 
ya (7, 200) + L(r,0) einst I gen? LO) =0 2... (69) 


pr (r, 2 
r? 


Die Entwicklung von y (r, 2) und - ee) 


und weiterer später benötigter Größen nach 
Potenzen von e ist zwar elementar, aber mühsam auszuführen. Der Rechnungsgang wird im 
Anhang geschildert. Man braucht dazu wesentlich die Hilfssätze und die Gl. (45a) und (48a). 
Man entnimmt den Formeln (A 12) und (A 17) des Anhangs: 


ap") (r; 2) 
BOTEN) = ln (ee Hl O0 (ein #P), So air le) N, (DAL 
Die Entwicklung von 9" (r, z0)) beginnt wegen V® erst mit einer Potenz e?” +1, Setzen wir (65) in 
(64) ein und vernachlässigen alle Terme in e?* +3, e?” +5, ,.., so ergibt sich als Bestimmungs- 
gleichung für Z(r, 0): 


L(r,0)= — in (r 2) Sa, ;@;(r?) &4n,; Konstante (66). 


j=2 
Diese Gleichung entspricht den Randbedingungen (28) bzw. (36) für die Funktionen Z, (r, 0) bzw. 
L;(r, 0) der zweiten bzw. dritten Näherung. Wie wir dort ausführlich dargelegt haben, so kann 
man auch von (66) aus auf die genaue Gestalt von Z (r, z) schließen. 


4n 


L(r, 2) = Z,5un. 20 ENTE Fe a RER NENRA TE 


Wegen (58) kann L(r, z) als eine Funktion L(r, 2) = Ly. (r, 2) geschrieben werden. Wie im An- 
hang im Anschluß an (A 12) gezeigt ist, geht die Funktion li,(r?) in (66) mit r—0 wie r!nach Null. 
Die Stromfunktion y* + (r, 2) der (n + 1)-ten Näherung erfüllt damit die Forderungen I, 
II, III, VII®+D exakt. Sie wird durch die Stromfunktion und die Gleichung der Meeresoberfläche 
der n-ten Näherung bestimmt. 
Für die Berechnung des Gliedes q(r?) e?” *?in (61) ist die Entwicklung der Geschwindigkeils- 
komponenten u” +2 (r, 2, FD) und w® + (7, 22 +D) erforderlich. Es gilt: 


1 [Oym (re) , Alın (rn 2) 


nr +1), nr +1)\ — ern tl 

U (13 Zu ) y | 02 02 a a") +alr 2) 2 +2 

109y9 (r,2)) 1 pm (r, 2) EUR 1 nn (r, 0) Rah] ng 
Ar ya an Hz An 7 an Hr lg2n+3 4 (00), 

Y 02 jr r 02? re 'r 02 x 7% ) \ 

1 Py (rn, 2, 1 0L,.(r, 0) 
MN) (n) ae re 2\ 720 +2 % gan +l 6) 2n +3) 
um (r,ay)t, a Nee +0(e 


Aus (67) folgt nun wegen (11): 


4n 
as 2) BR . Q6+D un 2). 


un 
Mit (14) hat man dann 


n(r, 0 
oL, oe 0) = Yoamsarıl)= — Kur) ei 33 EEK 


Fey! 
Setzen wir in (68) die Formeln (69) und (A 21) ein, so gewinnen wir: 


UM +N (y, zn +0) = u (r, are, are) nn. 2. (70). 
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Bei Berücksichtigung von (A 18a) findet man für das Quadrat von (70): 


4, (7?) lin+i (r?) 


y2 


een Here Zi 


[ur +D (7, Bar = [u (r, 2n))]2 ie 


Das mittlere Glied der rechten Seite können wir wegen (48c) als quadratische Form g%. +2 (r?) 


umrechnen: 
4, (r?) Dan +1 (r?) 


y? 


= an Kar?) RE u ne 
Dann wird aus (71): 

[ur HP (v2 Pe Tl I in sl Er Oz 
Wir bilden nunmehr w* =D (r, En 7»): 


1 oy) (r, 2) | OLyn (r, 2) & 


wir + (r, Zu z r, 2n+1 
DR S {Fr or or 2= en +q(r?) 


2n +2 
ent 


(n) (» >(n) 
a 9Y (r, ZM ) ex 0(& n + is) — pn) (r, m) - ) (e?" +2) : 
m or 


Für das Quadrat erhält man mit (A 30) 
[win +3, 2 D) JE = [wmn, zm Je + Ole +) 2. 2222. (MM). 


Die Bestimmung des Zusatzgliedes q(r?) "+? in der Gl. (61) für die Meeresoberfläche erfolgt | 
nun mittels der Forderung VI” +1), Danach muß für die (n + 1)-te Näherung die Gleichung | 


[ur FD (rer DYJa = [uote Hl, An HE E75) 3 


erfüllt sein, wenn man alle Terme fortläßt, die mit Potenzen e?"+#, &2"+6, ,,, behaftet sind. | 
Wir gehen mit (61), (73) und (74) in (75): | 


[u Cr, 1° + [u e, 0 )]° + gina?) 0?" +? +0 (er) = — 2g m —2galer) ern+E (76). 


In die linke Seite führen wir den Ausdruck (A 32) ein und berücksichtigen die Tatsache, 
daß die n-te Näherung 2%) der Forderung VI" genügt. Dann heben sich in (76) alle Terme in 


b 
| 
&2, et... , 82" weg.. Es bleibt: | 
ü 
{gin+2 (r?) + gin+a (r?)} er +? 0 (er +) = — 2ggl(r?)e®r+2, 
| 


Vernachlässigen wir jetzt alle Glieder in 2” +#, &®+6,..., soergibt sich als Bestimmungs- 
gleichung für q(r?): 


ET „2 vr Se 2 i / 
29 1(?) = Gin+2(r?) + Gin+2(r?) = Qan+2lr). © 


Damit ist das Zusatzglied g(r?) im Ansatz für die (n + 1)-te Näherung als eine quadratische Form 
J4n +2(r?) erkannt. (61) genügt wegen Hilfssatz 1 exakt der Forderung IV. Damit ist der Induk- 
tionsbeweis erbracht: 

Wir bemerken zum Schluß, daß das Zusatzglied q4n +2 (r?) &®*+? der (n—+1)-ten Näherung 
der Gleichung der Meeresoberfläche nicht allein aus den Größen der n-ten Näherung ableitbar ist. 
Es ist noch zusätzlich die Stromfunktion der (n + 1)-ten Näherung erforderlich. In die quadra- 
tische Form g%n +2 (r?) geht nämlich nach (72) die Linearform 13, +1 (r?) ein, die sich wegen (69) aus 
dem Zusatzglied L;,(r, 2) von y® +2 (r, z) berechnet. Ä 


/ Anhang. 


Al. Die Potenzen von 2{n) 


Die n-te Näherung für die Gleichung der freien Meeresoberfläche ist durch (59) gegeben. 
Obwohl die rechte Seite von (59) aus einer endlichen Anzahl von Gliedern besteht, wollen wir 
dafür formal die unendliche Reihe schreiben : f 


u = ul) + gel?) +... + Que REES - Duo: 2 
j=1. 
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Die quadratischen Formen g4;— (r?) sind also für 7 > n identisch Null zu setzen. Die k-te 
(kZ1) Potenz dieser Potenzreihe lautet dann: 


Bo Se u... (Al) 
o=k 
mit 
(= >, EN a EEE a) OR "do,—2(t?). er IA). 


01, 09,..- 9% > 
%,+0+°*-. E60; 
Wir betrachten einen Summanden 94,—2 (1?) 94,—2 (7?) - * 440,—2 (r?) von (A 2). Setzen 
wir für die quadratischen Formen g4.,—.(r?) ihre Ausdrücke (55) ein, so erhalten wir 


940,—2 (7?) un—2 (7?) 1, )=IB.. .i, (r?) 9; (r?) ©, (r?) @;, ER FE; ©;,(r?) @;,(r?) 
Kae BE oNstantey" a ae, el), 


Y 


Auf der rechten Seite von (A 3) steht eine Summe von Produkten mit genau 2% Faktoren 
©, 8,0, 9 Für die Summe der Indizes gilt wegen der Definition der 94.—s und (A 2): 


Rey netter et — ak... AA). 
Wir betrachten nun einen einzelnen Summanden 
Di, (r?) ©; (r?) @;, (r?) @;, (r?)...0; Bel) o;, (r?) 


der rechten Seite von (A 3). Wegen (A 4) gibt es mindestens einen Index > 1. Es sei i, dieser 
Index. Dann läßt sich auf das aus 2: (k— 1) + 1 Faktoren bestehende Produkt 


o;, (r?) @:,(r?) @;,(r?) . . ..@;, (1?) @;, (r?) 
der Hilfssatz 5 in der Gestalt (54) anwenden. Man findet: ' 
9; (r?) ©, (r?) ©, es ©, (r?) o;,(r) = 4 (r?) De a 
Damit hat man 
©, (r?) 9, meRer 80, (12) ©; (),=.4(r?) u <4Ao —k— % 


Für (A 3) kann man dann schreiben 


110, —2 (r?) 4410,—2 (Mrz Jia, — 2 (r?) = Quo -ı(n) 
und endlich für (A 1): 


Dean) 


A2. Die Entwieklung von x") (r, 2) nach Potenzen von e: 
Wir betrachten zunächst die Funktion 


oo F 2 
Lu( = Do) 209-2 SAN) ur... (a8. 
Enz 1 u en * 
Mit (A 5) erhalten wir dann 


Be me) 


j=1 o=25—1 / 


Umordnen der Doppelreihe nach Potenzen von g ergibt: 


[’ +0 
ST a l® 02. (a7 
22 wagt 
Aus der inneren Summe der rechten Seite greifen wir einen beliebigen Summanden heraus: 
1 ire 
an ra) 140—2;j(1?) l1<sj=s|5 


Zunächst ist klar, daß ein solcher Ausdruck mit r— 0 mindestens wie r* nach Null geht. 
Weiter können wir wegen (56) für ihn schreiben: 
- 1 
TOrF IT @2;—1(0?) 42; (0?) = lie (1?) . 
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ER 94 5 1 
Damit läßt sich auch der Faktor bei e?” in (A 7) als Linearform — > %o (#2),°), der 


@(r?), . . . » @40 (2) darstellen: 
> 


$ 


@s;—_ı(r? N 159% 
2 ze 1-2?) = —— Malr?) 


I: 


(A 7) lautet dann: 
L,(r,2w) = = I all) ER a Ne ner 


Die a;j,(r?) gehen Er r — O0 mindestens wie r er! Null. 
10 


Für die Entwicklung von y”Xr, en ) sind noch weitere Reihen erforderlich. Auf ihre Her- 
leitung soll hier nur kurz eingegangen werden. Man kommt zu den Resultaten, wenn man die 


gleichen Überlegungen anwendet, die zur Aufstellung der Entwicklung von Zo(r, zur) führten. 
So gilt zunächst für v2 2 


2% (r, 2) = Ar, 2 I bio+.(®) E20 — P bio +,(r?) A NEE): 


Die bio +. (? ") gehen für oz 1mitr 2 mindestens wie r* nach Null. Die Berechnung 
von L4u—«(r, ir) (u > 2) knüpft an die Formeln (58) und (A ” an. Setzt man: i 


4u—4 
er Nu—4,v biau102) nr a ) o21 


v=& 


so rechnet man nach: 
oo 
Lu DE) DS Ara)... 3 > (A110) 
o=1 


Führen wir in (A 10) ein: 
4u—4 


Lıu—.(r, 0) = =, Ku,» (r?) = C4u— lt?) 
so gewinnen wir ; 
De - 2 RL Era RT RS 
Als Linearform der bij, +,(r?) gehen die &s+44—4(r?) (co > 1) mit r—0 mindestens wie 


r“nach Null. Nach Induktionsvoraussetzung VII® gilt dies auch für Z,,_—4(r, 0) = 4—ı(r?). 
Damit gehen alle Koeffizienten es +44u—4(r?) (a = 0,1,...) mit r mindestens wie r* nach Null. 


Mit Hilfe von (A 8) und (A 11) sind wir jetzt Faktande, die Entwicklung von 


vr, & Di eLo(r, 2 270) we. (r ir) + +erTin,,Tule ir) + vr 4 en Dinsilh zn) 


nach Potenzen von eanzugeben. Eine leichte Rechnung ergibt: 
DR yon, 2) N Uzlr) eitı— ln ent 10er). 2.2... (AM). 
en j 


Die neuen Koeffizienten l4,(r?) sind als Linearformen der a4,(r?) und c4,(r?) Linearformen der 
Funktionen @, (r?), . . . , @gj(r?). Sie gehen sämtlich mit r — O mindestens wie r* nach Null. Die 
Summe der rechten Seite beginnt erst mit dem Index 7—=n. Nach Induktionsvoraussetzung V®) 


heben sich in y®(r, z5r) alle Glieder auf, die ‚mit 8, &, ..., e2”-1 multipliziert sind. 


A3. Die Entwicklung von [u®/(r, :{7)]® nach Potenzen von s: 
Unter Anwendung desselben Gedankenganges wie in A 2. findet man: 


OL, (r, zur EA TEESE! 
NA =2al)+ Ser)... (AB) 
9 Lofr, 231) 


ee 


Pen 


02° 


®) ) Diei in 6. definierten Linearformen I! sollen im Anhang auch mit anderen Buchstaben a), A v usw. 


bezeichnet werden. Der Buchstabe g ist jedoch ausschließlich den quadratischen Formen vorbehalten. 


( 
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Die Entwicklung von 


Lu?) . 
Ka Ser — 2 u, AP +V(r, 2) u22 
an der Stelle z— 2% "ergibt unter Benutzung von (A 9): 
oL4u.—4(t, cr SUr 
44 u- a RN e2° Der a N ANTD). 
o=0—d 
Eine einfache Überlegung zeigt: 
0° Li dr, 2a N 
N a —0O(1) u EAN NE SR (A 16). 
Für 
: (nyf. (n) I) (n) (n) 
oylr, u) _ „ob (r, zu LES 2 oL,(r, zu BR nn oLun—ı(r, 2m) 
02 ai 02 02 02 


errechnet man mit (A 13) und (A 15): 


(n)(, (n) = 
FE) _ ade Ziele)... (A) 


Damit lautet die N um(r, 27): 


ur ee N. ae) 


o=1 


2) BI On HA 188). 


Zur Bildung des Quadrates von (A 18) multiplizieren wir zwei beliebige Glieder der rechten Seite 
miteinander. Die beiden auftretenden Typen von Produkten sind: 


Rang g2c+2 


2 y2 2 y2 
und Fir rt Mara +(?) Bea MAG), 
r 
Wegen (48c), (48b) und fis, +1?) Barrılr?) = Yo ton .2(r?) lassen sich die Ausdrücke (A 19) 
als quadratische Formen 94, +2(r?) schreiben. 


(n 


Damit gilt für [u®(r, 267) die Entwicklung: 
[uw(r, 17)]? En B3 VE ee rt (A 20). 
0 


(RI 


Die Formeln (A 14) und (A 16) führen leicht zu 


02 yW(r, 3 
02? 


Ne ke ie sure (ABI), 


A4. Die Entwicklung von |wXr, z{7)) |” nach Potenzen von e: 
Mit (A 6) erhalten wir 


Tau ul) = le) 
r PS re] Lt Bar) dr An 
7 / 
Wegen (45a) und (A 5) wird hieraus: 
Ai Een em) 2 
= - Vy (r?) 2 ‚ao-a;tr?) e* )° 
>= 
Umordnung der Doppelsumme liefert: 
x E +o 
(n) 
ea.” PL TOT) ee) 
1 o=1 


Wir betrachten einen einzelnen Summanden der inneren Summe der rechten Seite: 


2;(r?) Q40—2;(r?) 
24 
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Dieser Ausdruck ist eine Summe von Gliedern mit drei Faktoren @; (r?). Von diesen drei Faktoren 
hat einer bestimmt einen von Null verschiedenen Index. Die anderen beiden Indizes können 
Nullsein. Nach Hilfssatz 3 läßt sich ein solches Produkt als Linearform der @,, ®g, ..., @o+1 
darstellen, wobei w, nicht auftritt. Setzen wir 


l+o 
D 340) guo_2;0) = Korıln) 
so haben wir für (A 22) die ee 
: Sal. =) - Du, EEE ER LEN; 


Die folgenden Entwicklungen werden wieder nur kurz skizziert. Man muß zur Ableitung den 
gleichen Gedankengang wie soeben einschlagen. Es en für9 22: 


1 9Q@(r, Bun a 2w = 0 
r .— r = Diele a >) 
58) liefert in Verbindung mit (A 24): 
1 oLuu_alr, 2 1 9L4.—.(r, 0 \ E 
lt zu) _ ET. (r, 0) I a) g2° W223. {Ab}. 
7 or 1 or Zi! | © 
Dabei wurde gesetzt 
1 9L4,_4,0) 1 92@(r,0) ,(r?) 
$ en — 2 N4u—4,v B; m Z, O4u—4,v : - (A 26). 


Wegen (45a)läßt sich (A 26) offenbar schreiben: 
1 la, (?, 0) 
r or 


=, leer 5 Nur ae az 


Nach Induktionsvoraussetzung geht nun L4„—4(r, 0) mit r mindestens wie r! nach Null. Damit 


verschwindet h a 0) 
r r 


mit r mindestens wie r?. Setzen wir in (A27) für r den Wert Null 


ein, so erkennt man, daß der Koeffizient von w,in l4„—s verschwinden muß. Wir können also 
für (A27) genauer eine Linearform l4„—s(r?) der @(r?), ..., ®4u—s(r?) schreiben. Damit 
nimmt (A 25) die Form an: 


EN 
en 2 = I herun-aln) g20 PER A ei 


Zur Aa der Entwicklung von 


a ENGE) 


(n) 
nl ZN —_ Ye) (n) gen KING eu e® oL,(r, zu) I % e 
wer, zu) = (r, zu) = 5; ar By r or | r Or © 
benutzen wir a 23) und (A 28). Das Ergebnis lautet: 
wr(r, ip) == B3 io +1(r?) gzeH+l RE NR EIER ge SC (A 29), 
o=1 ; 
OFEN AR ea Ta 


Zur Berechnung von [w®(r, 247 )]® multiplizieren wir zwei beliebige Glieder der rechten Seite 
von (A 29) miteinander: 
I4o, ar (7?) lie, Sr ı(r?) Br N 
Wegen 
1 1 = 
bio, +10?) la, + 1(7?) = Fate, +09 +2(r?) 
erhalten wir 


oo 


[wm (r,.20 )] = gas Her) rare re a: 


A5. Die Entwieklung von [u (r, 1)]’ + [we (r, m): 
Addition von (A 20) und (A 31) liefert: 


[uote IP + Lumen SF attaad)errtt 2. . (A32) 
Eingegangen am 6. November 1951. 2 


”- a eh ee ETTE a DEE 2 Zu. vr 
re . Re ! 3 u “ 4 
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Beiträge zur Theorie der achsensymmetrisch belasteten 
schweren dicken Kreisplatte 


Von Istvan Szabö in Berlin 


.. Auf einige Resultate früherer Untersuchungen zurückgreifend, wird noch einmal das Problem der 
achsensymmetrisch belasteten dicken Kreisplatte unter ewakter Berücksichtigung der Massenkräfte behandelt, 
wobei die zwischen Platte und Unterlage auftretenden Schubspannungen als Reibungskräfte im Sinne des 
Coulombschen Gesetzes erfaßt werden. 


Taking up again some results of previous investigations, the author discusses the problem of the thick 
eircular plate with awially symmetrical load. T’he own load is exactly taken into account, and the shearing 
stresses between plate and base are considered as friction forces in the meaning of Coulomb’s law. 


Revenant & quelques rösultats de recherches anterieures, l’auteur discute encore une fois le probleme de la 
plaque circulaire, Epaisse, chargee symmetriquement a l’aze, en considerant ewactement les forces de la masse 
et regardant les tensions poussant entre la plaque et la base comme forces de friction agissant au sens de la loi 
de Coulomb. 


B cBasu cC HEeROTOPEIMH pesy.IBTATaMU IIPesKHUX HCcIeNOBAHHi ee pas PaccMmaTpuBaertcH 
IPoÖ,IeMa OCecHMMeTpryecku HATPysKeHHOH TOICTOH KPYTOBOK INIACTUHLI IPM TOYHOM yuere 
MACCOBBIX CHJ; IIPH 9TOM, CABHTAIMIUME YCHJIMA, BO3HHKAMILIME MESKIY ILIACTUHOH HU OCHOBAHHEM, 
YYHTBIBAIWTCH KAK CHIBI TPeHMHA B CMbIc/Te 3aK0Ha Ky.ıoHa. 


1. Einleitung 


Das Problem der auf dem elastischen Halbraum aufliegenden oder anders gelagerten und 
achsensymmetrisch belasteten dicken Kreisplatte ist vom Verfasser unter Vernachlässigung des 
Eigengewichtes erschöpfend behandelt worden !). Die übliche Berücksichtigung des Eigen- 
gewichtes in der Weise, daß man es als gleichmäßig verteilt zur äußeren Last schlägt, dürfte 
bei dünneren Platten eine ausreichende, sicherlich aber eine sehr rohe Näherung bei Fundament- 
platten sein. Die Schwierigkeit bei der exakten Berücksichtigung der Massenkraft ist — wie 
bei den meisten physikalischen Problemen —eine mathematische und besteht darin, daßes allein 
mit Hilfe von Bessel- und Exponentialfunktionen — aus denen man üblicherweise die Lösungen 
aufzubauen trachtet — nicht möglich ist, die entsprechenden Differentialgleichungen zu erfüllen. 
Der Gedanke, weitere Funktionen, insbesondere geeignete Polynome als Partikularlösungen 
heranzuziehen, ist naheliegend und erweist sich auch hier als fruchtbar. 


2. Problemstellung für die freigelagerte auf dem Halbraum aufliegende Platte 


Die auf dem als elastischen Halbraum angesehenen Boden von dem Schubmodul @,z und 
der Poissonschen Zahl», aufliegende Kreisplatte habe die elastischen Konstanten @ und », 
den Radius «a, die Höhe h, das spezifische Gewicht y und werde durch die achsensymmetrische 
Druckspannung f(r) belastet (Bild 1). Setzen wir zwischen Platte und Boden stetige Berührung 
mit der Haftreibungszahl «, voraus, und bezeichnen die aus dem Bild 1 ersichtlichen Ver- 
schiebungen und Spannungen für den Boden mit dem Index B, so lauten die Randbedingungen 
folgendermaßen: 


nd en: SEE ON ee (D) 
een Salt) een. (2) 
ee ars re De er ee) 

tur. 0% 
„sh bw. 20: ist ar ae... 
eebzw er = Oimnda Ver Sa’ ist, we WE, ./.. (5) 
} {ul Zr Sa 
N chzw.. 2 = VO lstta > “ en en Ne RE N) 
ET RN ir 0 FR Pk ar N En 0) 
N N ae ee  (-) } 
Nach Einführung der Gesamtdehnung 
u, M ow 
= — ER ne (9 
= r ee or 02 9) 


1) Ing.-Arch. Bd.19 (1951), 8.128 und 8.342. 
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liefern die Ho o keschen Gesetze zwischen den Spannungen und Verschiebungen folgende Be- 
ziehungen: 


ow ‚O ' ou ‚Oo 
ne 2G MR = 2@G— er 
2 + Ta A ei ee 
/ LER ; , EN 2 RO 
u vo ou  0Ow | 
9 er 2 7 ee s u FITe re 
a «| rn" Be: i | 2 a) 
während sich aus den Gleichgewichtsbedingungen nach Bild 1 die Differentialgleichungen 
OB Ort Ar 0%, ,—,, 9 
a Fr RE. 
ergeben. 
al]. UM 
N! N IM _, 
YZ 


I- 


PS ISZIS DSDS DZ I SIDSSZREL ISIISISSHSTSST 
I Na 


! 
ae | 
| | 
'z, | | 
BR 2 j | 
4 3) | | 
pP: 79 | 
8 BET | | | 
| | j 
Ie® | 
| Bu U 
1 
| ö6, 
g rd 


Bild 1 


Setzen wir die Spannungen aus (10) in (11) ein, so bekommen wir nach Einführung des 
2 0? 


Laplaceschen Operators A = -H ar die Differentialgleichungen 


02? | r dr 


09 , y(l—2») RN 

IS 769} ’ 1 — — _— . . . . . D . . . . 

(1— 2») Aw- 2, . G 0 (12) 
u oo ch 

1—2»)| Au— — —= 3 A Sr a HE rn 

( (An n! r 0 (13), 
aus denen wiederum 

AG =.I: 3 Wr er a N AT 


folgt. Diese Differentialgleichungen gelten — wenn man y = 0 setzt — auch für den Boden 
und das Problem besteht nun darin, für diese Differentialgleichungen solche Lösungen zu er- 
mitteln, mit denen man den Randbedingungen (1) bis (8) gerecht werden kann. j 


U, a ke oe RT Di ER TEEN 
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3. Ermittlung der Verschiebungsfunktionen für Boden (Halbraum) und Platte 


Nach meiner ersten Mitteilung?) sind die Verschiebungen des Bodens, d.h. die den Dif- 
ferentialgleichungen (12) — für y= 0 — und (13) genügenden Funktionen 


ya 21 Dd(} 


oo ._ 


un- | 
6 


Hierbei bedeuten ®(A) und 'F(}) beliebige Funktionen und J, (Ar) und J, (Ar) die Besselschen 
Funktionen nullter und erster Ordnung. 

Die für die Platte gültigen Verschiebungen wollen wir aus zwei Anteilen zusammensetzen. 
Der erste Anteil besteht aus den für y = 0 gültigen Lösungen der Differentialgleichungen (12) 
und (13) und diese sind?): 


EN uN I. ee (15), 


oM 
F)tzg rt — 41, —A2ı) 


Jı(Ar)e 2 di . (16). 


=) IP nd ] jo Sreko]is er] oAr)+Az+B... (17) 
||? An N 2 | 

Uo | Pagani ’ ) e r- | | (18). 
| | Yenepnplunde Er ran) 


Hierbei sind A, B, P,Q, p, qund A beliebige Konstanten. Der zweite Anteil soll aus den Lösungen 
der Differentialgleichungen (12) und (13) bzw. (14) bestehen. Wie schon einleitend betont wurde, 
suchen wir Polynomlösungen und setzen aus diesem Grunde als Lösung von (14) 


PER AR ee ee el er LO) 
an, die in’ (12) eingeführt mit der Abkürzung 

A, 7 ee) Y 9 

NEE N vn) 

die Differentialgleichung 

Rasa 2 Ro SaegeL ag as Wele itelen . (21) 
ergibt. Gehen wir.mit dem Bene 

Bus ES" DE SUN BESSERE ee (22) 

k=0 n=1 


in (21) hinein, so bekommen wir 
Sk kN) + Nb,n(n—1)r"72 + Db,n?+20=0 


woraus sich nach Koeffizientenvergleich 4, = 0 für k > 2, b,—=0, b, = O für n> 2 und schließlich 


EEE aan re (28) 
ergibt, so daß gemäß dem Ansatz (22) 
ven tum: — (+26)? —-br v2... (24) 
ist. Setzen wir (19) in (13) ein, so erhalten wir 
u 
Au — ER. 0, 


aus der mit dem Produktansatz u= R(r) Z(z) für den Parameter Null die Differentialgleichungen 
BEE rennen. . (25) 

N NER Eu .. (26) 

hervorgehen. Nun ist (25) eine Eulersche Differentialgleichung mit dem für r = 0 regulären 


Integral R(r) =r, während das allgemeine Integral von (26) Z(2)= &+- 0,2 ist, so daß mit den 
beliebigen Konstanten a, und a, win der Gestalt 


is (27) 


2) Ing.-Arch, Bd, 19 (1951), S. 128, 


Yen 
TEEN 
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“erscheint. Das Bestehen der Gl. (9) ergibt mit (19), (24) und (27) | | 2 ar x E SR 


R, AA ee OO 5 

e: woraus nach Koeffizientenvergleich | s 

ü A= 2 -(CH2B)W 0 Ara . (28) 
Sa folgt, so daß für die Gesamtdehnung ER | 3 
Br - = 20 +m+2[m—(CH2b)]e 2.22.0202.) 


zu setzen ist. Nun wollen wir die Lösungen (24), (27) bzw. (29) so spezialisieren, daß sie den 
Randbedingungen (1), (3); (8) und der Randbedingung 2) für f(r) = 0 genügen. Berechnen wir 
die Spannungen, indem wir mit (24), (27) bzw. (29) in die entsprechenden Formeln von (10) I 
hineingehen, so liefert (1) 


t+2,=0 ..2. 000... see 0, 
womit auch (3) erfüllt ist. Aus (8) ergeben sich die Beziehungen 
v 


a zehn 


v 
Ä ’ “+7 9,4 m 0 RE EEE EN (32), 
a 2° | 
N während (2) mit f(r) = 0 auf i Be 
! | | - 
Be 3, 4,=0.: RR DE er BR 


führt. Die Auflösung von (28) und (30) bis (33) liefert mit Rücksicht auf (00) und a 


4; = a ar Te >% | = > 
40, = 0, Sage: “ 46(1+4»)’. ee von > 
so daß aus (24) und (27) die Lösungen LEE RA 
IT we | no 


hervorgehen, die bis auf die Übergangsbedingüngen R bis (7) Ah Randbedingungen ‚genügen. NIE 


Durch Superposition von (17) und (34) bzw. (18) und- an) wir die Sa Se ne 
‘ falls wir für a,+ B wieder B Schr in folgender Form: ER 


x. Be Pas, 
> w= | a e T 


K- 


a on a FB 
| 


er Auer 9 DA SER ER a) 4: ANSE in 


» 


. 
4. Erfüllung der ersten drei Rank 


Berechnet man aus (36) und (37) x gemäß in so 
 Jı(Aa) = 0, so daß sich für A — falls man mit i;die unenc 
BT ıd) BezeichzeER — als Zwei 3Y 


E- k=ı (39) 
; a4 
h a2 en EL SER I jr (A,r)+ Az+B— Teaznet+ vr?) 
F : 
& 2 eh [ 4 Pr BI 
| ee En 
| ee a ae ner 


Die Randbedingung (2) ergibt für die mit (39) und (40) nach m) berechnete Vertikalspannung 
die abe 


3 Be! 2 aD 5 (Pr + 1) — Ar Pr + 320.) A ee ag pri a Fl) 
re die mit rdr bzw. Ju(Aur) rdr multipliziert und zwischen r = 0 und r = a integriert wegen 
B, a a Bes L 
f Taerdr — 0, [ Elan Aldr)dr— | Ra 
: Y a 2 y L 0 für. %k em 
VER N ZU 
I REN na" Fr) dr ee a ae (AF) 
al hzys | 
% Br Ani DR R Aa ' _ nur 
a en. een, ERW (42) 
Be Kühe | 
RT H a8 er Nie Randbedingung &) liefert . 
RER EN = RE Kar ; nn E Be, ED ELETET . . (43). 


Een wir (a2) und an nach 27 2uN Q, auf und setzen diese Werte in (83) und (40) ein, so be- 
kommen wir RN f 


Y u . 


[e- rd Ahr + e’r? | + 


= 22 ze = a a (aa) 


el, 


ns 
is 


AR an + Art Bo gg Far EIG + vr?) 


57 A $ 1) 


a ii 
’ ri, * 


1 24) ehe = er | + 


(3. ragen ten I; ih a ER (45). 


tal ZIEcH 


Der. 
et 


De 
Fi 


“u. 
e 


5. Die Erfüllung der Randbedingungen zwischen Boden GE: und Platte Y R 
Die Übergangsbedingung (4) ergibt mit (15) und (16) bzw. (44) und (45) die Gleichung ' 


ERRELL N TIE 


296, I TeTL | % : | FE 
j . ö s “ , o 
FR | el I 5 (ae — Sin Ah) — 6). Er 
i a _ 0, (Aert h — Ein Ah) + Col di ‚| dien —yh er 3 
’ or üürrsaol f 


Mit Hilfe der Hankelschen Umkehrformel für n = 0 
[p(2) Iu(Ar) AdA = y(r) ; e(A) = /y MAlrrd...... (a7J: a 
0 i ' k Er 


und der Integralformeln 


| Jar) rdr = T Ilka), | Klar)Ilr)rdr enaarn V{ASyH 
ö .Ö 
erhalten wir aus (46) j ar 
Ya + ,,00=—g, 2a na) I a [7 
B. | | aa | 
\ on | ; A G54 ; a VE 2 


wobei zur Abkürzung 


a a -Einah + acafan CE 


eingeführt wurde. Br Ba 2 EN ; 
Die ne (5) liefert mit (15) und (4) die Gleichung 


RE: m [4 +24, Net erh 2 +. w Ra a ? 


[wo Ir) ar — 


_ die mit rdr bzw. J,(A,r) rdr multipliziert und von En bis I = inte 


N) . Ar F Ir is . 
zu der Beziehung Pur Er N 


> it 


ij vo Ale) u Eu 2 


[ RP (2) Jı(Aa) a 
ku, 


KR Le Ra TEN 
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führt, wobei 


j ? } a gs at 9 ent Ay,h 
Br seen) } ma —Ay Tr 2Mh)e nie] a 
3 -F A| [8 — 4» — 2A, h) ein’ + eh] + 40, (1 — 2») Sin A, | 
bedeutet. 
a Mit der Abkürzung 
A EÄREN BR EEE 
ET ae Ice ert]+ 
(55) 
+ 9 [(2Ah— 1)etr — ehr] — 40, (1 — 2») Sin ih! 
liefert die Randbedingung (6) 
= a = 
—_ fü s 
20 | voten en AJı(Ar) dA= _|e G>,M J (Ar) ur r } (56) 
0 AlUER a! 
woraus mit (47) — für n=1 — und (48) 
ve el - LLUDEND, 23 
72 Y()+ a)‘ 2 Te ut (57) 


folgt. 
Durch Elimination von ® (2) aus (49) und (57) bekommen wir 
| d Iula Tu(a 
pe rest a DD Ah } u ee) AS) Cart u a. 


(1—9)4 Bet — 08) 


ak SEEN 


Multiplikation dieser Gleichung mit Ida)“ on Integration zwischen A= 0 und A= © 
liefert unter Berücksichtigung von (52): 


pe: 100)", -.- Br MB— je Mi a | sr ei 


J L ’ 
283 ARa)dı 
SER — er a EN .: (A) 124 (50), 
ER t x x 
BER F Ri a0) d di, Be A yh\ (Itda)da 
4 | Man Sicht es leicht ein, daß die Ki vorkommenden uneigentlichen EN existieren. 
ix ae I Multiplizieren wir die Gl. (58) mi wi und integrieren zwischen 7—0 und A= 0, 


er so erhalten wir unter Heranziehung von (53): 


a e 30 a Re Sradliar 2 a 96, AL Sn 1 Dar (A: Ir a 


ee = a 


en 
en nn  uneigentlichen Integrale sind konvergent. 
me jetzt vhs Ranı nee (7) und bekommen aus ihr unter Beachtung von 


(60). 


„5 a AR 


Gag Sa a] rei). 


m], 
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Multiplizieren wir (61) mit r J,(}, r) dr bzw. r dr und integrieren zwischen r=0 und r=a, so be- 
kommen wir wegen der Orthogonalität der rechts stehenden Funktionen 


3 nr KR TRATEN 
k=1 
bzw. 
5 1—» yh 
2 Am Ph Te ee N E 
wobei 
On = | rJo (Auf) Jı(Aır) dr; 0, = [r Jılkar) REDEN 563, 
ö ö 


eingeführt wurden. 

Denken wir nun in (60) und (62) die durch (55), (53) und (50) definierten Abkürzungen A,, 
B, und C, eingesetzt und die so erhaltenen Gleichungen nach den Unbekannten 9, und q;. ge- 
ordnet, so kommen wir offenbar zu folgendem System von unendlich vielen linearen Gleichungen: 


2 oo x 09 
r Be Ay Pk ir > byr %.: di, ’ P2 X Pk je 2 Pr k Ik = :Yn » (n 7 t, 2, ER .) (64), 
=1 k=1 k=1 k=1 


wobei a7 Dr Kurs Par, 4, und y, feste Werte haben; wir wollen sie vorerst explizit nicht hin- 
schreiben, da wir folgendes feststellen: Die Auflösbarkeit von (64) vorausgesetzt, können wir 
aus (58) (2), dann aus (57) D(}) und schließlich aus (59) B bestimmen und damit sind die 
Verschiebungsfunktionen und aus ihnen die Spannungen für Boden und Platte berechenbar; 
d.h. wir haben mit den von uns verwendeten Fynktionen keine Möglichkeit mehr, die Rand- 


bedingung (8) zu erfüllen: Diese Funktionen sind eben nicht allgemein genug, um mit ihrer _ 


Hilfe einem solchen komplizierten Randwertproblem gerecht werden zu können. Die exakte 
Erfüllung der Randbedingung ist also nicht möglich, wir können sie aber im Sinne des de Saint- 
Venantschen Prinzips befriedigen und das wollen wir jetzt erledigen. 


6. Radialspannung und Radialmoment am Plattenmantel und die näherungsweise Erfüllung 
der Randbedingung (8) 
Am Plattenmantel (r = a) herrscht eine mittlere Restspannung 


h 
IL. La Be = Jul) a) 2yA 
,=7 [» (a,2) de = o| > 4A,— Ich + EEE, 


0 


a 


die mit dem Radialmoment 
h 


' a ) 

Mr E33 [o,(a, z) ger °,] 2d2 = Gm >! 2(1 = a Pı: (A, hre=ikh+ 2 2 Ar h Eoj Ar h— | 
Eu a PER, 65 
— /,h Sin A,h— 2 ©in A,h)+ qu(Ar hre?s" —2 4, hof A,h —A,h Sind,h + fe 9 

170 

+2 Ein A, h) — Cr [Ah Sin d.h +2 (1 — 60f A,h)] 2 (1 —2v De u | 
zusammen das statische Äquivalent der radialen Restspannungen bestimmen. ER, wir also 
den Plattenmantel durch — 0, bzw. — Mr, so wird dieser im Mittel kräftefrei sein. Die zuge- 


hörigen Verschiebungen sind ?) 


a 


( E 
) 3 Mr 2, (— 2) + (1 —»)r? re, — )\ (66) 


1 
G(1l-+ v)h?) f 


a8 ee h TB | 
= re N a 


und nach (9) ergeben sich die Zusatzspannungen 


0, Me AUSB EN, Pi (e- 5) EIERN 


s) Siehe Fußnote auf 8,359, 
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Aus den letzten Gleichungen folgt, daß von den Randbedingungen (1) bis (7) lediglich die aus der 
Randbedingung (5) fließenden Gleichungen (52) und (53) geändert werden. Aus wg = w-+ w 
fürz= h bzw. 2, —=0 ergibtsich an Stelle von (52) bzw. (53): 


[vana) 51 4An | a neepile | y’e) } 
| 


(69) 
Rn u: Perl) 
GA +») Gi+tnB 
bzw. AP) Ra)dr 1,(A,a) B IL, al EEG (70) 
j PEpey = aJ,(A,0) B, NETT ): 


0 


Dadurch treten jetzt entsprechend (69) und (70) an Stelle von (59) und (60) die neuen Glei- 
chungen 


FR y: RR v0,h 3M:b#®+(il—») a] 
ur, ea, | 2 @) 4 26 +) " 264 + »)h 
ER Ca)ydı 
ne ArAr Jo (Ar „Aue 2 (71) 
TE 0) 
T(da)da (1—» yh\ ( I2(Aa) di] 
—4(1— 5) 2% Jo (A: JE 2 — 2 ! er el: 1 Ar | 
er 0 0] 
ee ln). 
Jo @)Bı 26 1H9)% Gl v)AR | 
RS DAT il Ala)d 
N) ER) Fark ER 
+3G, >> (22 — A) (22 — 14) (73) 
| x f  32J2(Aa)dA eg yh\ (R(aa)dA |) 
Pa ee) FREE f 


0 
Fürn=1,2,3...liefern (62) bzw. (62a) und (72) das System der unendlich vielen linearen 
RR für die Koeffizienten p, und g;. Damit ist das gestellte Problem für die ersten 
sieben Randbedingungen exakt und für die achte im Sinne des de Saint-Venantschen 
Prinzips gelöst. 


7. Das System der unendlich vielen linearen Gleichungen für die Koeffizienten p; und q; 


Die beiden Gleichungen (62) und (62a) können wir zu einer einzigen zusammenfassen, 
wenn wir zu den positiven Nullstellen , >0(k=1,2,...) noch t,— 0 hinzunehmen (d.h. 
definieren) und 


einführen. Dann lassen sich (62) und (62a) unter Beachtung von (50) und (55) in folgender Form 
schreiben: 


»2 [PD (Ar) pr + PU) (Ar) al —Y® (Au)pn — 99 An) u = 


5 . (74) 
=D’ 49,9% (1— 29) SinAyh+2(1— 2») a2 16, JE (Ana) Cof ih (n=0,1,2...) 
k=1 
Hierbei bedeuten: | x 
| (A) = An l2 A —1)e + er"), | (75) 


YO (A) = a2 I8(Ana) (Anke ""— Sina) | 


EHRT Fl ZN ya Bl ale BE ER ae Fa En 
t \ “ wn N Se 
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Aus (72) erhalten wir, nachdem M, nach (65), B, nach (54) und schließlich für (71) und (72) 


RO) (I)de 22 (a)da F BR(i)- de 
ae ae a er ea 
h) f) 0, l) B (76) 
20) [AO _ am, AR | RE a 
AR—A) 6(t?- in) e (A? —Ar)(A?— An) (?— 1) (?—t7) 1% 
f) v 0 D j) 
eingeführt werden®): 
S (2) (2) (2) (2) J a 5 
Dt Mn (41) Pk un Py (> 2%) 4] ap (A,)Pa hr (A) On age: 2 (12 ya Sind, h + 
Peg 
Rs Dt a2), IH) un Sina +2 Co) Arh)]— 
eg! | Am) AnArh? (77) 


— bi al, o(Ar a)[(1 —2 v5) ar I» Sin, h+2 2 (1 — 5) Ir „&o] 4; A + 


G, kn | 
| en Be Er De 3(1 —v)4|, (n=1,2,3..,) 


Hierbei wurden gesetzt: 


5 6(1—v)J (A a) My? . 
Pr 4) = Een [Ant 4 2744 Co] A,h— Ash Sin 7, —2 SinAuh] — 
ne) 
G(l a 7 iA; a) a? I Y) (2 Ak h—1) Een erh 188 (78). 


+2 (1— 95) Jo (Ara) a In (A,he *" — Sin}, h), 


Jo (Aa) 


n 


yala)—: Na, aa near em (79). 


Durch die Transformation 
%»=2(1—2») CE 0R: U RZ EL Ce (80) 
gehen aus (74) und (77) folgende Beziehungen hervor: / or 


2 [PD (Ar) dx + PD (— 44) er) — Y (An) On — y) (— An) en = 
ur; 2 = 2 (1 Br v) 6% [2 On Sin 1, h— gi (/.)] + h (81) 
= = 


+2 a — 29), [y (A) + at Jo (An a) &ofA,h]), (n=0,1,2...) E 
= [9 (A). + PA) va) a I AR, 


RE En Ro & 
Do BE ne a [4,h ©inA,h+2(1 A ü . 


@ 
=, Fra) [A —2v5) aA) SinAeh+2 (1—vn) Ten CojAr h]— pp) "a4 


kn 
AED 
Arm 


Bezüglich der Lösbarkeit der Gleichungssysteme (81) und (82) mit den Unbekaunten dp 2 a Te gr 
und &, &, ... läßt sich folgendes aussagen: Wir dividieren: (81) und (82) durch y® (A) bzw i 
ya(A,) (n= 3 2,...). Dann gehen die Konstanten (d.h. die mit d, und. nicht Beh eE 4 


#) Bezüglich der allgemeinen und numerischen Auswertung dieser Integrale 8. die in Fußn 
S. 359 angegebenen Arbeiten, F 


4 


+27, a Pa Bn 2u—mAl, I Br JL: | 


* 


A Te 1 Re re a 
Ey IRUORT 
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Glieder von einem gewissen Zeilenindex an gegen Null. Von diesem Zeilenindex N an gehen 
für beide Systeme die Koeffizienten der ö, und eg; bis zu einem gewissen Spaltenindex K< N 
gegen Null. Folglich erhält man von der N-ten Zeile ab ein homogenes Gleichungssystem für 


ÖR, Öx+1, :.. und ex Exrtı »:.., welches durch öxk = ör+ı= E45 Eat An 0) 
befriedigt wird, d.h. für großes k können 

Pi: —— 2 (1 — 2») Ck3 I: == 0) 2 Sure ® E - > 5; E A € % R (83) 
gesetzt werden. Die ersten 9, Pa .:-.-; »K—-1, 91 I»: - » 9x—ı Werden in üblicher Weise aus 


den ersten (K — 1) Gleichungen vom System (74) und den ersten (K — 1) Gleichungen vom 
System (77) gewonnen. 


8. Die eingespannte Platte 


Die gesamte Last 
[2 
— 2n/flr)rdr + nahy = o,n(a — a)). BRENNER 8A) 
0 


werde bei z—= h auf den Kreisring 
rz(a?— at) gleichmäßig übertragen 


(Abb. 2). Die Randbedingungen G INN fir) 
In mm 


nee ln. tl (85), 2 
OA fr) Aa rr (8b), 

EEE ER l/h h 

az) | 

I er EEE 

TEEN ey. 290), ll | | 
Eat 3% (91). DR REN zZ 
Zunächst ermitteln wir eine den Dif- RB 


ferentialgleichungen (12)bis (14) genü- 
gende und die Randbedingungen (85) bis (87) (letztere für f(r) = 0, und (89) bis (91) Be [rie- 
digende Polynomlösung. Aus dem Ansatz (19) erhalten wir: 


Nr (ran) y 
Er aaenE + . . . . . . . . . . . . . (92) 
und 
+2 2NY P s 
Sy). ee a und (98). 
VE 


Den Differentialgleichungen (12) — (14) für y = 0 genügende und die Randbedingungen (85) — 
(87) bzw. (89) und (91) befriedigende Lösungen sind aus (44) und (45) sofort zu entnehmen: 
Man streiche die letzten Glieder und ziehe in (44) w(a, h) ab. Superponieren wir zu diesen Lö- 
sungen (93), so erhalten wir 


e- >, ei (3 —4» + 2iu2)e + e’*” ] E= 


TEL ES —2») [8 — 49 — 23,2)’ + ne Sul haar 4 


1—2»)y RN (94). 
Asien Bi. ran ++ 


en ie 4v—2A.h)e%" + + % Sim ande An + 
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EN > | Pk 3 — Ayz 2°] 


RM OD 
RI. 1 War leben > BAER 9 Ap2 RR | HR 
24781 29,) [3 49 + 2)22)e +e ] Ar er ıl un i 
Diese Lösungen sind nun noch den Randbedingungen (88) und (90) anzupassen. DO | 
Aus (88) ergibt sich mit (94) und (95): » 
be SU Ber ar A > ET LIOM hin Art 
2 it Son enihe Sin Ah) Be Sin A;h)+ 
| | 0 fürr<a (96). 
+aCofAhhlani—yh=ı L Be 
| x (a? — a?) RE z 


Multiplizieren wir diese Gleichung mit Jy(A,r)rdr und integrieren zwischen r=0 und r=a 
so erhalten wir wegen der Orthogonalität der Eigenfunktionen: 


) 
De De ne Ne BSBAREE LE. Ah Si In e 

30 —3,) (A, he SinA;h) DA) abe Sin A,h) +6,C0] A. h | 5 £ 
& La,Jı(2a,) ne : 

zarr,(® —a)@Is (Ara) * d f 


Aus der Randbedingung (90) kommen wir zu der Beziehung 
Pr (Arhe *" 1 Sin Ah) + qu(A,he*" + SinA,h)= 2(1 —2v),©inArh. . . (98). ° 


Die unbekannten p, und g; lassen sich aus (97) und (98) ermitteln: 


1—2» Aal h ß | 
= Sin a IE a (Ayh + Sin 24, h + Sin? h) — (Ach er" ı Ein, n| h | h 
RE (99). 


B ' 
6) E : 2 Br Eh 
de = Sea MR | 6% (A,he % —- EinA,h) — G, (At Sin 2%),.h — Sin?) ) 


Der Fall der schneideförmigen Stützung auf der unteren Plattenseite, d.h. a,>a ergibt aus (97) 
mit der L’Hospitalschen Regel | 
cd Ee 


a RR LONE 


p 
k 
; 
In diesem Falle hat man mit Rücksicht auf die Konvergenz der auftretenden Reihen zu beachten, F 
daß für «,> adie linke Seite von (96) Null zu setzen ist. 


9. Die freigelagerte Platte 
Die Randbedingungen (85) — (90) bleiben bestehen, während (91) durch 


U A m ER TE A E e (101) 


zu ersetzen ist. | 


Aus (44) und (45) sind die den Randbedingungen (85) bis (87) und (89) genügenden Lösungen: 
© i 
uv— ern u, NR 9 — AL? Ay? 
Da 9% [8 — 49 + 29,2)e Tender 


+ 418 —49 242) + E%] 4 (1— 29) Sin Az} Jo(Ar) + 
HAl@— Rh) 


ia ee (102) 


2 


21-8 —Ar +2) TH IB — Av —2uh er er] + 
Fe | 


+40 (1 — 2») &inArh} an 
Be 
u= I {mE 2) + + [B—4r + 2A Dee 2 
a (108 Ya 
AA | 1 IV RErRER 
Sc 0 DL ee F au) zr | | 


a NT Air 
N ee . 
“. 4 Ss n 
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Die Randbedingung (88) führt auch in diesem Falle zu der Beziehung (96), aus der sich wiederum 
(97) ergibt, während die Randbedingung (90) ebenfalls (98) ergibt, wodurch p, und q; berechnet 
werden können, ohne daß man noch eine Möglichkeit hätte, der Forderung (101) zu genügen. 
Diese Randbedingung kann auch hier nur im Sinne des de Saint-Venantschen Prinzips 
befriedigt werden. Zu diesem Zweck berechnet‘ man das statische Äquivalent der Rest- 
spannungen o, (a, 2) am Zylindermantel und erhält dafür wieder die nach (65) gegebene Einzel- 
kraft bzw. das durch (66) gegebene Kräftepaar. Belastet man den Zylindermantel durch —o, 
und — M7, so hat man dort ein Nullsystem. Beachten wir (68), so stellen wir fest, daß lediglich 
die Randbedingung (89) zu berücksichtigen ist bevor man die Verschiebungen durch Super- 
position von (67) und (102) bis (103) hinschreibt. Man erhält: 


w— Sina ar 2ER LEN HB — Ar — 2) te] + 
k=1 


Jo (Ar) 


+ 40,.(1— 2») Sin ),2} DIE 2, HA —h)- Ian [(®—z?) +» (a? — r?)] + 


} Le | ; . 
Ar Ga. 1ER { 3 Mr [2v (2 — h2) + (1 — ») (r? — a?)] — vh? 6,(h—2)} — 


na Ar toner. [a — Ar 2er rer) 
k=/1l 


) ; 2943 © er 
+40,(1 — 2») Sin A,h) I1.1- 2): 


u= I {mp8 Ar —2)e "+ er] 4 [8 —4r + 2ı2)e He ]— 
=1 


k 
Jı (4er) yv 
ER TURNIER, Be 
lea az arm” 


. 


ET HERRN 2. 
San) 12 Mn (ey) + 


Hierbei sind die p,. und g; durch (99), MR und o, durch (66) und (65) und schließlich A durch 
(41) gegeben. 


Eingegangen am 6. Dezember 1951. 


Eine neue Methode zur Berechnung der Strömungsfunktionen 
bei zeitlich veränderlicher Kontur”) 


k Von Guntram v. Gorup in Darmstadt 


\ ) Auf dem Weg über eine Funktionalgleichung werden die Strömungsfunktionen im Außenraum einer 
zeitlich veränderlichen Kontur unter der Voraussetzung einer meromorphen Abbildungsfunktion in ge- 
schlossener Form bestimmt. Als Beispiel für das Verfahren wird ein Joukowski- Profil behandelt. 


Using a functional equation, the flow stream functions for Ihe exterior of a contour that is variable with 
the time are determined in a closed form, whereby it is assumed that Ihe mapping function is meromorphic. 
The Joukowski profile is treated as an example. 


Par la voie d’une &quation fonctionelle les fonctions de courant dans l’espace exierieur d’un contour 
variable en &gard du temps sont determinees a la supposition d’une fonction de figuration meromorphe en 
forme serr&e. Comme exemple pour le prooe procede un profil & la manierede Joukowski est traite. 


Ilyrem, Benyıuum yepes O1HO PyHRIHOHAJIBHOE YP&ABHeHHe, ONLPeNEeIoTch PyHRIUM Teuenns 

BO BHESKHEM IPOCTPAHCTBE H3MEHFIOIETOCH BO CPeMEeHH KOHTYPA, B IIPeALOTOKEHHH MePOMOPPHOH 

. PyHkumm OTOÖpakennn, B BaMEHyToü dopMme. JLyıst 9T0rO crocoÖa PaccMaTpuBaeTchH B KauecTBe 
IpiMmepa, ee „RyKoBcKROro. 


1. Einleitung 


Ein unendlich langer Zylinder mit einer zeitlich veränderlichen Kontur (€ befinde sich in un- 
begrenzter ruhender idealer Flüssigkeit. Infolge der Veränderungen von (, die sich aus starren 
Verschiebungen und Drehungen, sowie aus Deformationen zusammensetzen, entsteht dann im 
Außenraum eine wohlbestimmte ebene Strömung. 


ER *) Ausschnitt aus der von der Technischen Hochschule Darmstadt 1951 genehmigten Dissertation des 
| Verfassers (Referent: Prof. Dr. C. Schmieden, Korreferent: Prof. Dr. A. Walther). 


a. 
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Eine exakte Methode zur Berechnung der hierzu gehörigen Strömungsfunktion wurde erst- 
mals von Steudingt) angegeben. Hiernach wird unter der Voraussetzung der Kenntnis der Ab- 
bildungsfunktion der Kontur die Strömungsfunktion als Potenzreihe bestimmt mit Koeffi- 
zienten, die sich nach Art von Fourier-Koeffizienten berechnen. 

Demgegenüber liefert die vorliegende Arbeit eine ebenfalls exakte aber außerdem direkte 
Methode. — Eine solche Methode liegt für den Spezialfall der starren Kontur-Verschiebung und 
-Drehung bereits in der — einen anderen Weg benutzenden — Arbeit von Morris?) vor, Wesent- 
lich neu ist'also nun die Ausdehnung auf sich deformierende Kontur. 

Der beschrittene Lösungsweg gründet sich auf folgenden Gedankengang: Die Strömungs- 
funktion A, deren zunächst unbestimmt bleibende Zirkulation gleich Null gesetzt werden mag, ist 
durch zwei Bedingungen festgelegt: 

1. Randbedingung; d.h.: die Normalkomponenten der Konturgeschwindigkeit und 
der Strömungsgeschwindigkeit müssen längs ( übereinstimmen. 

2. Regularitätsbedingung; d.h.: die Strömung A darf im Außenraum und auf © 
keine Singularitäten außer einer Senke im Unendlichen besitzen. 

Wir werden die Randbedingungen durch Benutzung eines Kunstgriffes zur analytischen Bil- 
dung konjugiert komplexer Ausdrücke in eine Funktionalgleichung für A umformen. Diejenige 
Lösung dieser Funktionalgleichung, die auch der Regularitätsbedingung genügt, ist dann leicht 
explizit anzugeben. 

Die Entwicklung und Anwendung von Formeln für Kraft und Moment auf die Kontur wird 
in einer später erscheinenden Arbeit durchgeführt werden. 


2. Kontur und Abbildungsfunktion 

Vorgelegt sei in der(physikalischen) komplexen z-Ebene die zeitlich veränderliche Kontur (©, 

von der wir voraussetzen, daß sie sich nicht überschneidet und — abgesehen von höchstens 

endlich vielen Spitzen — glattist. Ferner sei in der Bildebene £ eine feste durchweg glatte Kon- 

tur ©, gegeben, die den Punkt &= w enthalten darf. C', zerlegt die &-Ebene in zwei getrennte Ge- 

biete @,und H,. Bekannt sei schließlich die Abbildungsfunktion zwischen den Konturen (\, und 

C; sie soll so beschaffen sein, daß sie das Gebiet A, der ö-Ebene ausnahmslos konform in das 

Außengebiet 4 der Kontur C in der physikalischen z-Ebene abbildet. Wegen der zeitlichen Ver- 

änderlichkeit von (bei fester Bildkontur (', muß die Abhildungsfunktion selbst auch von der Zeit 
abhängen. Wir bezeichnen sie mit 

A 1 RE EEE nk 


Vorausgesetzt wird hierbei für das folgende: 
1. Auf Grund der Forderung der Konformität der Abbildung 4,<—>H und der An- 
nahme von höchstens endlich vielen Spitzen auf C’muß die Ableitung 
ns 
0, 
in 7, überall endlich und von Null verschieden sein, und darf auf C, nur endlich viele einfache 
Nullstellen haben. Für den Punkt {= (falls er in dem Bereich ®,—=H,-+ O, liegt) sowie für den 
‚Bildpunkt von z= ® ist die Formulierung dieser Bedingung noch zu korrigieren bzw. zu ver- 


schärfen. Wir können diese beiden Ausnahmestellen am einfachsten formal mit erfassen, wenn 
wie allgemein als Voraussetzung schreiben: 


02 regulär und +0 in H,sowie regulär und bis auf endlich viele ein- 


& fache Nullstellen + 0 auf (,. (2). 
Dabei soll sein 
z inder Umgebung aller endlichen Punkte, 
EIN WR f 
z in der Umgebung des Punktesz—=o, A (3) 


analog für E, 


2. z(&, t) und die Zeitableitung But 2(£,t) seien in der ganzen Ebene meromorphe 
Funktionen in der Variablen £. ot 


S.20—22. 


. %) Rosa Morris: Notes on Two-Dimensional Potential-Theory. II. Hydrodynamical Problems on the 
Motion of Cylinders. London, Philos. Mag., VII, Vol. 23 (1937), 8. 757. 


!) Steuding: Eine einfache Berechnung von Strömungsfunktionen. Z. angew. Math. Mech. 16 (1936), 
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Die feste Bildkontur O,, die bisher nur allgemein als ‚‚glatt‘‘ vorausgesetzt wurde, werden 
wir im folgenden von Fall zu Fall spezialisieren. Je nach dem gerade vorliegenden Problem er- 
weist es sich als zweckmäßig als ©, entweder die über co geschlossene reelle Achse (R.A. oder 
den Einheitskreis (E.K.)zu nehmen. Im ersten dieser beiden Spezialfälle benutzen wir für 
die Bildvariable den Buchstaben o, im zweiten 7°). Es ist also in den beiden Fällen: 


120, = RA. die Kurve Im ()=0, 
2. 0,= E.K., die Kurve url: 


Die Bezeichnung der Gebiete G,und H,in den beiden Fällen sind u.a. aus der Skizze (Bild 1) er- 
sichtlich. In dieser ist auch der Zusammenhang der beiden Abbildungsfälle über eine lineare 
Transformation angegeben, Alle in H 

o hergeleiteten Ergebnisse lassen C 

sich damit ohne weiteres auf z um- 
rechnen und umgekehrt. 


6-nicht 
physikalische 
7" Blätter 


(z) 


3. Regularitätsbedingung, 
NR und #’-Funktion 


Die Strömungsfunktion A der 
durch die veränderliche Kontur C 
hervorgerufenen zirkulationsfreien 
Strömung muß im ganzen Außen- zieh) zit) 
bereich B = H-+- C regulär sein mit 
Ausnahme einer logarithmischen 
Singularität mit reellem Vorfaktor 
in 2z=o. (Senke im Unendlichen). 
Folglich gilt dasselbe für das Ver- 
halten von A im Bereich B, der 
&-Ebene, nur liegt dort die gleiche 
logarithmische Singularität im Bild 
vonz=@. 


Das Verhalten der Ableitung 


oA 
A = —- leitet man aus dem von 


MR F 
E ab (wegen Z£ vgl. Formel 3). 


= Orl 
rlo)=- 67 

-1_.;1-1 
olt):iga 


er wird nämlich in ®, regulär sein 
0 Bild 1 

mit Ausnahme eines Poles 1. Ord- 

nung mit reellem Residuum im Bild von z=. Dasselbe gilt für A’, nur kommt noch zusätz- 
lich : Doppelte Nullstelle in & = (falls dieser Punkt zu ®,gehört), mit dem Sonderfall einer nur 
einfachen Nullstelle mit reellem Vorfaktor des Lineargliedes, wenn &= ® Bild von 2= » ist (das 
letztere trifft für die meisten der gewöhnlich verwendeten Abbildungen auf O(,=E.K. zu). 

Eine Funktion, die die für die Strömungsfunktion A formulierte Regularitätsbedingung er- 
füllt, nennen wir fortan eine R-Funktion, ebenso eine Funktion, die die Regularitätsbedingung 
für A’ erfüllt, eine X -Funktion. — Insbesondere ist also Aeine R-Funktion, A’ eine W-Funk- 
tion. 


4. Aufstellung einer Funktionalgleichung für die Randbedingung; Einführung der Operation * 
In diesem Absatz spezialisieren wir uns zunächst zweckmäßig auf den Fall C(,= R.A., 
schreiben also demnach wie in Abs. 2 vereinbart o statt £. 
Randbedingung: Sie lautet: Normalkomponente der Konturgeschwindigkeit z *) längs 


A 
.C = Normalkomponente der Strömungsgeschwindigkeit = 5) Jängs C©. 


3) Nur bei Untersuchungen, in denen — so wie bisher — O\, beliebig ist, schreiben wir Z als Bildvariable. 


„a 4) Man stellt fest, daß bei der Beschreibung der veränderlichen Kontur © durch eine veränderliche 
Abb.-Funktion bei fester Bildkontur ©, die Geschwindigkeit der Punkte von C gerade durch die partielle 


Yo _ Zeitableitung der Abb.-Funktion dargestellt wird. 


5) Der übergesetzte Querstrich bedeutet wie üblich die Bildung des konjugiert Komplexen. 
25 


A 


74 x 4 orup, Kine nous Methode zur Berechnung der Strömungsfunktionen Bd. 32 Nr. ern > 


a nn ee nn _ 


In Gleichungslorm : 


4 2). (2) Ja. 2 ee 


N 

f Dabei soll das eingeklammerte Gleichheitszeichen andeuten, daß Gleichheit nur auf © bzw. C, 
gefordert wird, 

[ (leie hung (4), die die Bildung von Normalkomponenten enthält, schreiben wir um in eine 
Gleichung zwischen Realteilen, indem wir auf beiden Seiten durch den in Normalenrichtung liegen- 


’ den Zeiger (— v#') dividieren : 
| Ne ( = Jim ( ). 
’ PH % 
oder: 
r j = 
i re ai RE N 
ug % 


a 


Operation*, Die in (5) vorkommenden Bildungen des konjugiert Komplexen, durch die 
die analytische Eigenschaft bei den betroffenen Funktionen zerstört wird, lassen sich durch einen 
Kunstgriff umgehen. — Die beiden Seiten von (5) brauchen, als Funktion von a und t betrachtet, 
nur längs (,= PLA. ti be reinstimmen, Kür reelle Werte von o kann man aber das konjugiert Kom- 
plexe einer beliebigen meromorphen Funktion f(o) als analytische Funktion bilden, wenn man nur 
in der Funktion » durch —v ersetzt, das Argument aber beibehält. — Wir bezeichnen die ent- 
stehende neue Kunktion durch Anhängen eines Sternes. Es gilt also: 


Ha) (=) Ha) m /{o)* 1 2 ea a 


Mat en * als einen Operator aullassen, der folgende Eigenschaften hat: 
. Assozintiv zur Addition und Multiplikation. 
z Mit Dilferentialion oder kuragration nach o oder einem reellen Parameter BE NaUuSER 
bar, - | 
>, Involutorisch, 


in 


We 


) 


Kunktionalgleichung: Bei Verwendung des Operators * geht (5) — bei gleichzeitiger 
Multiplikation mit 027° — über in die Funktionalgleichung 


[aa = rt irn] De N ERS REERTN 


(7) gilt gemäß der Herleitung zunächst nur auf R.A. Da aber beide Seiten der Gleichung nur 
noch analyt ische Kunktionen enthalten, so ergibt sich daraus nach dem Identitätssatz die 
Gültigkeit in der ganzen Ebene, Dies ist durch die Verwendung des KeRohzL CHEN Gleichheits- 
ee, zeichens stalt des eingeklammerten zum Ausdruck gebracht. 

Erweiterung auf 0, E.K.: Für den Abbildungsfall O,=E.K. verläuft die Herleitung! 
der entsprechenden Punktionalgleichung nur wenig anders®). Wir können sie uns aber auchleicht 
durch Umrechnung der Gleichung (7) verschaffen, Diese läßt sich wegen der Vertauschbarkeit N . 
der Operation * mit der Differentiation nach o in der Form schreiben &  -- 


y J a" a Y 
N, 94. KA Zen BL Pan ee N ar 
We a a 
(Die Ableitungen nach oa sind hier ausgeschrieben, um Verwechslungen mit A im folgenden 
auftretenden Ableitungen nach 7 zu vermeiden. j {org % AR PR 4 
Nun besteht der Zusammenhang zwischen o und r über eine analytische Paniion u RU 
} o= ol). : Ei RE 
| (ls ist hier nicht nötig diese Funktion speziell zu geben.) Demnach. RN 
0 


7: BEN N 


Da (8) in den Ableitungen nach o linear und homogen ist, hebt sich n 
same Nenner p’ heraus und es erscheinen statt der, aleituingen! nacl h 


®) SE ann Vortragsreforat, Z. ungow, up Mech, an (0), | 


It | h kr DIR: Eee 
% n FERN, d v N . n Bi er 


er a Ey a EN: SER Bey 3 © 
00 por roh ri % 6; Rn 4 . 


a 
I 


an a ee 
Et 
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Um dann die endgültige Form der neuen Funktionalgleichung zu bekommen, ist noch die 
Vertauschung von * mit > durchzuführen. Dazu muß man sich die Wirkung des Operators * 


aufeinein den Variablen rund tgeschriebene Funktion überlegen. — Es muß sein: 


HOMJROESIGE 
Da aber wegen |r| (=) 1 


gilt, so wird 


Dem. wen. lo), 


Damit findet man als Vertauschungsregel 


Are, al 
SE FRtzahge er ztll), 
Die Funktionalgleichung lautet daher in 7 geschrieben: 
N + Nrerid lie RB BUN (12), 


wobei nun, da Verwechslungen nicht mehr möglich sind, der Strich Differentiation nach 7 be- 
zeichnet. 


5. Lösung der Funktionalgleichung unter Beachtung der Regularitätsbedingung 


Wir beschäftigen uns zunächst wieder mit C&,= R.A., gehen also von der Funktional- 
gleichung (7) aus. — Eine triviale — jedoch die Regularitätsbedingung im allgemeinen nicht er- 
füllende — Lösung dieser Gleichung ist 2* 2’, wie ohne weiteres ersichtlich ist. 


Die homogene Funktionalgleichung (mit der rechten Seite Null, also 2 = 0, was eine 
ruhende, somit nur umströmte Kontur bedeutet) wird gelöst durch jeden Ausdruck der Form 


h(o) + h(o)* 
mit beliebigem meromorphem % (0). 


Die allgemeine Lösung der vollständigen Funktionalgleichung (7) ist demnach durch die 
Summe ?) 


BE HR) RO ae en. rar (13) 


gegeben. — h(o) muß nun so bestimmt werden, daß A’ die Regularitätsbedingung aus Abs. 3 er- 
füllt. Dazu betrachten wir die Lösung 2*z’. Sie hat Pole in beiden Halbebenen (also in den Ge- 


a BLM Ä 1 i 
bieten @,und H,)und besitzt im Unendlichen das Verhalten 0 e: Mittels Partialbruchzerlegung 
a spalten wir sie in einen Teil mit Polen nur in @, und einen zweiten mit Polen nur in H, auf. Wir 


fügen dann zu ersterem noch einen Term en ,‚ zum zweiten den entgegengesetzten Term hinzu 


(e = Bild von z=®; u reell, wird weiter unten festgelegt). Die ganze Aufspaltung sei mit 
BER SZ N nr (14) 
u 
0—c 
R(2*2’) — [8(2*z’)]* 
eine W-Funktion (siehe Abs. 3) ist. Und zwar aus folgenden Gründen: 
1. Die Hinzunahme des u-Termes bewirkt, daßin R — S* zusätzlich das Glied 
1 1 ) 1 “ 
BA Di 
( +—)=2r-+40|; 


Gero — 06 


bezeichnet. Durch die formale Hinzunahme des Termes wird erreicht, daß die Funktion 


vorkommt, mit dem durch passende Wahl von u das richtige Verhalten im Unendlichen, nämlich 


o(2). für den ganzen Ausdruck erzwungen werden kann. Es genügt, u rein reell zunehmen; 


?) Wegen der Linearität der Funktionalgleichung gelten alle bekannten Überlagerungssätze für die 
Lösungen. 


25* 
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376 Gorup, Eineneue Methode zur Berechnung der Strömungsfunktionen „4.39 Nr. na oe 


1 a * = 
denn offenbar müssen alle Terme O ei in R+S in 1. Ordnung sich gegenseitig zerstören 
(weil A o(=)) ; für die Terme mit rein imaginärem Vorfaktor folgt dann aber dasselbe für 

0” H : 


R— S*. 
2. S*istin B,regulär, Rdesgleichen bis auf einen Pol 1. Ordnung in o=c mit dem 
reellen Residuum w. 
Wir setzen nun in (13) 
h(g) = —=8l2#21) 2 232°. Do LEE 
Dann wird 


A = Rz) Als )]r Se 


Von dieser Lösung wurde aber bereits bewiesen, daß sie eine W-Funktion ist. Aist daher die 
gesuchte Strömungsfunktion für die zirkulationsfreie Strömung im Außenraum des veränder- 
lichen Zylinders. 

Bei Verwendung der Abbildung auf 0, = E.K. r) entsteht entsprechend: 


= Re) + Sara Sr dr) 


(auch direkt aus (16) durch formale Umformung zu erhalten). Damit die rechte Seite von (17) 
eine W-Funktion wird, muß — wie man leicht findet — die Zerlegung (14) in folgender Weise 
durchgeführt werden: 


Alle Partialbrüche mit Pol in @, werden zu R genommen, alle mit Polin ZA, zu 8, außerdem 
alle Konstanten zu 8. Eine Ausnahme erfordern aber die Glieder in a Diejenigen mit re- 


ellem Vorfaktor kommen zu R, diejenigen mit imaginäremVorfaktor werden zu gleichen 
| Teilen in R und S aufgenommen. 


6. Sonderfall der starren Verschiebung und Drehung 
Bestehen die Veränderungen der Kontur C'nur aus starren Verschiebungen und Drehungen, 
dann vereinfacht sich die Berechnung der Strömungsfunktion erheblich. 


Es sei nun z(o) die Abbildungsfunktion der in der mitbewegten z-Ebene festen Kontur (©. 
Die Abbildungsfunktion bezogen auf die ruhende z,-Ebene lautet dann: 


EI ER An ine . (18). 


Dabei ist & der Drehwinkel der 2- gegen die 2,-Ebene und s der Ort des Punktes z— 0 in der 
2,-Ebene, beides zur Zeit t. 


Setzen wir (18) in die Funktionalgleichung (7) ein, dann entsteht: 
X’ — X'r— — ijherz! —lazzt + Selig’ — gez, 8) 


Beachten wir hier die Vertauschbarkeit der Symbole ’ = = und *, dann können wir beider- 


seits nach o integrieren: 
X— X*t— — 1d22* + se 2 — sei%z*. 


Dies können wir auch als zwei getrennte Funktionalgleichungen für die beiden Bewegungsfälle 
schreiben: 


X, — Kews— 2 al ee (19) 
X — RK ana Eee) DE er 


v— se”'* ist der komplexe Geschwindigkeitsvektor des Punktes z — 0 bezogen auf die z-Ebene. 
Die Funktionalgleichungen (19) und (20) ER sich in derselben Form für den Abbil- 
dungsfall C,—= E.K. Das ist auch sofort einzusehen, da diese Gleichungen keine ect 
nach der Bildvariablen o bzw. enthalten. 
Die Lösung der Funktionalgleichungen (19) und (20) geht auf einem der Lösung von 7) 


entsprechenden Weg vor sich. Nur ist zu beachten, daß X, und X, jetzt R-Funktionen werden 
müssen. 


°) Zum Unterschied von der Strömung bei veränderlicher Kontur bezeichnen wir die Strömungs- 


funktion bei starrer Bewegung mit X statt A. 
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Als triviale Lösungen — die jedoch im allgemeinen keine R-Funktionen sind — erkennt man 
in den beiden Fällen: 
v2, für (19) 


Ir, u 
Sehr 1X 22* für (20). 


Zu ihnen sind wieder geeignete Lösungen der zu (19) und (20) gehörenden homogenen Funktional- 
gleichungen hinzuzufügen. So ergibt sich 


X, = v2 — (8(v2) + [8 (v2)]*) = R(vz) — [8@2)]* 
und schließlich 
ENT SA RSN 2) 1 u ARE a EEE Er 2%): 
Entsprechend: 
1% 1% Ba + + 
= —giseert SQ are) + [sl iäzer)| = —R[oiszer) + |8(Siäzee) ; 


woraus wegen [S(z2*)]* = R(zz*) das einfache Ergebnis entsteht 
N ee he (22). 


Die beiden für (21) und (22) nötigen Zerlegungen in R+ S sind hier — für beide Abbil- 
dungsfälle — leichter durchzuführen als bei (14). Nach Partialbruchzerlegung von 2 bzw. z2* 
sind nämlich — ohne weitere formale Hinzunahme von Zusatzgliedern — KR(z) bzw. R(zz*) so- 
fort R-Funktionen, wenn nur die einzelnen Partialbrüche je nach dem Ort ihres Poles auf R und 
S aufgeteilt werden. Mit dieser Überlegung ist gleichzeitig, da offenbar auch [$(z)]*, mit den 
an (, gespiegelten Singularitäten von S, eine R-Funktion ist, der noch fehlende Beweis dafür er- 
bracht, daß die Ausdrücke X, und X, tatsächlich Rt-Funktionen, und somit die gesuchten 
Strömungsfunktionen sind. 

Eine besondere Bedeutung der Formeln (21) und (22) liegt nun auch darin, daß sie sich 
wegen der Überlagerbarkeit von Strömungen und Strömungsfunktionen gemeinsam mit der 
Formel (16) bzw. (17) verwenden lassen. Man kann sich so eine beliebige Konturänderung passend 
aus einer reinen Verschiebung (Geschw. v) und Drehung (Geschw. &) sowie aus einer Deformation 
(Geschw. 2) aufgebaut denken. Die dazu gehörenden Strömungsfunktionen werden dann einzeln 
bestimmt und addiert. — Da die Berechnung des Anteils der Strömungsfunktion infolge einer 
in der gesamten Konturänderung enthaltenen Verschiebung oder Drehung über die Formeln (21) 
und (22) erheblich einfacher ist als über (16) bzw. (17), so läßt sich durch die Zerlegung in Teil- 
strömungen der gesamte Rechenaufwand erheblich abkürzen. 

Dem Schlußergebnis kann noch eine Zirkulation beliebiger Stärke J' überlagert werden. 
Für die beiden verwendeten Abbildungsfälle gehören dazu die Strömungsfunktionen 


IT 
= In. —— bzw. ——Int. 


2nı 0—6c BEA 


7. Beispiel: Joukowski-Profil 
In der z-Ebene wird durch die Abbildungsfunktion 


2-14, — a ns 03), 


Nora 


(Areell und >0,a=p—iqmit p> 0 und g>1) aus der reellen o-Achse ein Joukowski- 
Profil erzeugt?). Die obere o-Halbebene (H,) geht, da frei von Nullstellen der Ableitung 2’, kon- 


_ formin das Äußere des Profils über. o= i liefert den Punkt 2= », o = wgibt wegen 7’ = O 0) 


die Spitze des Profils, die in z = 0 liegt. 
Sind A und a Funktionen der Zeit, dann deformiert sich das Profil, bleibt aber immer ein 


Joukowski-Profil, insbesondere bleibt die Spitze dauernd erhalten. (Darin, daß für das stän- 
dige Vorhandensein einer Spitze keine zusätzliche Bedingung mehr aufgestellt werden muß, liegt 


'%) Diese weniger bekannte Art der Abbildung der Joukowski-Profile ist hier zweckmäßiger. — Die 
Gleichwertigkeit von (23) mit der üblichen Abbildung auf einen Kreis läßt sich durch eine kurze Rechnung 
" zeigen. 

ur 


Ber DEREN A Dan en U En £ I ah a ie a Pan NEE, ER 

N TERN RE DER ALHIeR u N ee > 
gr y ee Y ar EN A TERN EL Be . ae 

Fr " SIR . Ir F r Pu: 
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der bedeutende Vorteil der Abbildung des Profils auf R.A.) Zeitliche Änderung von 4 allein gibt 
ähnliche Vergrößerung und Verkleinerung des Profils; ändert sich nur a, dann variiert die Dicke 
und Wölbung. , 

Es soll hier nur der Fall eines veränderlichen a behandelt werden. (Der andere Fall verläuft 
in der Rechnung ähnlich, das Ergebnis ist zu dem hier hergeleiteten zu addieren.) . 


Aus (23) bildet man 


i AA 1 
a, Rare 
z da De: 
1 1 
+ Al ge ea. 
Der sich ergebende Ausdruck für 2*z’ wird nun in Partialbrüche zerlegt. — Das liefert: 
KE, | - 1 1 Re | 
Kol 2 EN 
| ae | (a -a®(o —a%  (a—-a)"o-a 
B 1 1 2 1 1 1 2 1 
2 7 Br Ber et 4 
ae R aloe ao i Ga ae 
1 1 2 1 | 
(a — a) (o—a)? (a—a”’o—al' 


Hier sind die einzelnen Terme, je nach dem Gebiet, in dem ihr Pol liegt, bereits aufgeteilt in zwei 
# _ pzw. un : 
—ı o—V 


zu addieren. Aus der Forderung, daß R — S* eine W-Funktion sein soll, erhalten wir dann _ 


Gruppen. Um die Anteile Rund 8 zu erhalten, ist zu den beiden Gruppen z 


BE Ba 1 
Wk 

Somit liefert (16): . 

1 1 ar a | Fat Roh 1 1 1 | 
(a— a) (o—a)? (a—a)”o—a (a—ao—i (a—-A)o+ti 

ul 1 Bl 1 re Ar 
3 Abe ner, RT ETR NEN E: (+ a)? co —a 

1 1 2 1 1 1 1 1 | 

aa lo aaa (-aori Waozil 


1-4 | 


X und bei Einsetzen von «= p— ig und Integration: 2 
43 g. ko —ajl,  2A2s Seo Eee au 1 1 
Der En o®+1 For NGFE TI Boa GEolotitoza 
Es ist nicht schwierig nun dieser Strömung infolge der Deformation des Profils eine aus Vere .) 


schiebung und Drehung resultierende zu überlagern. — Es seien v uns & die Geschwindigkeiten Bi, a 
(vgl. Abs. 6). Dann findet man direkt: RE 


— 


X, =i40- = er 


ee: 5 A Pr oti 2 5 . . i & . 

Be 5 5 SE 
= Für die Berechnung von X, hat man 2 2* zu bilden, in Partialbrüche zu zerlegen und den reg 

Anteil zu nehmen. Dies ergibt nach (22): g ae ae 


AN Fi je In 
= 2 a el a en on 7 
% Zule tale rer ur 
Es ist dann die gesamte Strömungsfunktion: 


w=A-+X, ubee 


\ 


Eingegangen am 24. Dezember 1951. 
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Über die Charakteristiken des vollständigen ebenen 
_ Plastizitätsproblems 
Von Hilda Geiringer, derzeit Rom 


Ein allgemeines kompleties Problem für ein isotropes ideal plaslisches Material wird so angesetzt, daß 
Spannungscharakteristiken und Geschwindigkeitscharakteristiken nicht zusammenfallen, und weder diese 
noch jene die Winkel der Hauptlinien halbieren. Für diesen Ansatz werden die Richtungsfelder der Charak- 
teristiken, so wie die sämtlichen Kompatibilitätsbedingungen in einfacher Weise abgeleitet; u.a. mit dem 
Ziele einer Analyse des ‚„‚speziellen‘‘ Problems, bei dem diese Kurvennetze zusammenfallen. 


A general complete problem for an isotropie ideally plastic material is pul in such a way Ihal Ihe charac- 
teristics of tension and celerily do nol coincide and neither Ihese nor those halve the angles of the main lines. 
For this purpose the fields of direction of ihe characteristics as well as all conditions of compatibility are 
derived in a simple way, among olhers with Ihe purpose of oblaining an analysis of the „special“ problem 
at which these nets of curves coincide. 


Un probleme general complet pour un materiel isotrope idealement plastique est tellement Pose que les 
caracteristiques de tension et les caracteristiques de vitesse ne coincident pas et que ni celles-ci ni celles-lü 
ne partagent en deux les angles des lignes principales. Pour ce but les champs de direction des caracteristigues 
ainsi que toutes les conditions de compatibilite sont derives d’une manitre simple, parmi d’autres avec 
le but d’une analyse du probleme „special“, auquel ces reseaux de cowrbes coincident. 


Oömeii nmo1Hol MPoÖTeMe AH U3OTPONHOTO HNEANIBHOTO ILMIACTHYHOTO MATEPHarıa IIPHnNaeTcH 
TakoÜ BUN, YTO XapakTepnucTuku HAmpsUKeHHÜu U XapaKTepuHcTuku CKOPOCTei He COBHANAMT, 
Hu Hu Te, HU ApyTue He ABAIAOTCH ÖHcceRTpucaMM YTJIOB TIIABHBIX ımHmü. uk 9Toro Ipmema 
BBIBONATCH IIPOCTBIM CIHOCOOOM NOANA HAIPABIEHHH XAPpaKTepHuCTUK, a TAKIKe BCE YCIOBHA COB- 
MECTHUMOCTH, — M. IIP., C. Ie/IbIO aHaımaa „YACTHOH“ IIPOÖ,TeMEI, IPH KOTOpoÜ Tu ceru 
EPHBBIX COBHANAHOT. 


Einleitung 


Das fast ausschließlich studierte und in Anwendungen zugrunde gelegte vollständige 
Gleichungssystem eines ideal plastischen ebenen Körpers besteht aus sechs Gleichungen für 
die folgenden sechs Unbekannten: o,, o,, T,„—= sind die Komponenten des ebenen Spannungs- 
tensors 2; v,, v, die der Fließgeschwindigkeit und keine nicht negative Proportionalitätsfunktion, 
die von Punkt zu Punkt variieren mag. Für diese gelten 


005, ,.0€ \ 0 00, 
Er AR ET (2) 
a aa A TE On I A re a Re Bis (b) 
a 2 v, u 2 Ov, k Be =) 
EN 2 RR EST Er ee ae 
Bezeichnet man mit J die Einheitsmatrix, mit s= — p= (o, + 0,)/2 die mittlere Spannung, 
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so sind die Größen links in (ce) gleich den Komponenten des ‚reduzierten Spannungstensors 
2’—=2+n»J. Die Ableitungen, die rechts in (c) auftreten, sind die Komponenten des Tensors 


der Deformationsgeschwindigkeit E und (c) kann somit kurz in der Tensorform 


geschrieben werden. 

Die drei Tensoren, 2, 2’, und E haben in jedem Punkt dieselben Hauptrichtungen. Die 
Linien, dieinjedem Punkt die Winkel dieser Hauptrichtungenhal- 
bieren werden Gleitiinien genannt. Sie besitzen eine Reihe bemerkenswerter Eigen- 
schaften !). Im folgenden werden wir sehen, daß diese Gleitlinien in gewissem Sinne die Rollen 
von drei Kurvennetzen in sich vereinen; wir bezeichnen diese als: 1. Die 45° Linien, 2. Die 
Charakteristiken der Spannungsgleichungen, 3. die Charakteristiken der Geschwindigkeits- 
gleichungen. 

In den letzten Jahren wurde ein allgemeineres Problem behandelt, bei dem an Stelle 
von (b) eine weitgehend willkürliche Beziehung f (o,, o,, 7) = 0 tritt ?). Betrachtet man die 
Gl. (a) zusammen mit dieser allgemeinen ‚‚Fließbedingung‘, so zeigt sich zunächst, daß die 
45° Linien (so nennen wir jetzt die Linien, die die Winkel der Hauptrichtungen in jedem 


1) Siehe z.B. H. Geiringer, Fondements mathematiques de la theorie des corps plastiques 
isotropes. (M&morial des Sciences Math@m. 86) Gauthier Villars, Paris 1937. Oder R. Hill, The mathematical 
theory of plastieity. Oxford, Clarendon Press 1950. 

2, Für einige Literaturangaben siehe H. Geiringer, Proc. Nation. Acad. of Sciences, Vol. 37, 
No. 4, p. 214—220. April 1951, sowie eine im Druck befindliche Arbeit desselben Autors, in Proceedings, 
First U. S. Nation. Congress of Applied Mechanics. Chicago 1952. 


‘= 


.. daher das folgende allgemeine System an, wobei die Gl. (1) mit (a) übereinstimmen: 


EN PORN N REUTE. ETERNERT TEN VODERN ER. AL IN ELTERN Wine SOCKEN 

NW Ir $ rk N 45 r; h TANTE EAN .EE H,V o DR”. . ’ \ N ve j 

e WE % , er ERPE, “a Ra de en Bit NT, PR FU) a nt “ se 
ee Sen 
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Punkte halbieren), im allgemeinen mit den ‚„Spannungscharakteristiken“ nicht mehr überein- 


stimmen. Dies ist nur dann der Fall, wenn für f, identisch in den Argumenten, = -F Io 
® 

gilt. Ferner werden die „Spannungscharakteristiken“ mit den ie 
teristiken‘“ nur dann zusammenfallen, wenn in (c), links, an Stelle von 2” der ‚‚derivierte 
Tensor‘ von f tritt, d.i. ein symmetrischer Tensor, den wir mit Grad f bezeichnen, und der 
die Komponenten 8 ER AR besitzt. 
0032 02.11.00, 

Wir wollen in dieser Arbeit das hier angedeutete allgemeine Problem behandeln, in dem 

die genannten drei Kurvenscharen im allgemeinen gesondert auftreten, ($ 1), und die Grund- 
eigenschaften dieser Kurven feststellen. Als solche betrachten wir hier: A) Richtungsbedingungen 
der Charakteristiken ($ 2 und $ 5), B) Kompatibilitätsbedingungen entlang den Charakteristiken 
($4), C) gewisse einfache Eigenschaften der 45° Linien ($ 3). Abgesehen von dem Interesse an 
diesem allgemeinen Problem scheint es uns, daß es für das Verständnis des speziellen Problems 
(a), (b), (c) nützlich ist wenn man die Eigenschaften seiner Gleitlinien nach obigen Gesichts- 
punkten analysiert. Es wurde Gewicht darauf gelegt, die Beweise möglichst kurz und durch- 
sichtig zu gestalten, einerseits um die mathematischen Eigenschaften unseres allgemeinen 
Systems ins Licht zu setzen, anderseits um gewisse Zusammenhänge, und gewisse Unterschiede 
klarzustellen. 
Die Richtungsbedingungen der Charakteristiken folgen sachlich aus der Charakteristiken- 
theorie partieller Differentialgleichungen. Es schien uns jedoch von Interesse, die elegante 
geometrische Deutung dieser Richtungsfelder, die auf der Betrachtung der Einhüllenden ge- 
wisser Familien von „M ohrschen Kreisen‘ beruht, in die Daiskelluns: einzubeziehen. 


1. Das allgemeine Problem der ebenen Plastizitättheorie 


Wie man leicht sieht, kann man den in (c’) eingeführten Tensor & als derivierten Tensor 


der in (b) definierten quadratischen Form, q, erhalten, da z.B. n 1 = 0,, usw. daher 
1 ROFH 
De nr Grad qg. Man kann also in (a), (b), (ec) die Funktion qg der Spannungen als „‚plastisches 


Potential‘ ansehen, von dem die Beziehungen (c) abgeleitet sind. Es erscheint nützlich in der 
Aufstellung eines allgemeinen Problems, dem eine allgemeine Fließbedingung f(o,, o,, 7) = 0 
zugrundegelegt wird, den entsprechenden Ansatz wieder von einem plastischen Potential ab- 
zuleiten. Anderseits wollen wir nicht von Anfang an dieses Potential, }, mit der Fließfunktion, f, 
identifizieren. 

Da wir uns hier mit einem isotropen Material beschäftigen, müssen wir annehmen, daß 
sowohl f wie A Invarianten des Spannungstensors sind, oder, was dasselbe ist, Funktionen 
der Hauptspannungen, o, und 05; (wir wollen im folgenden o, > o, voraussetzen). Wir setzen 


00, | 0% O7 00y_ (1) 
Omi OT a y 
gr f (or, 0, )= Flo, 5,)=0 } 
Gradh= kE,..' mit «AN 2 Del A Re 
wobei 
- ‚ h(o,, Gy; r)aulo,, 05) 
und 
SRad ; Aa 
Gradh = 00, Au 
On 
2 07... 00, | 
Hier geht nun zunächst (3) durch Elimination von % in die zwei Gleichungen 
oh oh oh 


do, 90, = 3:98 ri 


Ir 


0%, 0. ,.00, 8 my 


über, wobei in üblicher Weise &, — Er Y=Y =, in 9 = Ay‘ ‚ Andergeits Ka Ko 
(1) und (2) in vielfacher Weise ein System von ZIWe i en a, Ber zwei eig one 
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so daß dann unser gesamtes System aus vier partiellen Differentialgleichungen besteht. — Be- 
merken wir noch, daß die Annahme f # h nicht notwendig eine komplizierte oder „ausgefallene“ 
Annahme darstellt. Wenn man z.B. das System (a), (c) mit irgendeiner von (b) ver- 
schiedenen Fließbedingung kombiniert, (z.B. mit einer der Bedingungen, die aus der drei- 
dimensionale Mises’schen oder aus der dreidimensionalen Saint Venantschen Fließ- 
bedingung, für 1, (0, d.h. für ‚plane stress“, folgen), so ist = g und f > 


2. Charakteristiken des kompletten Problems 


Wenn wir nun darangehen, die Charakteristiken zu bestimmen, so zeigt es sich, daß eine 
bemerkenswerte Symmetrie in der Struktur unseres Systems am besten hervortritt, wenn man 
zunächst in (1), (2) gar nichts eliminiert, sondern nur (2) nach x (oder y) differenziert. (Es ist 
leicht zu sehen, daß so eine. Differentiation nur eine triviale Charakteristik, parallel einer der 
Achsen hinzufügt, aber sonst nichts ändert.) Wir haben dann im ganzen fünf Gleichungen von 
der Form 

18 


T Ar. Au. N | | 
22 (a Ey + bir a): - 0) LTD), 


wobei die u; für o,, 0, T, d%,, v, stehen. Die entsprechenden Determinante fünfter Ordnung, die 
die charakteristischen Richtungen A: u = dx: dy liefert, ist dann 


ua —Ab;.|| = 0 


deren allgemeines Element u a;; — A b;, in der ?-ten Zeile und k-ten Spalte steht. Es ist solorl 
zu sehen, daß in dieser Determinanlte in unserem Falle an vierter und fünfter Stelle der ersten 
drei Zeilen Nullen stehen, und ebenso an erster, zweiter und dritter Stelle der beiden letzten 
Zeilen. Sie zerfällt alsoin das Produkt zweier Determinanten und es ist somit 
gezeigt, daß wir sämtliche Charakteristiken unseres kompletten Systems erhalten, wenn wir das 
System (1), (2) einerseits, das System (3) anderseits gesondert betrachten. 

Führen wir die Abkürzungen 


NE u 
ee aa 


ein, so ist die charakteristische Determinante von (1), (2) durch 
vu —ı1A 0 


erh ra) ee (4) 
nf uh u 
gegeben; also, abgesehen von y durch die Differentiation eingeführten Schar z = const, Io; 
halten wir für die Charakteristiken der SE TUE OST hungen (l), ( 


die Beziehung A?fı + Aufs TR Varlak = BEN Setzen wir A= cos y, 4 =siny, wo y den Winkel eine | 
Charakteristik mit der &-Achse bezeichnet, so wird 


Beeren (4 = . ART Re Em. F 
BR y+2 337 u een Wale Ra ie St on an. (A): 


Diese Gleichung hat zwei reelle Lösungen, wenn fi fa— 3/4 0. Bezeichnet man die charak- 
teristischen Richtungen mit m, so sieht man, daß die linke Seite von (4°) nichts anderes ist als 
die Komponente d,,„ des Tensors D= Grad f in der m-Richtung, (Komponente in der m- 


| — 
Richtung des der m-Richtung zugeordneten Vektors d,„.) Die Richtung der Spannungs- 
charakteristiken ist also durch die Gleichung 


Er ee. (5) 
definiert. Führt man die Hauptwerte des Tensors Grad fein, nämlich aus (2) = =f,,(t=1,2), 
i 0; 


und bezeichnet mit o den Winkel einer Charakteristik mit der ersten Hauptrichtung des Tensors 
Grad f, so erscheint (4’) in der Form 


F\ cos?p + F, sin®p = 0 N al 1) 
woraus 
Re ergL te i (6) 
2 
’ 
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folgt. Dies definiert zwei reelle, voneinander verschiedene Richtungen, falls F, F,< 0. Eine 
elementare direkte Ausrechnung zeigt, daß die Hauptrichtung des Tensors Grad f mit denen 
von & zusammenfallen; es gilt somit für die Winkel y, y’ der Charakteristiken mit der &-Achse, 
wenn ® den Winkel der ersten Spannungshauptrichtung mit der x-Achse bezeichnet 


vd dEn ee rweldg 
Betrachten wir nun die Geschwindigkeitsgleichungen (3) und wählen, 
‚20h 
mit En =— usw ZB 
ha &, eig; h, &y —( > hs &, an Kı Y— 0 a BE A kr (35% 
so lautet die zugehörige charakteristische Determinante 
ha ut h,? 
s en N ET N N mu A ke): 
hku+hl —hıu ET A as 
Es tritt also jetzt neben (4’) die Gleichung 
oh h 
2.008284 52 005 0sin 8 +, sind —0 EEE EHE) 


welche zwei reelle Wurzeln hat, falls h, h,— h?/A <0. Die Geschwindigkeitscharakteristiken 
sind also die Richtungen, nennen wir sie n-Richtungen, für die die Komponente @,,, — (Grad h),n 
von @ — Grad h verschwindet 


(Grad 0 ee 
und wenn wir (8) in der Form 
H,&o2 ve H, sn y—0 N. ee (a 
schreiben, (y Winkel der ersten Hauptrichtung mit einer Charakteristik), so folgt 
H 
Fan ı ERS 
an? y H, (10) 
und 
ö=d-+y, ‘=-d—y we BR (11) 


für die Winkel der Geschwindigkeitscharakteristiken mit der z-Achse. Es ist klar, daß =», 
d.h.Spannungs-undGeschwindigkeitscharakteristikenidentisch, 
dannundnurdann, wennh=f. 

Es ist sehr bemerkenswert; daß die so verschieden aussehenden Gleichungsgruppen (1), (2) 
einerseits, und (3) anderseits, die analogen Gleichungen (4) und (8) ergeben. Wir wollen aber 
noch auf einen andern Umstand als den der Struktur der betreffenden Gleichungen und 
Determinanten hinweisen: Die Gleichungen (1), (2) sind nicht linear, da die Koeffizienten in der 
differenzierten Gleichung (2) die Spannungen enthalten, d.h. die a bhängig Variabeln; 
die Koeffizienten enthalten aber nur die abhängig Variabeln und gar keine unabhängig 
Veränderlichen; dementsprechend treten in (4) nur Spannungen auf. Hingegen sind die Ge- 
schwindigkeitsgleichungen (3) linear, insofern dort die Koeffizienten keineswegs die ®,, v, 
enthalten, sondern wieder nur Spannungen (die wir nun als gegebene Funktionen der x, y be- 
trachten). Daher treten in (8) wieder nur Spannungen auf. 

Der Winkel @, (bzw. ) wird dann und nur dann gleich 45°, also orthogonales Charak- 
teristikenfeld, wenn — F\/F,—1 (bzw. — H,/H,—=1) oder, F,+F,= 0 (bzw. H,-+ B,—0). 
Wir nennen kurz eine Funktion f(o,, o,, 7)—= F(0,0,) orthogonal, wenn 

of 0] 1.200 QM 
Du ae a Nr (12) 


identisch für alle Werte der Argumente, und haben somit das Resultat: Nurfüreine 
orthogonale Fließfunktion (einorthogonales Potential)stimmen 
die Spannungscharakteristiken(dieGeschwindigkeitscharakte- 
ristiken) mit den45° Linien überein. Wenn wir im folgenden von dem ‚‚allge- 
meinen“ Problem (1), (2), (3) sprechen, wollen wir, wenn nicht anders bemerkt, voraussetzen, 
daß weder [noch h orthogonal sind. Falls beide orthogonal sind, unterscheidet sich das Problem 
nicht wesentlich von dem speziellen Problem (a), (b), (ce). 


3. Die 45° Linien. 


Für jeden symmetrischen Tensor gibt es zwei stets reelle, die Winkel der Hauptrichtungen 


halbierenden Richtungen, sagen wir «4 und », so daß bezüglich dieser die Werte in der Haupt- 


Zu ee > 


Aa a 
Be 
F air r 
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diagonale des Tensors einander gleich sind, und die in der Nebendiagonale, ihr Maximum er- 
reichen. Die ersteren sind gleich der halben Summe, die letzteren gleich der halben Differenz 
der beiden Häuptwerte, also z.B. für den Tensor E: &,,—= &,—4(& + 8,) und dies ist ge- 
mäß (c) im Falle des speziellen Problems gleich Null, während Y,, gleich dem Größtwert aller y 
im betrachteten Punkt. Diesen Eigenschaften verdanken die ‚‚Gleitlinien‘“ ihren Namen. Im 
klassischen Fall des Problems (a), (b), (ec) ist also bezüglich der 
45° Linien, die wirhier Gleitlinien nennen, die Normaldehnung 
&uu—=&öw—=0undsowohldie Schiebungy„ wie die Schubspannung 
erreichen für diese Riehtungenihren Größtwert. 

Wir sehen nun sofort, daß dies im nicht orthogonalen, allgemeinen Fall anders ist. Zunächst 
bemerken wir, daß es auch hier in jedem Punkt ein und nur ein Paar von 45° Linien gibt, da, 
wie oben erwähnt, die Hauptrichtungen für die Tensoren X und Grad A, somit auch für X und E 
übereinstimmen. Bezeichnen wir die 45° Richtungen wieder mit u, v so gilt wieder &,„= &» 
—=$(&,+ &), dies ist nun aber nicht gleich Null, da hier k(&, + &) = H, + H,#0. Daher: 
Parterdes allgemeinen Problems (1). 2),.@) haben in jedem 
Berk eedrera ursdies45 Linienvals Achsen.bezogene Schiebung y„; 
sowie Schubspannungz, ihr Maximum, während die beiden Deh- 
Dunsene,„undss„iimmer noch einander'gleich, aber nicht mehr 
alench. Nullssind. 

Wir werden sogleich sehen, daß letztere Eigenschaft nun von den 45° Linien auf die Ge- 
schwindigkeitscharakteristiken übergeht. In der Tat entnehmen wir aus (3) und (9) daß 


I RR STATE N prrea s-(13): 


Die Dehnunglängs einer Geschwindigkeitscharakteristikist Null. 
Diese Aussage gehört allerdings systematisch schon in den nächsten Abschnitt. 


4. Die Kompatibilitätsbedingungen 
In jedem Fall, auch wenn f=h, also die beiden Charakteristikennetze identisch, existieren 
zwei Gruppen von je zwei Kompatibilitätsbedingungen; wenn f=h gelten längs jeder 
Charakteristik zwei Bedingungen, deren eine sich auf die Spannungen, die andere auf die 
Geschwindigkeiten bezieht, entsprechend der Doppelwurzel A/u der charakteristischen Dete- 
minante fünfter Ordnung. 
Um die Spannungsrelationenin einer Form zu erhalten, die der des klassischen 
Falls ähnlich ist, führen wir im Anschluß an v. Mises °) als Spannungsveränderliche 9 und sein, 
wobei %, wie bisher, den Winkel zwischen erster Hauptrichtung und &-Achse bedeutet, während s 
irgend eine symmetrische Funktion von o, und o, ist, die sich zur parametrischen Darstellung der 
Fließbedingung F(o,, 0,) —= 0 eignet. Sei also 
LE SEN he ern (A) 
A R RE ER OBEN S: A 
eine Parameterdarstellung von # = (0, und führen wir, mit ee die Abkürzungen 


(06730! Oil ne 
a a p Ber ua Re U PR inne hei (15) 
o 0, 


ein, und bezeichnen schließlich mit 0/0£, 0/07 Ableitungen in Richtung der ersten und zweiten 
Hauptrichtung, so erhält man an Stelle von (1) und (2) die Gleichungen °) 
08 0% 0 9 
DE Pen m To 
Die Charakteristiken von (16) müssen natürlich mit den in (4) usw. gefundenen zusammenfallen. 
In der Tat ergibt sich aus (16) wenn dieselbe Bedeutung wie bisher hat 


16) 


ph. Ace F, ‚ 

tano=—- = = —— ER Et ey (16.); 
5 RN F, 
wie in (6). Bemerken wir noch, daß das System (16) „‚reduzibel“ ist (in einer von R.Courant 
eingeführten Bezeichnung), da es aus zwei homogenen Gleichungen besteht, die in den Koeffi- 
zienten nur die abhängigen Variabeln enthalten (in unserem Fall nur s allein). Man kann daher 
die Kompatibilitätsbedingungen auch finden, indem man aus (16) ein neues lineares Sy- 
stem ableitet, in dem s und ® die unabhängigen, &£ und n die abhängig Variabeln sind. Diese 


3) Reissner Anniversary Volume. Ann. Arbor. Mich. p. 415—429. 
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neuen Gleichungen haben Charakteristiken /, I” in der s— ® Ebene, welche die Bilder der 
Charakteristiken €, 0’ der ursprünglichen Gleichungen (16) in der x, y Ebene sind. Als Resultat 
folgt die Richtungsbeziehung in der s — Ö Ebene 
dd 1 
ds Vpg 
aus der die Kompatibilitätsbeziehungen entlang einer (, oder C’ folgen: Bezeichnet man mit 
dm, dm’ Linienelemente einer (' oder 0”, so sind diese Beziehungen 


dd 1 ds dd 1. ds 


dm Ypq dm’ dm’ Yo dm AN ee 12) 
oder mit R 
—— m —@(s pr 
En (8) 
9 — G(s) = const, Ach, 0; 89 +-G(e) = Const} entlang.On er... 2% (17°) 


Im Spezialfall der Bedingung (b) folgt mit s= (+ 0,)/2K, 5 = Ks— K, = Ks-+K, 
pqg = (0,— 0,)°/070,—=4, @(s)=4s, und somit die bekannten Beziehungen: 9 +p/?2 K= 
const, entlang der Schar der ersten, bzw. zweiten ‚‚Gleitlinien“, Beziehungen, die als Ausgangs- 
punkt der sog. Hencky-Prandtlschen Gleitliniensätze angesehen werden können. Diese 
Hencky- P randtlschen Sätze sprechen also geometrische Eigenschaften der Sp annungs- 
charakteristiken des klassischen Problems aus. 

Noch einfacher ist die Aufstellung der Kompatibilitätsbeziehungen des 
Geschwindigkeitsproblems (3). Diese sind nämlich bereits in (13) ausgesprochen. 
Aus (13) ergeben sich die entscheidenden Gleichungen für die Geschwindigkeiten: Aus 


LT cos + + 22) sin 8 008 +2 "sin? 3 


_ 84 50# sin d)cos 0 + (5 cos 545, sin ö)sin 0 


folgt, wenn d/dn Differentiation in der n-Richtung bedeutet, und die n-Richtung mit der z-Achse 
den Winkel ö, die u den Winkel ö’ einschließt 


: dv, dv, ; ax 
Enno 08.0 Sr "sinö—0, En = 7,60 08 Ö nn sind = DR ner 


eine an sich interessante est, (s. 2)) ‚oder auch, wenn wir mit K, K’ die Geschwindigkeits- 
charakteristiken in der x, y-Ebene bezeichnen 


dv,/dv, = — cot ö, entlang K; dv,/dv, = — cot RAN er Ro (19% 
wobei ö und “ durch (11) gegeben sind. Und weiters aus 


% a eh li cosö-+v», sin 0) (Fo s0+ Zrsin 8) + (cos 0 — 0,500) 


folgt Kat (18), wenn man nun mit v(oder »’) die Richtung normal zu n (oder n’) bezeichnet, 
in, vn; 
d dö d dö’ 
4 )) rl _— — ee Ina?) — Un’ =, = Pe N EEE... 20 . 
dn (7 se) v, dan Fo 0, dm (0 ) &y dn’ 0 ( ) 


Also entlang K oder K’ beziehungsweise: 
dv, TE V, dö = 0 ’ don’ AU Uyr do’ = 0 en > » (21). 


Diese wichtigen Gleichungen reduzieren sich im orthogonalen Spezialfall 4), wenn die Richtungen N 
und n’ beziehungsweise mit »’ und mit — v zusammenfallen, auf die von mir 1930 angegebenen 
Grundgleichungen des Geschwindigkeitsproblems. 
0%, i 4 
nz v0: ehr en (22) 
Wir gehen auch hier auf weitere Ausführung und Diskussion der Beziehungen (20), (21) 
nicht ein. 


4) Geiringer, H. (Pollaczek), Berichte III. Intern. Kongreß für Techn. Mich, Stockholm (1930), 
Bd. 2, 185—190. | n 
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Bemerken wir noch, daß sich zeigen läßt, daß, wenn zwei verschiedene Spannungsfelder 
entlang einer (Spannungs)-Charakteristik zusammenstoßen, die Spannungen selbst bei diesem 
Übergang sprungweise Unstetigkeiten nicht aufweisen können; dies ist nur für deren Ab- 
leitungen möglich. Hingegen ist es zulässig, daß an einer (Geschwindigkeits)-Charakteristik sich 
die tangentiale Komponente der Geschwindigkeit beim Übergang sprungweise ändert. 


5. Charakteristische Riehtungen und Mohr sche Kreise 


Wir wollen nun eine Eigenschaft der charakteristischen Richtungen zeigen, die auch zu 
deren Definition dienen könnte. Die sogenannten „Mohrschen Kreise“, die sich im ebenen 
Fall auf einen einzigen Kreis reduzieren, bilden ein bequemes und sehr suggestives Mittel zur 
Darstellung eines symmetrischen Tensors. Der Vollständigkeit halber geben wir kurz die Tat- 
sachen an. Jeden Tensor (doch denken wir an den Spannungstensor 2) ordnet bekanntlich 


jeder Richtung » einen Vektor m zu (Bild 1b), der sich in eine Normalspannung p,, = 0, in 
Richtung von v» und eine darauf senkrechte Schubspannung 7 zerlegen läßt. Wenn man nun 


jeden dieser Vektoren Dr durch einen Punkt darstellt, dessen kartesische Koordinaten, bei be- 
liebig gewähltem Anfangspunkt, o und z sind, so ist nach Mo hr der Ort all dieser Punkte für 
alle möglichen Richtungen » ein Kreis (Bild la) mit Zentrum (4 (o,—+ 0,), 0) und Radius gleich 


Bild 1 


3 (0, — 0,). Um dies zu zeigen, bezeichnen wir mit e den Winkel der v-Richtung und der ersten 
Hauptrichtung. Dann gilt bekanntlich 


o— 0,c05S?e + 0,5in?e T—= (09,— 0,)cosesin e | 
+0 H— 0 GH— 0. : DES 
Fr u ae 5 ® cos 2e = ,—sin2e | (23) 


woraus sofort die Behauptung, nämlich 


a a (Sei) Var ee 


folgt. Der Hauptwert dieser geometrischen Darstellung ist, daß sie sogleich (s. Bild 1a) zu jedem 
Punkte, d.h. zu jedem o, z, den gemäß (23) zugehörigen Winkel e liefert, und umgekehrt. In der 
Tat sehen wir, daß der in (23) eingeführte Winkel e mit dem im Bild so bezeichneten zusammen- 
fällt. Für den Winkel im Bild gilt nämlich 


T Be! 
tane= Be N (24) 
und daraus 
tanz, - 9 . (25) 
Pa 10! 


und, immer für den Winkel in (25), cos?e — nn 


2 — 2 ; man sieht aber sofort, 
| ; SEEN a 01 
daß diese Werte in (23) eingesetzt diese Gleichung identisch befriedigen. — Soweit sind dies wohl- 


bekannte Tatsachen. 


ee re + = 
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Nun betrachten wir die eindimensionale Familie vonM o hrschen Kreisen, deren Konstante, 
o, und o,, eine gegebene Gleichung @ (o,, 05) = 0 befriedigen. Diese Kreise haben eine Ein- 
hüllende, die mit jedem Kreis der Familie (mindestens) einen reellen oder imaginären Be- 
rührungspunkt gemein hat. Im Fall reellen Kontakts mit einem Kreis, Q, hat der Kontaktpunkt 
Koordinaten o, r und dazugehörige Winkel e und n. Wir wollen nun zeigen, daß, wenn wir 
G(o,, 0,) mit unserer Plastizitätsfunktion Fidentifizieren, dann 
n=%—e=gist, wo pderin (6) eingeführte Winkel der Spannungscharakteristik mit der 
ersten Hauptrichtung ist, während im Falle, daß für @ das plastische Poten- 
taäal ZHtritt,fürden entsprechenden Berüuhrungs puUakten mare 
wo y der die Richtung der Geschwindigkeitscharakteristiken definierende, in (10) eingeführte 
Winkel ist. Der folgende Beweis ist wohl auf die einfachsten Elemente reduziert: 

Benutzen wir, wie in (14), o,— 01(8), 0, 0,(s) als Parameterdarstellung von F(o,, 0,)—0, 
und gleichzeitig s als Parameter in der Darstellung unserer eindimensionalen Familie von Kreisen. 
Dann ergibt Differentiation von 


®— 0 (5 + 0) Hy + Ü= ON ne nern 2 
nach s 
(5-5 %)= 0,0% 040, =. 2 mE. Re Euro 
Diese beiden Gleichungen definieren die Enveloppe. Aus (26) entnehmen wir (0,— 0)/(o — o,), 
was nach (25) gleich tan?n ist: - 
tan de So © 
06-4 © FR, 


und somit ist, gemäß (6), unser Satz bewiesen; natürlich kann genau so bez. H (o,, 0,) geschlossen 
werden. 

Nicht reeller Berührungspunkt entspricht dem Fall nicht reeller Charakteristiken (ellip- 
tischer Fall). Betrachtet man eine allgemeine Fließbedingung, so ist es möglich, (man denke, um 
ein konkretes Beispiel zu haben, etwa an die durch = +03 — 0,9, — 4 K?—=0 definierte 
Schar von Kreisen), daß zunächst für eine stetige Aufeinanderfolge von Kreisen auf jedem Kreis 
zwei reelle Berührungspunkte liegen (entsprechend dem hyperbolischen Fall), die dann allmählich 
zusammenrücken, bis sie für einen bestimmten Kreis zusammenfallen (parabolischer Fall), und 
sodann folgen Kreise, die ganz innerhalb der Einhüllenden, ja sogar innerhalb des letzten 
„‚parabolischen‘ Kreises liegen (elliptischer Fall). (Man betrachte z. B. obiges Beispiel mit 
5—=s—V(4— 82)/3, = s+lV(4— 8)/3. Verfolgen wir die einzelnen Mohrschen Kreise 
für 0 <s <2; beginnend mit einem Kreis, dessen Mittelpunkt im Nullpunkt ist und radius 


gleich 22 ‚ rücken die Mittelpunkte nach rechts, während die Radien kleiner werden, bis, für 


y3. 
s—/|3 , der Mittelpunkt die Abscisse 4/V3, hat und der entsprechende Radius gleich 1//3 


ist; die nun folgenden Kreise, die man für V3 <s < 2 erhält, liegen ganz innerhalb dieses letzten - ) 
Kreises. Dieser letzte Kreis hat nur einen Berührungspunkt, der auf der o-Achse liegt, mit € 
der Enveloppe gemeinsam, während die Kreise mit kleinerem s je zwei, symmetrisch zur o-Achse 

gelegene, Berührungspunkte aufweisen. Für —2 <s <( ergibt sich eine zur -Achse symme- 

trische Konfiguration. 


Der obige Beweis ist unabhängig von einem von Hill (l.c. 1), p.295ff). Hill zitiert 
an dieser Stelle einen früheren Beweis von J. Mandell (1942), der mir nicht bekannt ist, 
und erwähnt in diesem Zusammenhang Äußerungen von Prandtl1 (1920). Da diese sich jeden- 
falls nur auf den speziellen Fall unserer Bedingung b), Tas, = const, beziehen, so erwähnen 
wir, daß in diesem Fall der in Rede stehende Satz selbstverständlich wird; denn die Ein- 
hüllende ist dann eine zur o-Achse parallele Gerade; der stets reelle Berührungspunkt liegt stets 
im höchsten Punkte des betreffenden Kreises, wo e= n — 45° ist, und anderseits ist in diesem 
orthogonalen Fall9=45°. In Möhrs Bruchtheorie ist der Winkel, der dem Berührungs- 
punkt mit der Hüllkurve aller ‚‚Grenzspannungskreise‘‘ entspricht, der Bruchwinkel. In einem 
interessanten Artikel hat C. Torre) gezeigt, daß für diesen Winkel Gleichung (6) gilt. Dies 
folgt für uns natürlich aus unserem Beweis, der in gleicher Weise gilt, was auch die Bedeutung 
von @ (o,, 0,) ist, ob diese eine Fließfunktion ist, oder ein plastisches Potential, oder die 
Grenzspannungskreise definiert. Für unsliegt aber das eigentliche Interessein der Iden- 
tität der in der Mohrschen Theorie ausgezeichneten Richtungen 
(definiert durch Berührungspunkt mit Hüllkurve) mit den charakteristischen 


5) Z. angew. Math. Mech. 31 (1951), S. 275ff. 
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Richtungen unseres Problems. Und wieder erscheint hier der bemerkenswerte 
Parallelismus zwischen den so verschieden gebauten Gleichungen des Spannungs- und des 
Geschwindigkeitsproblems ! 


Zusammenfassung 


In einem allgemeinen ebenen Plastizitätsproblem, bei dem Fließfunktion f und plastisches 
Potential } nicht zusammenfallen, und die Orthogonalitätsbedingung (12) weder für f noch 
für Ah zutrifft, spielen die folgenden drei voneinander verschiedenen Paare von Kurvenscharen 
eine Rolle: Die stets reellen ‚45° Linien“, sowie die (bei Zutreffen bestimmter Ungleichungen 
reellen) „‚Spannungscharakteristiken“ und ‚‚Geschwindigkeitscharakteristiken“. Die Richtungs- 
felder der beiden Charakteristikennetze sind, in voller gegenseitiger Analogie, einerseits durch (5) 
und (9), anderseits mittels der Hüllkurven gewisser M o hrscher Kreisscharen definiert. Entlang 
den 45° Linien haben Schubspannung und Schiebung ihren Größtwert; entlang den Spannungs- 
charakteristiken gelten die der Grundgleichung der Hencky Prandtlschen Beziehungen 
entsprechenden Gleichungen (17); entlang den Geschwindigkeitscharakteristiken ist die Dehnungs- 
geschwindigkeit gleich Null, eine Beziehung die zu den Hauptgleichungen für die Geschwindig- 
keiten führt. Im speziellen klassischen Problem (a), (b), (c) fallen die drei Netze in dem Netz 


der sog. ‚‚Gleitlinien‘‘ zusammen. 
Eingegangen am 14. Januar 1952. 


BUCHBESPRECHUNGEN 


Dr.-Ing. M. Strscheletzky, Hydrodynamische 
Grundlagen zur Berechnung der Schiffs- 
schrauben. Mit einem Geleitwort von Prof. 
Dr. Betz. 257 S. mit 57 Abb. Karlsruhe 1950. Verlag 
G. Braun, 

Der Verf. bringt eine Zusammenfassung der von 
ihm unter Zugrundelegung der PrandtIschen Trag- 
flächentheorie durchgeführten Untersuchung des 
dreidimensionalen Strömungsfeldes einer Schiffs- 
schraube. Dabei werden alle Einflüsse, wie der der 
Zähigkeit, der Form und Zahl der Flügel, der Strahl- 
kontraktion, der Nabenform usw. in praktisch aus- 
reichender Annäherung berücksichtigt. Darauf auf- 
bauend wird an Hand zahlreicher Berechnungsbei- 
spiele in Tabellen- und Kurvenform ein Verfahren 
zum Entwurf und zur Berechnung der in bezug auf die 
kinetischen Strahlverluste, den Profilwiderstand und 
die Kavitationssicherheit günstigsten Form der 
Schiffsschraube beschrieben, wobei nicht nur die 
Skelettform des Schraubenflügels vollständig be- 
stimmt, sondern auch die tatsächliche Form der 
Flügelprofile durch Einführung entsprechender Quell- 
Senken-Systeme berücksichtigt wird. < 

Bei der von den allgemeinen Grundlagen der Theorie 
der Schiffsschrauben ausgehenden lehrbuchartigen 
Darstellung des Stoffs wird auch den Belangen des 
Praktikers weitgehend Rechnung getragen, so daß 
diese Neuerscheinung allgemein sehr zu begrüßen ist, 
wie dies auch in dem besonderen Vorwort von A. Betz 
zum Ausdruck kommt. 


Berlin, K. Krienes. 


Prof. Dr. Oskar Becker und Prof. Dr. Jos. E. Hof- 
mann, GeschichtederMathematik. (Ge- 
schichte der Wissenschaften II. Naturwissenschaften.) 
3408. Bonn 1951. Athenäum-Verlag. Preis geb. 
10,— DM. 

Der knappe zur Verfügung stehende Raum hat die 
Verfasser, deren erster die Geschichte der antiken 
Mathematik, beginnend mit der Zeit, in der klar er- 
kennbare schriftliche Äußerungen vorhanden sind, 
deren zweiter die morgen- und abendländische 
Mathematik bis zum Beginn des zwanzigsten Jahr- 
hunderts dargestellt hat, gezwungen, eine gewisse 
Vertrautheit mit den in Frage kommenden Pro- 
blemen vorauszusetzen. Ein ausführliches, sorgfältig 
angefertigtes Register ermöglicht dem Leser eine 
weitere Verfolgung ihn interessierender Fragen. Über 


die Art wie die. Verfasser ihre Aufgabe aufgefaßt 
haben, sagt das Vorwort: ‚Die Aufgabe der modernen 
wissenschaftsgeschichtlichen Forschung auf mathe- 
matischem Gebiet muß darin bestehen, je nach der 
Eigenart der vorliegenden Probleme sowohl den 
pragmatischen wie den morphologischen Standpunkt 
zur Geltung zu bringen; den ersten aus methodisch- 
systematischen, den zweiten aus stilgeschichtlichen 
Gründen. In der sinnvollen Vereinigung beider Auf- 
fassungen, die sich nicht etwa kontradiktorisch, son- 
dern polar gegenüberstehen und demnach nur als 
Bliekpunkte von verschiedenen Standpunkten aus 
gewertet sein wollen, sehen wir einen wesentlichen 
Teil der Aufgabe, die wir uns mit der Abfassung einer 
— wenn auch nur skizzenhaften — Gesamtgeschichte 
der mathematischen Wissenschaften von den An- 
fängen bis zum Beginn des 20. Jhrh. gestellt haben.‘ 
Möge diese interessante Darstellung viele zur weiteren 
Beschäftigung mit der Geschichte der Mathematik 
anregen. 

Dresden. Willers. 

J.L. Krames, Darstellende und kine- 
matische Geometrie für Maschinen- 
bauer. Zweite erweiterte Aufl. mit 306 Abb. und 
2 Tafeln. Wien 1952. Franz Deuticke. Preis brosch. 
12,— DM; geb. 15,— DM. 

An technischen Hochschulen entwickelt sich gegen- 
wärtig notgedrungener Maßen immer deutlicher die 
Tendenz, das Studium der mathematischen und 
naturwissenschaftlichen Grundlagenfächer hinsicht- 
lich des Stundenausmaßes stark zurückzudrängen. 
Die Vertreter dieser Fächer sehen sich demnach ge- 
zwungen,ihren Vortragsstoff einschneidend zu kürzen. 
Die Wahrung eines hochschulmäßigen Niveaus ist 
dabei meist nur dadurch möglich, daß man die Vor- 
lesungen, getrennt nach den verschiedenen technischen 
Fachrichtungen, deren besonderen Erfordernissen 
tunlichst anpaßt. Für das Gebiet der Darstellenden 
Geometrie, die bei der genannten Entwicklung be- 
sonders stark in Mitleidenschaft gezogen wurde, stellt 
das vorliegende Werk den sehr erfolgreichen Versuch 
dar, den für die Studierenden und Ingenieure des 
Maschinenbaues erforderlichen Stoff in einem Band 
von mäßigem Umfang von Grund auf in leicht ver- 
ständlich und doch wissenschaftlich gediegener und 
strenger Form aufzubauen. Es war dabei ein besonders 
glücklicher Gedanke, die Grundlagen der kinematischen 
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Geometrie, die heute an fast allen Technischen Hoch- 
schulen von dem Vertreter der Darstellenden Geo- 
metrie für Maschinenbau vorgetragen werden, mit 
einzubeziehen. Dem Verfasser, einem hervorragenden 
Vertreter der auf kinematischem Gebiet bekannten 
Wiener Schule, kamen dabei seine langjährigen 
Unterrichtserfahrungen an Technischen Hochschulen 
zu gute. 

Der Zweckbestimmung für den Maschinenbau ent- 
sprechend steht inhaltlich die Behandlung der krum- 
men Flächen in zugeordneten Normalrissen im Vorder- 
grund, während z.B. die Perspektive nur gelegent- 
lich gestreift und die Gesetze der axonometrischen Ab- 
bildung nur etwa in jenem Umfang entwickelt werden, 
wie sie zur Ausbildung der für Maschinenkonstruk- 
teure wichtigen Fertigkeit im axonometrischen 
Skizzieren gebraucht werden. Die in dem Band auf- 
genommene Grundlage der kinematischen Geometrie 
der Ebene und des Raumes mit ihren Anwendungen 


auf die Aufgaben der Verzahnungstheorie setzen den 


Leser in den Stand, die bewegungsgeometrischen ° 


Eigenschaften der im Maschinenbau vorkommenden 
Getriebe zu erkennen und konstruktiv bearbeiten zu 
können, 

In der vorliegenden zweiten Auflage waren nur ge- 
ringfügige Verbesserungen und Ergänzungen gegen- 
über der ersten Auflage erforderlich. Dagegen ist 
erfreulicherweise eine Sammlung zahlreicher Auf- 
gaben zum gründlichen Einüben des dargebotenen 
Lehrstoffes hinzugekommen, wobei vorwiegend ma- 
schinentechnische Objekte herangezogen wurden. 

Nach der günstigen Aufnahme der ersten Auflage, 
bei der gediegenen Ausstattung, die die zweite Auflage 
durch den Verlag erhielt und bei dem mäßigen Preise 
ist zu hoffen, daß das Werk auch weiterhin die Ver- 
breitung findet, die es verdient. 


Freiberg (Sa.) W. Sehmid. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 
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Technische Hausmitteilungen des Nordwestdeutschen 
Rundfunks. Sonderheft: Unterlagenfür UKW- 
Netzplanungen. 448. mit 39 Abb. Hamburg 
1952. H. H. Nölke Verlag. Preis brosch. 6,50 DM. 


Paul Günther, Zur Gültigkeit desHuy- 
gensschen Prinzips bei partiellen 
Differentialgleichungen vom nor- 
malenhyperbolischen Typus (Berichte 
über die Verhandlungen der Sächsischen Akademie 
der Wissenschaften zu Leipzig. Mathematisch-natur- 
wissenschaftliche Klasse, Band 100. Heft 2). 43 8. 
Berlin 1952. Akademie-Verlag. Preis brosch. 5,— DM. 


Alfred Müller, DieSchaubarkeitin der 
Axonometrie (Berichte über die Verhandlungen 
der Sächsischen Akademie der Wissenschaften zu 
Leipzig. Mathematisch-naturwissenschaftliche Klas- 
se. Band 100, Heft 3). 22 S. mit 3 Tafeln. Berlin 
1952. Akademie-Verlag. Preis brosch. 3,— DM. 


Prof.Dr.H. Brandt, Über Stammfaktoren 
beiternären quadratischen Formen 
(Berichte über die Verhandlungen der Sächsischen 
Akademie der Wissenschaften zu Leipzig. Mathe- 
matisch-naturwissenschaftliche Klasse, Band 100. 
Heft 1). 248. Berlin 1952. Akademie-Verlag. Preis 
brosch. 2,25 DM. 


Dr. Milan Vidmar (Prof. a. d. Techn. Hochschule 
Ljubljana), Die Gestalt derelektrischen 
Freileitung. (Lehr- und Handbücher der 
Ingenieurwissenschaften. Band 21). 199 S. mit 49 Abb. 
Basel 1952. Verlag Birkhäuser. Preis brosch. 16,65 Fr. 
geb. 19,75 Fr. - 


Dr. W. Meyer zur Capellen, Instrumentelle 
MathematikfürdenIngenieur. (Fach- 
bücher für Ingenieure.) 383 S. mit 190 Abb, Essen 
1952. Verlag W. Girardet. Preis geb. 27,80 DM. 


NACHRICHTEN 


Prof. Dr. Pöschl 70 Jahre alt 


Am 6. September vollendete der o. Professor der 
Mechanik und der angewandten Mathematik an der 
Technischen Hochschule Karlsruhe Theodor Pöschl 
das siebenzigste Lebensjahr. Er studierte an der 
Technischen Hochschule seiner Geburtsstadt Graz, 
an der er auch 1907 promovierte und sich 1910 bei 
Wittenbauer für ‚Allgemeine und technische Me- 
cehanik“ habilitierte. 1912 wurde er an die Deutsche 
Technische Hochschule Prag zunächst als a. o. Pro- 
fessor berufen, wurde aber schon 1916 zum Ordinarius 
ernannt. Wiederholt war er in Prag Dekan und Rek- 
tor. Dort schrieb er außer zahlreichen, wertvollen 
wissenschaftlichen Abhandlungen, die zum größten 
Teil in der Z. angew. Math. Mech. veröffentlicht sind, 
eine Reihe bekannter und viel gebrauchter Lehr- 
bücher, so das zweibändige Lehrbuch über Technische 
Mechanik, zusammen mit Terzaghi das Buch über 
Berechnung von Behältern, zusammen mit Witten- 
bauer das über graphische Dynamik; ferner über- 
arbeitete er die Wittenbauersche Aufgabensammlung 
zur Technischen Mechanik. Fast .alle diese Bücher 


erlebten mehrfache Auflagen. Zu erwähnen ist weiter 


seine Mitarbeit an dem Handbuch der Physik und an 
dem Lehrbuch der Physik von Müller-Pouillet. 


1928 folgte Pöschl einem Ruf an die Technische 
Hochschule Karlsruhe. Dort wurde ihm von 1937 bis 
1945 durch die nationalsozialistische Regierung die 
Lehrtätigkeit verboten. Diese Zeit benutzte Pöschl 
zu eingehenden Forschungen auf dem Gebiet der 
Werkstoffprüfung und der Plastizitätstheorie, aus 
denen u.a. die in Zusammenarbeit mit der Firma Leitz 
gebaute Mikrozerreißmaschine hervorging. 1945 über- 
nahm er wieder seinen Lehrstuhl und wurde für das 
Jahr 1946/47 zum Rektor seiner Hochschule ge- 
wählt. In diesem Jahr fand die von ihm vorbereitete 
Tagung der GaMM in Karlsruhe statt. 1949 ließ er 
eine ausgezeichnete Einführung in die analytische 
Mechanik erscheinen. Möge ihm, der kurz vor seinem 
siebenzigsten Geburtstag noch an der internationalen 
Mechanik-Tagung in Istambul teilnehmen konnte 
die bisherige geistige und körperliche Frische und die 
Freude an wissenschaftlicher Arbeit lange erhalten 
bleiben. Willers. 
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